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responsible for the Flight Test Plan". 

Fourth line, change from "meansur ators" to "mensurator s" . 

, (3) 

For Figure 5-43 change caption from "(Camera Looking in Y 
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For Flight 16, under Stilwell Parameters, add new data in column 
to the right of 8, 1500 ft, 1/200 af 22, new data should read: 

No. of Runs 4 

Altitude 10, 000 ft. 

Camera Parameters 1/200 at f 22 

Sixth line, change from "P " to "P " • 
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FOREWORD 


This report contains the results of the Space Geodesy Aircraft Experi- 
ment awarded Raytheon Company under Contract No. NASW-1932 by the Geodetic 
Satellite Program Office, Office of Space Science and Applications, National 
Aeronautics and Space Administration. 

The experiment was conducted by the Equipment Division of Raytheon 
Company under the direction of Mr. Myer Kolker as Program Manager with 
Mr. Emile Genest as Project Engineer and Dr. Charles J. Mundo, Jr. , 
as Program Scientist. • 

Successful implementation of this effort was due largely to Mr. Jerome 
D. Rosenberg, Manager, Geodetic Satellite Programs, NASA OSSA, who ■ 
initiated the effort, and provided the necessary initial direction and guidance, 
and to Mr. Allen Selser, Technical Monitor, NASA Wallops Station, who continued 
the direction and guidance through the conclusion of the experiment. 

The primary objective of this experiment was to obtain measurements 
at normal incidence of radar backscatter and waveforms for varied ocean 
conditions. 
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ABSTRACT 


Measurements of the radar pulse shape and cross section per unit area, 
cj°, at vertical incidence from various ocean conditions were made during 
the 1969-1970 Winter approximately 120 miles east of Norfolk, Va. The 
radar equipment consisted of an X-band transmitter, receiver and antenna sys- 
tem generating pulses of ten through one hundred nanoseconds. The reflections 
were received on a high speed oscilloscope inside an aircraft flying at 10, 000 
feet. Ocean truth was provided: (1) by two cameras located in the aircraft, 

one of which obtained pictures which were subsequently processed to provide 
two-dimensional Fourier transforms of the ocean surface; (2) by a NASA 
ship on location which provided measurements of ocean and atmospheric con- 
ditions, and; (3) by a second aircraft with a laser profilometer which pro- 

o 

vided precise measurements of the ocean waves* The results indicate a a 
ranging from 8 to 21 dB and a noticeable trend for a versus wind speed* The 
waveform of the return compared favorably with the* expected waveform* 
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SPACE GEODESY AIRCRAFT EXPERIMENT 
EXECUTIVE SUMMARY 

During December, 1969 and January, 1970 sixteen test flights were 
flown out of NASA's Wallops Station, Virginia with the purpose of measuring 
radar pulse shapes and radar cross sections per unit area, a°, at vertical 
incidence from various ocean conditions. 

The radar equipment consisted of an X-band (9 GHz) transmitter, re- 
ceiver, and antennas (standard gain horns) capable of generating pulses as 
narrow as ten nanoseconds. The transmitter was a low level local oscillator 
with a TWT amplifier (12 watts) after the pulse modulation. There were two 
receiving antennas (standard gain horns of 22 dB with a 12° beamwidth) that 
were cross polarized with respect to each other permitting measurements 
of both direct and cross polarized return energy. The transmitting antenna 
(8° beamwidth) was mechanically bo re sighted to the same axis as the receiv- 
ing antennas. The receiver was a fixed gain superheterodyne type. A wide 
band oscilloscope was used to indicate the sea echos which were photograph- 
ically recorded. The number of pulses that were integrated, displayed and 
photographed varied from one to 278. Most of the flights were flown at 
10, 000 feet using a 20 nanosecond pulse width which resulted in pulse limited 
returns (rather than beam limited). 

Ocean truth data was provided by (1) a down-looking camera which 
photographed the ocean corresponding to the radar scope photographs, (2) 
a second camera which obtained photos which were subsequently processed 
to provide two-dimensional Fourier transforms of the ocean surface, (Stilwell 
process), (3) a NASA ship on location at the test area which supplied measure- 
ments of ocean and atmospheric condition, and (4) an ONR aircraft with a laser 
profilometer which yielded precise measurements of the ocean waves. 

Approximately 600 frames from ten of the flights were selected for 
further examination. A point coordinate mensurator was used to precisely 
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measure the return pulse peak and noise level, and digitally record these 
points on to punched cards. All the pertinent radar and environmental data 
were recorded onto cards. These cards were used as the input data for a 
processing program which calculated, by using a form of the radar equation, 
the radar cross section, cr°, for each frame. For each flight and for a 
selected number of runs, the mean, the standard deviation, and frequency 
distribution were then calculated. 

The spread of values of cf° ranged from approximately 8 dB to 2 1 dB 
for wave heights up to 12 ft and wind speeds up to 26 knots. In comparing 
<r° for each flight with the sea condition, there appears to be correlation 
between o-° and wind speed. Since the capillary waves are directly related 
to wind speed, the data indicates that at X-band o-° is probably more dependent 
on capillary wave formations than any other ocean parameter. The data also 
showed no significant change in or° when the aircraft was flying at various head- 
ings with the sea direction. Expected pulse shapes, calculated for various 
pulsewidths, altitudes and cr°, compared favorably with the actual pulse shapes. 
The rise times were linear and equal to the pulsewidth and the decay times 
were related to altitude and beamwidth. 

It is recommended that the analysis be further extended to measurements • 
of rise and decay times along with amplitude distributions as a function of time. 
The Stilwell process should be further developed by analyzing all the available 
photographs. The relation between cr° and capillary waves should be further 
investigated by implementing a method of measuring the capillaries and per- 
forming further aircraft experiments. These aircraft experiments should 
also include measurements of. correlation between pulses, altitude bias from 
wave heights, waveform sampling techniques, and candidate data processors. 



SECTION 1. INTRODUCTION & SUMMARY 


The Space Geodesy Altimetry Aircraft Experiment Program supports the 
goal of developing a satellite radar altimeter which will be capable of measur- 
ing the distance between the orbit and the ocean surface to accuracies of ± 1 
meter or better. At the present time, there is a lack of empirical data on 
the way in which radar pulses are backscattered by the ocean at vertical in- 
cidence. 

During this Aircraft Experiment empirical data was taken which will aid in 
developing a valid electromagnetic model of the ocean at X-Band as a function of 
pulse-length, sea-state, and altitude. Radar returns from a variety of ocean 
conditions were observed, recorded, and analyzed. This data will be useful for 
designing a satellite radar altimeter capable of measuring heights above the ocean 
to accuracies of ± 1 meter or better at altitudes of 1000 to 1500 kilometers and at 
orbital velocities of 7 to 8 kilometers per second. In addition, the data will also 
be useful in interpreting the data from the GEOS C Radar Altimeter Experiment 
in which accuracies of ± 5 meters are expected and in which consideration is being 
given to recovering the shape of the return waveform. 

The objectives of the flight tests were to measure the pulse shape as 
reflected from different ocean conditions and to attempt to find a correla- 
tion between the return pulse shape and the ocean parameters. Specifically, 
the intent was to obtain a measure of cr° (radar scattering coefficient) and 
H(<d) (impulse response) for various states of the sea or ocean. 
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This final report describes the equipment used for making the measure- 
ments (see Figure 1-1). It describes both the radar equipment and the auxil- 
iary equipment for making the necessary measurements of ocean and aircraft 
characteristics. Briefly stated, the basic radar system consisted of three 
antennas, a receiver rack and a transmitter rack. The antennas were mech- 
ically installed, and bo’resighted in the NASA C54 aircraft. The receiver and 
transmitter units were mounted in standard nineteen inch racks in the aircraft. 
Three operators were required to operate the radar equipment, auxiliary data 
recording equipment, and the three cameras (oscilloscope, ground viewing, 
and ocean spectra). 

The data gathered in this experiment consisted of approximately 10, 000 
frames of oscilloscope photographs, 10, 000 frames of ocean photographs at 
vertical incidence, and 150 frames of ocean photographs at 45° incidence. 600 
frames of oscilloscope photographs were reduced in order to obtain <r° and pulse 
shape measurements. The process used in reducing the 600 frames of data is 
described in Section 4. 0. 

This report also attempts to relate cr° to wave height, wave period, 
wind speed and roll angle. The a° measurements are made with narrow pulse 
returns as would be encountered in a satellite radar altimeter. The rela- 
tions between pulse shape, altitude, pulse width and wind speed are derived 
and compared to actual measurements. Various methods of obtaining ocean 
parameters are described and compared. These include "eyeball" estimates, 
laser profilometer measurements, and two dimensional Fourier transforms 
using photographs. The "eyeball" estimates of wave height, and wave period and 
anemometer readings of wind speed were obtained from the NASA Range Recoverer 
ship (Figure 1-2} on location in the measurement area. The laser profilometer 
measurements were obtained by ONR personnel who used a Cessna 310 aircraft 
with the laser equipment aboard (see Figure 1-3} in the same area. The ocean 
spectra photographs were obtained from photos taken by a side looking camera on 
the NASA C54 aircraft used in the experiment (see Figure 1-4, 1-5). 
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Figure 1-3. CESSNA 310 Aircraft With Laser Equipment Aboard 
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Figure 1-4. Ocean Spectra Camera (Internal) 
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Summarizing the conclusions and recommendations: 


Conclusions 

o 

1. o- appears to have a functional relation to wind speed but none 
to wave height, period, or direction. 

2. - Variations in <r° over the looking angles of 0 to 4 degrees are 

small. 

3. The distribution of the amplitude of the return signal and hence cr°, 
appears to be Rayleigh. 

4. Cross polarized returns at normal incidence were 30 dB down 
from direct polarized returns. 

5. The return waveform and cr° are unrelated to sea direction. 

6. The leading edge of the return waveform was a consistent ramp 
with little fluctuations. 

7. The trailing edge of the return waveform contained large fluctua- 
tions for the individual pulse but the average value had a predictable 
decay functionally related to altitude, pulsewidth and beamwidth. 

8. The average characteristics of the measured return waveforms 
agreed with the computed waveforms. 

9. Ocean parameters were consistently obtained by visual observa- 
tions. 

10. Implemented techniques for measurement of ocean parameters 

included Stilwell photography and laser profilometry. The former 
is simpler but it requires further development. The latter is 
more complex and provides one dimension spectrum. 
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Recommendations 


1. The same analyses conducted on the chosen 6 00 frames of data be 
extended to the remaining frames to fully exploit the data from 
this experiment. 

2 . Measure rise times and decay times and correlate with theoretical 
values. 

■t 

3. Analyze the distribution of pulse amplitude as a function of time 
using a densitometer, 

4. The Stilwell process should be further developed by processing all 
the available photographs and extending the analysis of the results. 
This includes, at least, analyzing frames from the same flight for 
self consistency, correlation of Stilwell and laser profilometry 
data from flight 5, and correlation with surface M eye ball 11 reports. 

5. The relationship of c r to surface wind should be further investigated 
including flying over areas where the surface has no capillary waves 
(no wind) present and includes conditions of no waves present and 
swell waves present, 

6. Altitude biases as a ftmction of sea state should be investigated. 

This includes ground range calibration and an independent accurate 
source of aircraft altitude. 

7. Correlation measurements between pulses as a function of time and 
frequency shift should be performed to compare with satellite 
measurements . 

8* A waveform sampler should be implemented in addition to the 
oscilloscope and camera. 

9- Candidate data processors should be flown and compared. 
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Methods for ground truth measuring of capillary waves should be 
investigated. 

In planning future aircraft experiments, major emphasis should 
be given to methods, availability, and confidence level for the 
ground truth measurement of the ocean surface conditions. 
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SECTION 2. TEST SYSTEM 


2. 1 System Parameters 

Table 2-1 gives the system parameters used in making the measure- 
ments. Much of the data was taken at 10, 000 ft with 20 nsec pulses because: 

(1) sufficient S/N was available at this altitude; (2) pulse limiting conditions 
prevailed; and (3) the higher altitudes created operational problems because 
of oxygen requirements. 

The parameter values used were dictated by various practical and 
operational considerations. The altitude of 20, 000 feet is the maximum limit 
that aircraft available for this operation can fly. The pulsewidth values are 
those which future radar altimeter designs will be considering and also those 
which can be generated with present available equipment and a minimum of 
in-house design effort. The transmitted power of 12 watts is obtained by 
using a TWT amplifier unit packaged complete with its power supply. Higher 
power units proved to be proportionately more expensive and impractical for 
this application. The antennas are standard gain horns which are readily 
available and have well controlled characteristics. The transmitted wave- 
length is approximately 3 cm or X-band. The entire X-band from 8 to 12 GHz 
was capable of being used although the experiment concentrated at one ar- 
bitrarily chosen frequency (9 GHz). X-band was used because previous studies 
have indicated that future GEOS radar altimeters will probably operate in that 
band. 

The oscilloscope camera has a shutter speed of l/50 of a second and 

the prf was determined so that one pulse or up to 278 integrated pulses can 

be displayed on an oscilloscope and photographed. This permits an evaluation 
I 

of individual pulse returns as well as the integrated effect of many pulses. 
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Table 2-1. 

Altitude (ft) 

Pulse width (nsec) 
Transmitted Power (w) 
Antenna Gain (dB) (two way) 
Frequency (GHz) 

Roll and Pitch (deg. ) 

Radar Output 
Auxiliary Data 


Pulse Repetition Rate 
Bandwidth (IF) 

Antenna Beamwidth - Transmit 
Antenna Beamwidth - Receive 
Detector 

Signal- to -Noise 
Prime Power 
Size-Receiver 
Size Transmitter 
Weight per Rack Unit 


System Parameters 

1500, 5000, 10000, 15000, 20000 
10, 20, 100 
3, 12 
47. 3 
9.0 

0, 3, 5 

Scope Photographs 

Vertical Photos 
Side Looking Photos 
Roll & Pitch Recordings 
Ship Reports on Sea & Weather 
Laser Profilometer 

0 to 5 KHz 
200 MHz 
8 degrees 
12 degrees 
Square Law 

18 dB or greater 
1200 watts 

19- inch rack 4 ft high 
19- inch rack 4 ft high 
300 lbs. 
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A wide IF bandwidth (200 MHz) was deemed desirable in order to obtain 
clear pictures of the expected fast rise times. Bandwidths wider than 200 MHz 
would have- made the system much more complex and expensive. 

The antenna beamwidths were selected so that the system remained 
pulse limited at the operating altitudes and yet provided sufficient gain. 

The detector operated as a square law device and permitted convenient 
calibration between power and voltage. 

Signal- to-noise was kept high in order to see the effect of ocean param- 
eters on the pulses independently of the thermal noise. 

The prime power, size and weight for all the necessary equipment were 
not critical parameters. 

2. 2 Radar System Components 

The equipment (see Figure 2-1, 2-2) can be broken down into three 
functional and physical subsystems. They are the transmitter, receiver, and 
antenna. The transmitter and receiver subsystems consist of rack-mounted 
units of general purpose equipments with a minimum of in-house designed 
and fabricated units. 

The antenna assembly (Figure 2-3) consists of two receiving standard 
gain horns cross polarized with respect to each other and fed by a waveguide 

switch selecting either receiving horn. A standard gain transmitting horn 

/ 

completes the antenna assembly. The E -plane vector of the transmitting horn is 
perpendicular to the line of flight. 

2. 3 Transmitter 

The type of transmitter selected for this application is a low level 
local oscillator with a TWT amplifier after the pulse modulation (see Figure 
2-4). This is preferred over the magnetron type because it permits wide 
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RECEIVE 

ANTENNA 


TRANSMIT 

ANTENNA 



Figure 2-1. Altimeter System Block" Diagram 




Figure 2-3. Antenna Assembly 



Figure 2-4. Transmitter Chassis RF Source and Modulator 



bandwidth or fast rise time pulses. It also permits variable pulsewidths and 
pulse repetition rates to be used, thereby giving flexibility to the experiment. 

The equipment shown in Figure 2-2 consists of various pieces of general 
purpose lab instruments such as the klystron power supply, pulse generator, 
TWT amplifier, and variable attenuator. The other units were designed and 
built by Raytheon. This approach of using as much off-the-shelf equipment 
was used in order to minimize the design effort and have equipment avail- 
able on a short schedule. It also has proved to be the most economical ap- 
proach and provides flexibility in parameter variations. 

The local oscillator is a low level klystron (2k25 or equivalent) which 
is readily powered with standard laboratory voltage power supplies. The 
same power supply is also used for the receiver local oscillator. The kly- 
stron is a very stable oscillator which also can be mechanically adjusted over 
the entire X-band. 

The requirement for fast rise times (5 nsed and narrow pulses (20 nsec) 
necessitates the use of a varactor switch which in turn requires an avalanche 
transistor switch driver. These are shown in Figure 2-5. The-rise time is 
controlled by the switching speed of the avalanche transistor and the varactor 
diode. These have been tested and shown to produce better than 5 nsec rise 
times . The pulsewidth is controlled by the length of the delay line used in 
the circuit of Figure 2-5. The delay lines are various lengths of eoaxial cables 
(1. 5 nsec delay/ft). 

The switch driver is triggered by a standard pulse generator such as 
those manufactured by Hewlett Packard or Data -Pulse. This will permit 
selection of pulse repetition rates up to 5 KHz and as low as one pulse (by 
manual trigger). The pulse generator is also used to provide the synch sig- 
nals for the oscilloscope and to provide the triggers for the camera controls 
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The pulses ar.e amplified in a TWT amplifier which has the necessary- 
bandwidth to preserve the pulse shape. The amplifier can provide a minimum 
of 10 w peak CW power. 

A variable attenuator (0 to 30 dB) on the output of the transmitter permits 
the use of the system at low altitudes without imposing a requirement for large 
dynamic range on the receiver. 

2. 4 Antenna and Microwave Assembly 

The transmitted and received signals are fed'to the antennas through 
..waveguide runs in. the aircraft. Waveguide is required in order to keep the 
losses to a minimum. 

A mechanical antenna switch connects either of two receiving antennas 
to the receiver. The two receiving antennas are cross polarized with respect 
to each other permitting measurements of the direct and cross polarized re- 
turn energy during flights. The receiving antennas are standard gain horns 
of 22 dB gain with a 12 degree nominal beamwidth. They have low sidelobe levels 
(less than 20 dB). Switching between antennas is remotely controlled from the 
operator's console. 

The transmitting antenna has an 8 degree nominal beamwidth. It is mechan- 
ically boresighted to the same axis as the receiving antennas. The H-plane of the 
transmitting antenna is parallel to the line of flight. 

The isolation between the transmitting and receiving antennas is re- 
quired to be in excess of 60 dB in order to avoid transmitter leakage. Tests 
were run on the antenna prior to aircraft installation (see Figure 2-6) with 
a simulation of the aircraft skin area in order to assure proper alignment, 
isolation, and beam characteristics. 

2. 5 Receiver 

The receiver used in this experiment (see Figure 2-7) is a fixed gain 
superheterodyne. type. The RF (9.00 GHz) received by the receiving antenna 
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Figure 2-6. Antenna Mockup (Front View) 


NOT REPRODUCIBLE 




is converted to an intermediate frequency of about 120 MHz in a "crystal mixer 
with no RF preselection. The local oscillator signal needed to perform this 
mixing action is obtained from a low power klystron oscillator (Yarian VA- 
203B/6975). Both this local oscillator klystron and the klystron used to pro- 
duce the transmitted signal are powered from the same power supply. 

The intermediate signal produced- by the mixer is amplified in a wide 
band IF amplifier (see Figure 2-8). 


The IF amplifier provides sufficient gain to produce an output noise 
signal after video detection of approximately 10 mV. This means that with 
a noise figure of 10 dB and a video rectification efficiency of 50 percent, the 
rms noise signal from the IF amplifier will be 20 mV (IF noise power of 8 p-W). 


The gain (RF input power to IF output power) needed to produce this 
size signal is therefore 


A 


N 


IF 

KTBF 



( 2 - 1 ) 


Gain A in dB = 60 dB. 


( 2 - 2 ) 


These gains are necessary in order to produce useful signal levels at 
the input to the oscilloscope to be used for recording the signals returned 
from the sea surface. With a S/N of 20 dB, the voltages at the scope input 
will be 100 mV and, with 20 mV/cm deflection, produce 5 cm of deflection on 
the scope face. The IF amplifier has a useful dynamic range of greater than 
40 dB. The unit requires power supplies of ± 12 volts dc. Both the IF and a 
video output are available on this unit so that other kinds of detectors could be 
used if needed (square law, coherent, etc. ). The physical size of this mixer 
IF amplifier is shown in Figure 2-8. 
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A wide-band oscilloscope (Tektronix model 585) is used to indicate the 
return from the sea. This unit has a frequency response sufficiently wide 
to obtain pulse rise times of 3 to 4 nsec. A Beattie Coleman camera was 
utilized to obtain photographic records of the resulting waveforms. Special 
films (Kodak 2485) were used to permit recording single pulse returns as well 
■as multiple returns. 

2, 6 Ground Support Equipment 

The purpose of the Ground Support Equipment (GSE) (see Figure 2-9) 
is to permit field maintenance, calibration, and sparing of the radar system. 

The sparing is accomplished by providing units in the GSE which are identical 
to the onboard equipment. The receiver and transmitter units in the GSE have 
already been described in the previous sections. 

For purposes of maintenance, a power meter, a VSWR meter, and a 
Polaroid camera have been included as part of the GSE. 

Preflight calibration of the onboard equipment is accomplished by dis- 
connecting the aircraft antennas and attaching the GSE at the coax' termina- 
tions. 

A special purpose for which the GSE was used was that of runway calibra- 
tion checks. The equipment was used as a transponder by attaching receiving 
and transmitting antennas. 

The equipment is conveniently packaged in a mobile rack unit (see 
Figure 2-10) for use and for transporting. 

2.7 Component Specifications 

Specifications for transmitter and receiver are shown in Tables 2-2- 
and 2-3. 
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Table 2-2. Transmitter Specification 


Transmitter LO: 

Low power klystron oscillator 
Power 

Tunable frequency 
Output connections - 
Convection cooled. 


20 mW 

8.5 to 9.6 GHz 
UG39/U 


Diode Switch: 

B andwidth 

Connections 

Switching Time 

Incident power 

Max I Loss 
ms 

Isolation 

TWT Amplifier: 

Gain 

B andwidth 
NF 


8.5 to 9.6 GHz 
Waveguide - UG39/U 
5 nanosec 
50 raW cw 
1.0 dB 
> 25 dB 

>30 dB 

8.5 to 9.6 GHz 

30 dB 


Max power out at gain of 30 dB 10 watts 
Integral Power Supplies rack mounted. 
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Table 2-3. Receiver Specification 


Mixer IF Amplifier: 

RF input 
RF input conn. 

IF frequency 
Bandwidth 

Noise Figure 
IF output 
RF to IF gain 
Video Output Required 

video det. efficiency > 50% 
video bandwidth 


1 dB IF compression pt 


9.0 GHz 

Waveguide (UG39/U flange) 

120 MHz 
200 MHz 

~10 dB 

Z =50 ohm 
o 

• '■■'•56 to 60 dB 


~ 100 MHz 

into reasonable resistive 
load 50 ohm — T kohm 

0 dBm 


Oscilloscope Requirement: 

1. Delayed sweep greater than 5 psec delay time 

delay jitter <0.5 nanosec’onds 

2. Rise. time < 5 nanoseconds 

3. Vertical sensitivity adjust, min. of 10 mV/cm 

4. Sweep speeds up to 5 nanoseconds/ cm. 


Camera {Oscilloscope): 

1. Solenoid activated shutter 

2. f/16 lens or better, shutter speed 1/50 

3. Capable of using at least 100 ft of 35 mm film 

4. Data storage on film 

5. Simple to load. 
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2. 8 Output Data 

The outputs of the system are pulses photographed from the face of an 
oscilloscope. A single pulse display or an integrated pulse display can be 
provided and photographed. The transmitted and received pulses can be both 
displayed on an expanded time scale and photographed. 

The oscilloscope photographs also contain a counter reading and clock 
which aid in the data reduction. 

2.9 Ground Camera 

Photographs of the ocean were taken at the same time that radar 
pictures were taken via the oscilloscope. These photos permitted a check on 
the sea conditions and more importantly helped to explain any anomalies in the 
data, such as returns from ships and clouds. 

The ground camera (Flight Research Model IV) is a 35 mm roll-film 
camera. The camera shutter is controlled by the circuits that control 
the oscilloscope camera at the operator's console. Data synchronization is 

y 

obtained through use of a counter and clock which are photographed on the film. 

2.10 Ocean Wave Measurements 

The main requirement of auxiliary data was to obtain some quantitative 
measure of the ocean- characteristics to be compared with the radar data. 

Ocean wave heights, directions, and wavelengths were the quantities of interest. 

Wind direction and wind speed indicators were furnished by the NASA ship 
stationed in the test areas. Experienced personnel aboard the ship also provided 
some "eyeball" estimates of wave height and period. 
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Photographic techniques were used to provide an optical Fourier transform 
of the photographed ocean and thereby provide the directional energy spectra of 
the ocean. This technique was developed by D. Stilwell of NRL and is referred 
in this report as the "Stilwell" or "Ocean Spectra" process. A side looking 
camera (see Figure 1-4 and 1-5) provided high resolution photographs for use in 
this technique. The camera uses nine inch film of low graininess. It is 
mounted to look out the side of the aircraft at an angle of approximately 45 
degrees from the vertical. Photographs were taken while climbing to the operating 
altitude and at the end of the run when descending. Sets of photographs were taken 
at angles of 0, 90, 180, and 270 degrees with respect to sea direction. 

2. 1 1 Aircraft Data 

Roll and pitch data was continuously recorded on a strip chart-recorder. 
Aircraft altitude, cabin temperature, outside temperature and weather were 
monitored and recorded by the operators. 
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SECTION 3. FLIGHT OPERATIONS 


3.1 Source of Data 


The data collected on these flights came from many sources, 
section lists the sources and the data supplied by each. 

3.1. 1 Raytheon 

A) Radar Data 

1. Transmitted power 

2. Antenna beamwidth, gain, and polarization 

3. Transmitted frequency 

4. Pulse width 

5. System losses 

6. Number of pulses per frame 

B) Oscilloscope Data 

1. Oscilloscope settings 

3.1.2 NASA Aircraft Data 

1. A/C Altitude 

2. A/C Heading 

3. A/ C Speed 

4. A/C Location 

3.1.3 Strip Chart Data 

1. A/C Roll Angle 

2. A/C Pitch Angle 

3. Oscilloscope Camera Synchronization 

4. Time of Day 


This 
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3. 1. 4 NASA Range Recoverer Ship 


A) Ocean Truth Data 

1. Wave direction 

2. Wave height 

3. Wave period 

4. Wind velocity (speed and direction) 

5. Air and Sea temperature 
6* Air Pressure 

3.1.5 Oscilloscope Film 

1. Received power 

2. Pulse shape 

3. 1. 6 Ocean Film (Flight Research Camera) 

1. Ocean photograph for each oscilloscope photograph 

3. 1. 7 Dr. Stilwell - Ocean Specta Film (K-17 Camera) 

1. Ocean spectra data 

3. 1. 8 ONR - Laser Profilometer 


1. Ocean spectra data 



3. 2 Summary of Flights 

The value of the flights is determined primarily from the availability of 
useful data. Table 3-1 is a summary of the radar and environmental parameters 
that were recorded for each flight. In the following paragraphs each flight is 
briefly summarized. Included also are various observations that were made 
during the flights. 

3.2.1 Flight 1 - Nov- 24- 19 6 9 

This flight, the first of three shakedown flights, was flown over coastal 
waters near Wallops. From a height of 6000 feet, 8 Stilwell photographs were 
taken and 4 radar runs were made using Flight Pattern 3, (Figure 3-3) the 
rectangular pattern. However, it was observed that sea direction is difficult 
to obtain accurately at this altitude, and that morevoer any estimate is likely 
to be off by 180°. The strip -chart was not recording for the entire flight and 
the ground camera worked only for a short while. Also, according to the 
NASA FPS-16 radar, the C-54 altimeter (SCR-718) could be off by 400-600 
feet at 6000. ft. This flight lasted two hours. 
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Table 3-1. Summary of the Radar and Environmental Parameters (Cont. ) 


RADAR PARAMETERS 

FLIGHT 4 
12/15/69 

FLIGHT 5 
12/17/69 

FLIGHT 6 
12/18/69 

FLIGHT 7 
1/6/70 

FLIGHT 8 | 

1/8/ 70 ^ 

FLIGHT 9 
1/9/70 

Pulse Width 

♦ 100 ps 

. 100 ps 

. 020 ps 








. 020 ps 





— 

NO DATA 
TWT FAILURE 












Polarization 

1 



l > t 


cross j 

direct 

cross 

direct 

direct 

direct 

direct 




direct 



direct 


50. 47 
27&, 500 

1 

2 

50 

278 

Pulse/Frame 

i 

50 

50 

50 J 


1 

50 

50 

50 

1 

Frame Rate 

NA 

1/sec 

1/sec 



NA 

l/sec 

1/ sec 

1/ sec 

NA 

NA 

NA 









* 

P. R. F. 

NA 

2500 

2500 

2500 


NA 

25CJ0 

2500 

2500 

NA 

VAR. 

NA 

NA 

2500 

1390 

Peak Power 











12 watts 

12 watts 

















Attenuation 











0 dB 

0 dB — 






* 







Sweep Speed 

. 050us/cm 









,500 

. 050 ps/cm 

. 05 ps/ ' 

cm *j 












Amplitude Set. 

. 0 50 v/ cm 

1-4 

. 200 v/ cm 
21-25 

TTO5 

v/cm 

vj cm 

• 5 v/ cm 

. 1 v/cm 

, 200 v/ cm 

. 100 

.010 

, 050? 

1* 9 ps 

. 100 v/cm 

v)cm 

. 1 

.05 

.05 

.Time Delav 

3. 1 us 

10 us 

10 us 

2, 8 us 

20 us 

3. 2 us 

36. 0 us 

36. Qus 

54. 0 ps 7 











ENVIRONMENTAL 

PARAMETERS 






Sky Condition 

Partly Clou 

An - 


Partly Clou 

‘fl V 

Partly Clou 





Cioudy 

Cloudy 




Cloudy 

.ay 










Cabin Temp, 







! 




75 *F 

70 °F - 



_ 

70 *F 



* 







* 





Water Temp, 

NA 



NA 


NA 





58 °F 

53 *F - 




6 1 *K 




* 

Wind Direction 






330 ° 





310* 

350° -J 




290* 










* 





Wind Speed 











3 knots 

19 knots 




30 knots 










* 




Sea Direction 

290 9 




320 ° 


330° 




» 

320 * 

340 * - 

— 


— ■* 

290* 

Altitude 

1500 ft 

5000 ft 

5000 ft 

1500 ft 

10,000 ft 

1500 ft 

15,000 (1 

15, 000 

20,000 ft 

1500 

15,000 ft 

no, 000' 





ft 




A/C Velocity 




155 

165 

160 

160 

160 

150 

165 

198 

gOfl . 



186 


165 






No. of Runs 

8 

12 

6 

4 

6 

4 

6 

6 

6 

4 

6 

6 

6 

6 

6 


]j;07*x73» 

37 °07' x 73*28' 


37*07 ' x 
73*38* 

Location 

7*07* x33°38' 


37*07’x73*38* 

37*07* x73*28» j 




Air Surface Temp. 






4 jor 





9. 6 *C 





-2. 3 *C 

40. 4 r 


* 


! * 






Wave Heights 






3-5 ft 





3 ft 

12 ft — 




7 ft 

0 It 














Wave Period 




_ 


NA 





5 see 





NA 

o sec ■■ 


' * j 

sec 






STILWELL PARAMETERS 
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No. of Runs 

8 

8 

8 

W 

£ 

0 

£ 

4 

4 

Altitude 

1500 

1500 1 

1500 

1000 

1000 <<t fl6 

Camera Parameters 

fl 1 Focal Length 

Ho at£U 

I 

(<? fl 1 Laser Data Available 

250 * f16 

1 and 
200 


lla 3er~ Da ta _ A va i 1 a b 1 e Laser Data Available 



Table 3-1. Summary of the Radar and Environmental Parameters (Cont. ) 


RADAR PARAMETERS 

FLIGHT 10 
1/20/70 

FLIGHT 11 
1/21/70 

FLIGHT 12 
1/22/70 

Pulse Width 

.010 [AS 

— 





• 020 us 






- tr- 

Polarization 













direct 






direct 






Pulse/ Frame 

1 

2 

10 

50 

148 


1 

2 

10 

50 

148 

278 

Frame Rate 

NA 

NA 

NA 

1 

1/2 


NA 

NA 

NA 

I 

1/2 

1/2 

PRF 

NA 

NA 

500 

o 

o 

m 

<\i 

or 

4^. 

O 

O 


NA 

NA 

500 

,2500 

5400 

13900 

Peak Power 













3 watts 





3 watts 





* 

Attenuation 






w 







C.\J UD 





< > 

0 dJb 






Sweep Speed 

. 100 

_ 

. 05 



H u 

050 







pis/cm 





< o 
p n 

pis/ cm 






Amplitude Set. 

.010/ 





w 
n o 

.05 



. 02 




cm 





£ s 

v/ cm 





- 

Time Delay 

10 |AS 




* 


21 a s 






ENVIRONMENTAL 












PARAMETERS 











Sky Condition 













wiouciy 






Clear 







70 ° IT 













iU r 






70 °F 






Water Temp. 

32 °F 






NA 















* 

Wind Direction 

ncn ° 






330° 

















Wind Speed 

3 knots 






15 knot 















"" ■*“ 

Sea Direction 

TOf a 






265“ 






Vi (JI1C 







i 

i 




Altitude 

5000 ft 






10,000 ft 

















A/C Velocity 

160 






164 knots 


— 





knots 












# of Runs 

4 

4 

4 

4 

4 


4 

4 

4 

4 

4 

4 

Location 

Tangier Sound 

i Tangier 

Long Island Sound 

Air Surface Temp. 

27°F 







NA 






Wave Heights 

Calm 

— 






2 ft 






Wave Period 

Calm 

— 





NA 






STILWELL 











PARAMETERS 










/ 

1 

# of Runs 

4 

4 


NA 

NA 


W 

Altitude 

1500 ft 

5000 ft 


NA 

NA 


z — 
o 

— i 

Camera Parameters 

250 @f8 

m@ {8 


NA 

NA 

\Z 



Table 3-1. Summary 


of the Radar and Environmental Parameters {Cont. ) 


RADAR 

•PARAMETERS 

FLIGHT 13 
1/26/70 

FLIGHT 14 
1/27/70 

FLIGHT 15 
1/28/70 

FLIGHT 16 
1/29/70 

Pulse Width 
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. 02 [is 





m 



MB 

Bl 



T 
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B 


NM 

EH 
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■be 
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wmtm 
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* 
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direct 

- 

- 


- 
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E&9H 

direct 


- 

- 

- 



- 

- 

r~ 

- 


Pulse /Frame 

1 ' 

m 

aa 

BE 



1 

B 

m 

m 

IB 


xwm 

1 

sa 

S!l 

IEE 
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HJ 

B 

M 

m 

148 

Frame Rate 

NA 


B 

E 

mi 

l 

NA 


D 

m 

□ 

D 


NA 


B 

E 

BE 

n 

NA 

E 

B 

m 

a 

1/2 


NA 

m 

m 

0! 

m 


NA 

Soo 


% v 
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NA 

B3 


g 

2 


NA 

SB 

m 



3H 
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3 watts 







m 






12 watts 





■M 
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* 

mm 


Sm 
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ti 

H l 
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Parameters 
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, 

3° 
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* 


[MB, 
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wm 
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mm 
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Si 
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Cabin Temp 

75 °F 






70 °F 
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Water Temp 

58 °F 
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m 

m 

4 ! 

Location 
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m 
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rz 
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5 sec 
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3.2.2 Flight 2 - Dec-11-1969 

This flight, the second shakedown flight, was flown over approxi- 
mately the same area as the first flight. Upon arriving at the test area a 
Stilwell pattern was flown at 1000 feet without regard to sea direction. The 
C-54 then climbed to 18, 000 feet for the radar runs. The climb took about 
45 minutes, which was considerly shorter than had been predicted. Since the 
C-54 was above 10, 000 feet, oxygen was required by the crew for the first 
time. The radar return signal was strong, but at approximately 17, 000 feet 
the oscilloscope developed synchronization problems and a trace could not be 
obtained. All radar runs, therefore, had to be aborted. A Stilwell pattern 
was flow at 18, 000 feet and then again at 1000 feet. The purpose of the second 
Stilwell pattern at 1000 feet was to observe what effect condensation on the lens 
had on lens resolution. 

3.2.3 Flight 3 - Dec -13 -196 9 

This was the last shakedown flight. Since there was no usable data 
from flight 2, a similar flight plan was used except that only the 100 nanosecond 
pulse width was attempted (It was extremely difficult to change pulse widths 
while using oxygen). The radome was removed for this and subsequent flights 
in order to improve the return signal level. Stilwell patterns were flown in 
the following order: 1000 ft. , 20, 000 ft. , 1000 ft. Some radar data was taken 
at 1000 feet and a full Radar Pattern 2 was flown at 20, 000 feet. It was difficult 
to see any return signal using cross polarization, so direct was used 
during the entire flight. The sea direction could not be determined, so that 
arbitiary headings were flown during the data runs. 

3.2.4 Flight 4 - Dec -15 -69 

This was the first data flight. It was flown over the Atlantic Ocean 
about 120 miles east of Norfolk at 37° 07* N, 73° 38' W (hereafter known as 
GEOS test area) with ocean truth data being furnished for the first time by the 
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NASA Range Recoverer ship. ONR personnel were to fly with us in their Cessna 
with the laser profilometer but they had to cancel due to aircraft problems. How- 
ever, it had been decided previously to fly at a low altitude, 5000 feet, so as to 
make coordination between the two aircraft easier. Alow cloud ceiling forced the 
cancellation of the 5000 foot Stilwell pattern. Two Stilwell patterns were flown at 
1500 feet, one before the radar runs and the other after. At 5000 feet Radar 
Pattern 1 (Figure 3-1) (12 runs) and Pattern 2 (Figure 3-2) (6 runs) were flown. 
Some single pulse per frame data was taken during the Stilwell patterns. The re- 
sults confirmed that individual pulses are visible on the film. 10 dB of attenuation 
was deliberately put in the receiver so that the return signal would not be saturated. 
The roll a nd pitch gyros were operating and the strip- chart was calibrated 
and operating. The counter of the ground (Flight Research Camera) camera 
jammed although the camera itself operated. 

3.2.5 Flight 5 - Dec-17-1970 

This flight was flown at 10, 000 feet in the GEOS test area. Due to low 
cloud cover Stilwell Patterns were only flown at 1500 feet. ONR personnel flew 
their Cessna and took laser profilometer data. The Range Recoverer was on 
station and relayed ocean truth data to the C-54. An attempt was made to change 
the pulse width but during the change a connector pin was bent, making further 
data runs impossible. Therefore, only one Radar Flight Pattern 2 (Figure 3-2) 
was flown at 10, 000 feet. 

3.2.6 Flight 6 - Dec-18-1970 

ONR personnel accompanied the C-54 again with their laser profilometer. 

The low cloud cover over the GEOS test area again forced cancellation of the 
high altitude Stilwell pattern. Radar Flight Pattern 2 was flown twice at 15, 000 
•feet and once at 20, 000 feet. There did not seem to be any noticeable return 
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Figure 3-2. Flight Pattern 2 : 



during any of the cross polarization runs, but there was a strong return 
signal when direct polarization was used. The transmitted frequency was 
measured to be 9. 00 GHz. 

3.2.7 Flight 7 - Jan-6-1970 

The Range Recoverer was on station in the GEOS test area, providing 
the only source of ocean truth data for this -flight. There was a Navy block 
to all aircraft in the area from 0 to 5000 feet with a cloud cover somewhat 
below 5000 feet. All Stilwell patterns were therefore cancelled. Also, 
because of the cloud condition, the ground camera was not used. 

Pattern 2 was flown once at 15, 000 feet with the pulses per frame varied 
during each run. Although usable data was obtained, the procedure of 
varying pulses per frame during each run was time consuming causing each 
run to cover a large area of ocean. Therefore, on following flights pulses 
per frame were kept constant during each run and only changed between runs. 

Upon returning to Wallops a system calibration check was attempted by 
flying over the GSE on the runway. The GSE picked up the aircraft signal 
for a very short time while the GSE signal was not received in the aircraft. 

3. 2. 8 Flight 8 - Jan-8-1970 

Before flying to the test area, another attempt was made at a system 
calibration check using the GSE. Four passes over the runway at 5000 feet 
gave very marginal returns but did permit a verification of the system calibra- 
tion. Because of the low ceiling and the trouble receiving the signal from the 
GSE, th e 10, 000 foot passes were cancelled. The calibration runs scheduled 
upon return from the test area were cancelled due to stormy weather. After 
arriving at the GEOS test area (37° 07' N, 73° 38'W), it was learned that the 
Range Recoverer went in toward the shore during the night because of stormy 
seas. The ship was at that time 80 miles west of the test area, at 37°N, 75° W. 
Since the ship was in coastal waters, it was decided to make the radar runs at 
the test area rather than over the ship. 
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The ship's ocean truth and weather reports were used as estimates of the 
conditions in the test area. Because of the very low overcast and surface 
fog, only one Stilwell pattern was flown at 1000 feet. Flight Pattern 2 was flown 
4 times at 10, 000 feet with the pulses per frame varied with each pattern. 

3.2.9 Flight 9 - Jan- 9-1970 

Two Stilwell patterns Were flown at 1000 ft upon arriving at the test 
area. The C-54 then climbed to 10, 000 feet for the radar runs where it was 
discovered that the TWT had failed. The transmitted and received signals 
were lost, and thus the radar runs were cancelled and the C-54 returned to 
Wallops. 

3.2.10 Flight 10 - Jan 10-1970 

The 12 watt TWT that failed was replaced with a 3 watt TWT. In order 
to get calm water conditions, this flight was flown over Tangier Sound in 
Chesapeake Bay. The Range Recoverer was in the sound and reported a calm 
surface although the surface appeared from the aircraft to have a small wave 
pattern. There was a considerable amount of ice along the-shore of the Sound, 
although the center of sound was ice free except for a few small patches. Two 
Stilwell patterns were flown, one at 1500 feet and one at 5000 feet. Radar 
Pattern 3 (Figure 3-3) was flown 5 times at 5000 feet. 


3. 2. 11 Flight 11 - Jan-21-1970 

This flight was scheduled over Tangier Sound,but during the previous 
night the Sound became completely frozen over. The Range Recoverer, still 
in the Sound, was unable to get to any other possible test areas that day. 
Therefore, the flight was cancelled and the C-54 returned to Wallops. 
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3.2,12 Flight 12 - Jan-22-1970 


This flight was scheduled over Long Island Sound in an attempt to 
get calm water conditions. Sea direction, wave heights, and wind velocity 
estimates were based on observations from the aircraft and on radio conver- 
sations with the Bridgeport, Conn. Air Control Area. There were no Stilwell 
photos taken on this flight and all the ground film was later accidently double - 
exposed during flight 13. Radar Pattern 3 was flown six times with a different 
number of pulses per frame for each pattern. The 3 watt TWT was again used. 

3.2.13 Flight 13 - Jan-27 -1970 

With the Range Recover again on station, the test flights resumed 
in the GEOS test area, 37° 07*' N, 73° 38'W. Three Stilwell Patterns were 
flown at 1500 feet, and low cloud cover forced the cancellation of the 10, 000 
foot Stilwell Pattern. As mentioned above, the ground film was accidently 
double -exposed on the film of flight 12. Radar Pattern 3 was flow six times, 
five times with direct polarization and once with cross polarization. There was 
no signal return when the horns were cross polarized. 

3.2.14 Flight 14 - Jan-27-1970 

The 12 watt TWT was repaired and reinstalled for this and the 
remainder of the flights. This flight was over the GEOS test area with the 
Range Recoverer supplying ocean truth data. Three Stilwell Patterns were 
flown six times each. For one pattern the aircraft was kept at a constant 3 roll 
angle which caused a slightly reduced return signal. The cross polarized mode was 
tried again with no success.’ Before returning to Wallops, the C-54 flew over 
Tangier Sound and found it to be ice-bound, which precluded any further testing 
in that area. The ground camera was jammed for the entire flight. 
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3.2.15 Flight 15 - Jan-28- 1970 

Radar Pattern 3 was flown six times on this flight over the GEOS test 
area. One pattern was flown with a constant 5° roll angle, which again caused a 
slightly reduced return signal. Three Stilwell patterns were flown, two at 
1500 feet and one at 10, 000 feet. 

3.2.16 Flight 16 - Jan-29-1970 

This last test flight was flown over the GEOS test area. Three Stilwell 
patterns were flown, two at 1500 feet, and one at 10, 000 feet. Radar Pattern 3 
was flown six times. On the last two patterns extended range color film (E. G. G. 
Type XR), rather than black and white film, was used. 

3. 3 Typical Data Sheets 

The data sheets presented in this section are the actual data sheets for 
flight 14. The original data sheets for the other flights are similar to the ones 
presented here. The data was recorded during the flight by the Raytheon 
controller. He received this data from the Raytheon radar operator, the 
aircraft pilot and/or co-pilot, the NASA Range Recoverer, and the ocean 
spectra (K-17) camera operator. 
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3.3.1 Flight Data Sheet 


KWH 



FLIGHT DATA SHEET ( 1 ) OF 


DATE 1/27/70 
Take-off 9:45 A.M. 
Land 1:30 P.M. 


LOCATION 

NAME OF RECORDER 

Wallops 9-D 

J. Bartlett 

STATUS OF OTHER ELECTRONIC EQUIPMENT 



EQUIPMENT STATUS ; 
PULSE WIDTH: 
POLARIZATION: 
PULSES/FRAME: 
FRAME RATE: 
BEAM WIDTH: 
PEAK POWER: 
FILM ROLL #: 


ENVIRONMENT : 

SKY CONDITION: 


(a) (20) 

(a) direct 
(see below) 
(a) 1/2 . 

(a) 8° 

(a) 10 

(a) c 


(b) 50 

(b) cross 

(b) 1 

(b) 16° 

(b) 20 

(b) BW _V 


(c) 200 (nanosecs) 

(see below) 

(c) 2 (per sec) varied 

C.<fo watts 


Uniformly Overcast 
Cloudy 

Partly Cloudy y 

Puff Ball Clouds 

Clear 

Rain 

Fog 


VISIBILITY RANGE: 


HUMIDITY: 

CABIN TEMPERATURE: 70° F 
OUTSIDE AIR TEMPERATURE: 

WIND VELOCITY AND DIRECTION FROM A/C: 


THIS SHEET PERTAINS TO RUNS - HJ - L __ 

total number OF runs for flight Run # Pulse/Frame Polarization 


1, 2, 3, 4 
5, 6, 7, 8 
9 , 10 , 11 , 12 
13, 14, 15, 16 
17, 18, 19, 20 


21, 22, 23, 24 


dir 

dir 

dir 

dir 
cross 
r 









3.3. 2 Run Data Sheet 
3900-3906 -*• Transmitted Pulse 


plight # 
14 


FRAME # 


RUN I START E 


1 3919 

2 3951 

3 3983 

4 4016 

5 4048 

6 4080 

7 4112 

8 4144 

9 4176 

10 4218 

11 4260 

12 4302 

13 4344 


Sheet 1 


1/27/70 


J. Bartlett 


AIRCRAFT 


HEADING 


VELOCITY 


ALTITUDE 


DIRECTION 

HEADING 


LOCATION 


LATITUDE LONGITUDE 


3948 

3980 

4012 

4045 

4077 

4109 

4141 

4173 

4215 

4257 

4299 

4341 

4383 




170° 

080° 

350° 

260° 

170° 

080° 

350° 

260° 

170° 

080° 

350° 


177 kts 
177 kts 
177 kts 
177 kts 
177 kts 
177 kts 
177 kts 
177 kts 
177 kts 
177 kts 
177 kts 
177 kts 
177 kts 


10, 000 ft 
10, 000 ft 
10, 000 ft 
10, 000 ft 
10, 000 ft 
10, 000 ft 
10, 000 ft 
10, 000 ft 
10, 000 ft 
10, 000 ft 
10, 000 ft 
10, 000 ft 
10, 000 ft 


350° 
350° 
350® 
350° 
350° 
350° 
350° 
350 0 
350° 
350° 
350° 
350° 
350° 


37° 07'N 
37° 07'N 
37° 07'N 
37° 07'N 
37° 07'N 
37° 07‘N 
37° 07"N 
37° 07'N 
37° 07'N 
37° 07'N 
37° 07 ''N 
37° 07'N 
37° 07'N 


73° 38'W 
73° 38'W 
73° 38'W 
73° 38'W 
73° '38'W 
73° 38»W 
73° 38'W 
73° 38'W 
73° 38'W 
73° 38'W 
73°’38'W 
73° 38'W 
73° 38'W 


00 




































3. 3. 3 Ocean. Truth Data Sheet 


FLIGHT # 

RUN (S) #(S) 

DATE 

/70 

LOCATION 

14 

AH 

1/27 

Wallops 9-D 

J. Bartlet 

t Time: 10: 

10 AM 

PARAMETER 

VALUE 

SOURCE 

SURFACE WIND 
DIRECTION 

360° 

Range Recoverer 

SPEED 

12 kts 

n 

WATER SURFACE 
TEMPERATURE 

54. 6°F 

it 

AIR SURFACE 
TEMPERATURE 

43°F 

ii 

WAVE HEIGHTS 
AVERAGE 

5 ft 

n 

(DOMINANT) 

WAVE DIRECTION 

350° 

Tt 

AIR PRESSURE 
VISIBILITY 
WAVE PERIOD 
DEW POINT 

1020. 3 mb 

8 miles ' • 

5 seconds 
36° 

11 
tl 
11 
1 1 
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3. 3. 4 Ocean Spectra Data Sheet 


OCEAN SPECTRA CAMERA (STILWELL) 


FLIGHT 
# ' 

FRAME 

# 

ALTITUDE 

A/C 

HEADING 

SEA 

DIRECTION 

TIME 

COMMENTS 


n— . 






14 

1 + 2 

1, 500 Ft. 

350° 

350° 

10:30 

Before Climb 


3+4 


260° 


10:32 

£16 @ 1 




o 



200 


5 + 6 


170 


10:33 



7 + 8 


080° 


10:34 



9 + 10 

10, 000 Ft. 

350° 


11:46 

f22 @ J_ 




o 



200 


11 + 12 


260 


1 1:47 



13 + 14 


170° 


11:48 

Broken Clouds 


15 + 16 


080° 


11:49 



17 + 18 

1, 500 Ft. 

350° 


12:02 

After Climb 


19 + 20 


260° 

i 

12:03 

f22 @ 1 


i 


0 



200 


21+22 | 


170 . 


12:04 



23 + 24 


080° 


12:05 

Heavy Haze 


TAS = 142 Kts 

TAS = True Air Speed 
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SE CTION 4 . DATA P RO CESSIN G 


4. 1 Radar Analysis 

The data from the flight tests in the form of scope film, ocean film 
(ground camera), ocean spectra film (K-17 camera), laser profilometer output, 
weather and ocean condition reports from the Range Recoverer, and strip- 
chart recordings were gathered, correlated, and analyzed with the objective of 
finding the average a 0 , (radar cross-section per unit area) for various ocean 
conditions. The o' was calculated for each of the selected frames and average 
o' 1 s were calculated for each flight. 

4. 1. 1 Procedures 

Approximately 600 frames from ten of the flights (flights 4, 5, 6, 7, 8, 

10, 12, 14, 15, 16) were selected for processing. All the pertinent environ- 

mental and radar information was recorded onto cards for computer processing. 
Using a form of the radar equation, the radar cross section, cr° , was then 
calculated for each frame, for each flight and for a selected number of runs. 

The mean, the standard deviation, the frequency distribution, and the cumulative 
probability distribution of a° were then calculated. The results are presented 
in Section 5 and Appendix A. 

4.1.2 Analysis 

The radar cross section, cr° , was calculated from the radar equation, 
using the relations derived below. 

Derivation of Equations for q° 

The radar equation can be written in the form: 
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(4-1) 


P- = 


R 


G 

X 

(T 

h 

L 


2 2 

P_ G X <r 

T 

3 4 

<4irr h L 


where 

= received power 
= transmitted power 
= peak antenna gain 
= transmitted wavelength 

- target cross section 

- altitude 

= combined system lossed 

1. waveguide loss 

2. attenuation (receiver calibration correction) 

3. antenna pattern loss 

If the return is from an area such as terrain, then 
«r = tr® A 


(4-2) 


where 


cr° = radar cross section per unit area 
A = area of terrain illuminated 


The area illuminated is calculated from the geometry of the radar 
system (see Figure 4-1). 



Figure 4-1. Geometry of the Radar System 
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The area. A, will depend on whether the radar system is pulse limited 
or beam limited. 

From the geometry, if 



c 


then the system is beam. limited, and if 



c 


then the system is pulse limited, where 

t - = transmitted pulse width 

h = aircraft altitude 

c = speed of light 

Q = incidence angle 

A. Pulse Limited Case 

The area illuminated for a pulse limited system is: 
A = ITCTh 

The relationship of the range to altitude is, 
r = h sec P sec R 

where 

P = pitch angle of aircraft 

R = roll angle of aircraft . 
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cr 


angles of the aircraft during 
less, and it can be assumed that 


(4-3) 


lystem leads to an illuminated area 


(4-4) 
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Estimation of Average <J C 


The value of the received power, P , as expressed in equations (4-3) 

K, 

and (4-4) is an average of the peak power return. For these test flights, where 

f the majority of the oscilloscope frames recorded at least 50 return pulses, the value 

of P desired for each frame is the average of the power of the pulse peaks 
& 

recorded. Since a procedure of calculating the return power for each pulse 
and then averaging these values was not feasible, an estimate of the average 
peak power was made for each frame. This was possible since each frame 
represented integration of 50 or more pulses. Figure 4-2, the oscilloscope presenta- 
tion, illustrates the maximum and average pulse peaks. Point 4 is the average 
pulse peak and point 3 is the maximum return pulse peak. As can be seen by 
this figure, point 3', can be measured rather easily, while point 4 can only be 
estimated. 

Figure 4-3 shows a single pulse representation. Point 2 is the peak 
power return. The value of P R derived from point 2 is equivalent to the P R 
derived above using the average return power for multiple pulses. 

Statistical' Groupings of g” Data 

Average <r°'s were found in various groupings of cr° . This data is 
presented in Section 5. Standard statistical equations were used to find the 
mean and standard deviation of cr 0 (as a ratio). 


of frames 
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The mean and standard deviation of cr° were later converted to dB for presen- 
tation purposes. 


The histograms of cr° are based on the mean and standard deviation with 
each bar of the histogram being one-half of a standard deviation wide. The 
number of <r° values were totaled in each of the ten ranges of the histogram 
for each grouping. The ten ranges are as. follows: (let S = ) 


o\° < o- 6 -2. OS 

1 


cr° -2. OS5 0-.' < - 1. 5S 

1 


<r° - 1. 5S < o\° < <r° - 1. OS 

l 


<r° - 1. OS < tr.° < £r° - 0. 5S 

l 


cr° - 0. 5S < <r ° < (r° 
• i 


c r° < cr° <ir H. 5S 
i 


o-° f 0. 5S < tr.° <-cr + 1. OS 

l 


o- 0 + 1. OS < o-.° < <r f 1. 5S 
i 


tr ° f 1. 5S < <r.° < o- + 2. OS 

i 


<T° f 2. 0 < 0-.° 
l 
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The cumulative probability distributions were found by integrating the 
frequency distributions, i. e. , the histograms. If f((r° ) is the frequency dis- 
tribution of cr° , then the cumulative probability distribution is: 

<r° 

F(cr° ) = ^ f(x) 

x=-oo 

4. 2 Selection of Flights and Frames 
4.2.1 Flights Selection Criteria 

Of the sixteen test flights, ten were selected for analysis. Those not 
selected included the three shakedown flights (flights 1, 2 and 3) which had 
ver.y little ocean truth data, flights 9 and 11 which were aborted before any 
radar data was taken, and flight 13 which had poor multiple pulse oscilloscope 
film quality. Frames were selected from all of the other flights in order to get 
radar data for a variety of sea conditions. 

4. 2. 2 Frame Selection Criteria 

Approximately six hundred frames were selected for processing from the 
ten selected flights. From each flight, the frames were selected in such a 
manner that the effect on the return power of varying aircraft altitude, pulse 
width, number of pulses integrated, aircraft heading, antenna polarization, 
etc. , could be studied. The frames that were eliminated from consideration 
fell into one of three categories: 

1. if the return signals were saturated or were masked by the noise, 

2. if the roll and pitch angles of the aircraft were excessive (greater 
than Si® ) at the time of the recording of the return pulse, and 

3. if there were any possibilities the film could be misinterpreted. 
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4. 2. 3 Point Measurement of the Scope Film 

In order to calculate the radar backscatter per unit area, <r° , an accurate 
value of the return power, P , was needed. The return power was measured 
from the scope film using a point coordinate mensurator, a measuring instru- 
ment with a 24" x 24" viewing screen and a magnification capability of 20 Xto 
30 X. The procedure was to precisely locate each point of interest under the 
reference crosshair of the instrument, and automatically record the coordinates 
of the point on a punched card. The measured coordinate values are in microns. 

Using Figure 4-2, the oscilloscope representation of multiple pulses as a 
reference, the following points were measured for each frame with multiple 
pulses. 

Point 1 - The top 'of the noise level at the start of the return pulse. 

Point 2 - The bottom of the noise level at the start of the return pulse. 

Point 3 - The absolute return pulse peak. 

Point 4 - The average return pulse peak. 

A. Multiple Pulse Procedure 

The value of the average return signal power, P- , was obtained in the 

R 

following manner. First, the y-coordinate of point 2 was subtracted from the 
y-coordinate of point 4 . By knowing the amplitude setting of the oscilloscope, this 
value could easily be converted to volts. However, this represents both signal 
and noise returns and to obtain just the signal return, a minimum receiver 
noise level in volts had to be subtracted from the original result. This mini- 
mum noise level was found by subtracting the y-coordinate of point- 2 from 
the y-coordinate of point 1 and then dividing this result by 2. The final value 

of P was then converted to dBm by means of a receiver calibration curve 
R 

(see Figure 5-8). 
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B. Single Pulse Procedure 


For those frames with only single pulses, the following two points were 
measured (see Figure 4-3): 

Point 1 - The start of the return pulse 

Point 2 - The peak of the return pulse 

P was fo un d for single pulse by subtracting the y- coordinate of point 1 
from the y-coordinate of point 2 and then converting this value to volts and 
finally to dBm by using the receiver calibration curve. 

4.2.4 . Card Format 

All the pertinent information concerning each frame has been stored on 
punched cards. There are three sets of cards. The first set, the flight cards, 
includes information that was constant throughtout each flight; the second set, 
the run cards, includes information that varied during the flight but was con- 
stant during each runj and the last set, the oscilloscope cards, includes the measured 
points from the oscilloscope film. 

A. Flight Cards 

The following information is included on each of the flight cards : 


1. Flight number 

2. Location code 

37° 07’N x 73° 38' W (120 mile's east of Norfolk, Ya. ) 
. Tangier Sound (in Chesapeake Bay) 

. Off shore near Wallops Island, Va. 

Long Island Sound (Middle Ground) 

3. Sea Direction (° ) 

4. Wave Height (ft) 

5. Wave Period (sec) 

6. Wind Speed (knots) 
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7. Wind Direction (° ) 

8. Water Temperature (° F) 

9. Air Pressure at sea level (mb) 

10. Peak Power Transmitted (watts) 

11. Receiver Attenuation (dB) 

12. Beamwidth (° ) 

13. Ocean Spectra Data Available (Stilwell) 

. yes 
. no 

14. Laser Profilometer Data Available 

. yes 
no 

1 5 . Sky C onditi on 

. Glear 
. Cloudy 
. Partly Cloudy 

16. Date of Flight 

B. Run Cards 

The run Data Cards included the following information 

1. ' Flight number 

2. Run number 

3. Aircraft Altitude (ft) 

4. I Oscilloscope Sweep Speed (nsec/ cm) 

5. | Oscilloscope Amplitude Setting (millivolts/ cm) 

6. Number of Pulses per Frame 

7. Transmitted Pulse Width (nsec) 

8. Antenna Polarization Code 

. Direct 


Cross 



9. Aircraft Heading (° ) 

C. 1 Oscilloscope Cards 

l The Oscilloscope Data Cards have the following data. 

1. Flight number 

2. Run number 

3. Frame number 

4. The x- and y-coordinates of the points of interest as described 
| in Section 4. 1. This is the output of the Point Coordinate 

Mensurator. 

4. 3 Computer Processing Program 

A data processing program was developed which used the flight, run, 

| and oscilloscope cards as injiut to calculate the radar cross section per unit area, 

(r° , and the return signal power, P^, for each frame. It further calculated 

for each flight the average cr° and its standard deviation, a histogram of the 

individual cr 0, s, and a cumulative probability distribution of cr° . This was 

done for the <r° based on the average peak return signal and for the <r° based 

on the absolute peak return signal (see Section 4. 1). This program (Figure 4-4) was 

written in FORTRAN IY for use on a CDC 6600. 

4.4 Results of Processing Program 

All of the output from the processing of the flight test data is presented 
in Appendix A in computer printout form. Further analysis and groupings of 
this data are found in Section 5. 
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Figure 4-4. Computer Processing Program 
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Figure 4-4. 


Computer Processing Program (Cont.) 
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A. 5 Distribution of Pulse Peaks 


The theoretical distribution of pulse peaks for a number of pulses is a Rayleigh 
distribution. This is actually the circular normal distribution in polar coordinates. 
The Rayleigh frequency distribution is given as: 



where 

a = [(x-p. ) 2 + (y-|l )] 1//2 ; the distance to the origin in a bivariate 

x y 

normal distribution 

(j ■■= a = O (x, y refers to rectangular coordinates) 
x y 

(See Figure 4-5.) The mean of the distribution equals 1.253 times the standard 
deviation. 

To show that the pulse peaks appear to have a Rayleigh distribution, the 
distributions of pulse peaks were found for a few frames from each selected test 
flight. This was done by using a microdensitiometer/isodensitracer which 
automatically scans and measures the density of points in a film transparency 
and plots the values as a quantitative two dimensional density map of the 
scanned area. The following curves (Figures 4-6 through 4-15) are the distri- 
butions of pulse peaks for one frame from each flight. The numbers on the 
vertical axis of the densitometer tracings correspond to the densities of each 
step of a "21-step wedge" (Kodak Photographic Step Tablet No. 2). The hori- 
zontal axis corresponds to the amplitude of the return signal, with the amplitude 
increasing from left to right. The point where the distribution starts to rise out 
of the noise level corresponds to the minimum return signal. 

As can be seen from the density tracings the distributions appear to be 
Rayleigh. Further analysis is required to definitely show this relation. The 
film transfer functions should be considered at that time. 
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4. 6 System Error Analysis 


The equation for <j° has been shown to be 

P (4TT) 3 h 3 L 
o R 

& - 22 

P „,G h ttct 
T 

Since the nominal roll (R) and pitch angles (P) are zero, the applicable 
error equation is 



To obtain the errors listed in Table 4-1 the calibration errors of the 
equipments (oscilloscope, power meter, etc.) were used. An additional factor 
was also included for human error in making the measurements. The major 
contributors to the error budget are the received power measurement (AP r/ R R>* 
the altitude measurement (Ah/h), and the transmitted power measurement' 

(AP j/Pj.). The power measurements involve the use of: (1) an oscilloscope 

(accuracy better than 5%); (2} power meter (accuracy better than 2%); and human 
error (estimated to be 3%) in reading the oscilloscope and/or photograph. The 
altitude measurement error (5%) is that which was specified by manufacturers 
of the aircraft altimeter. 


The error contribution (standard deviation) expressed in dB is listed in 
Table 4-1 for each error with a notation describing the sources of each error. 

The RSS error for <j° obtained from these elements is therefore 



o 


cr 


RSS « 0.9 dB 
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Table 4-1. FLIGHT TEST ERROR CONTRIBUTIONS 


ERROR 

PERCENTAGE. 

OF 

ERROR 

l<r VALUE 
(dB) 

SOURCES OF ERROR 

aP r 

(D “p^ 

P R 

10% 

0.4 

VIDEO SCOPE CALIBRATION AND 
OPERATION. INTERPRETATION OF 
POINT COORDINATOR MENSURATOR. 

(2) 

2% 

0. 1 

MEASUREMENT OF SYSTEM LOSSES 
AND IN CALIBRATION. 

(3) 3f 

3 x 5% 

0. 6 

MEASUREMENT OF A/C ALTITUDE BY 
THE-A/ C ALTIMETER. 

ap t 

(4) _i 

P T 

10% 

0.4 

MEASUREMENT OF TRANSMITTED 
POWER BEFORE FLIGHT AND SCOPE 
MONITORING OF CHANGES DURING 
FLIGHT. 

(5) 2^_ 

2 x 4% 

0. 3 

MEASUREMENT OF ANTENNA GAIN 
AT 9. 0 GHz. 

(6) 2 ~ 

2% 

0. 1 

MEASUREMENT OF FREQUENCY. 

(7) 

4% 

0.2 

MEASUREMENT OF TRANSMITTED 
PULSEWIDTH. 


RSS 


22 % 


0.9 




SECTION 5. DATA EVALUATION 


5. 1 Average Radar Cross Section (a°) 

As described in the previous section,-0 was obtained by averaging the 
multiple trace oscilloscope photographs and by computing the average of individual 
pulses displayed on the oscilloscope. The oscilloscope was used to obtain 
1, 50, 147, and 278 pulses per frame. 0° values obtained by averaging in- 

dividual pulses were about two dB less than when averages were taken of 50 or 
more pulse traces. 

The plotted values of 0° on a per flight basis (Figures 5-1, 5-2, 5-3, 5-4) 
involve the averaging of more than 18 frames per flight with each frame repre- 
senting 50 pulses. 

Plots of 0° were made for each flight in order to relate 0 values to 
ocean conditions. The numbers shown on the Figures 5-1, 5-2, and 5-3 
next to the points refer to flight numbers and hence to the day which the data 
was taken. 

The spread of values of 0^ ranges from 8 dB to 21 dB. As can be seen 
from the plots of 0° it would be questionable to claim any trend or functional 
relation to ocean parameters. There seems to be little or no relation between 
the 0° value or the change in 0 value to wave height or wave period. If there 
is any relation it exists with respect to wind speed (Figure 5-3). A trend line 
is shown on Figure 5-3 and would apply if a functional relation exists. 

Flights 8, 10, and 12 in Figures 5-1 and 5-3 are flights in which the 
relation between wind and wave height were not related as would be the case 
in a fully developed sea. Flight 8 had very high seas (greater than 12 ft) but the 
wind had started to die down. On the other hand flights 10 and 12 were made in 
inland waters with the wind blowing an insufficient time and fetch to produce a 
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Figure 5-4. Average cr° Vs Roll Angle{Pitch - 0°) 






fully developed sea. It is felt therefore that these three points on Figures 5-1 
and 5-3 reinforce the argument that o° is functionally related to wind but not to 
wave height because they improve the fit of the data in Figure 5-3 (a vs wind) 
but they destroy the fit of the data in Figure 5-1 (cr vs wave height). 

The ocean data used for making the above plots are the "eyeball" measure- 
ments obtained from the NASA ship on location when taking radar data. More 
definitive ocean parameters are possible from the Stilwell photos taken for each 
flight but the process of obtaining the necessary two dimensional Fourier trans- 
form from these photos in the Stilwell process requires further development 

(see Section 5-3). 

If a relation exists between and wind speed, it probably exists because 
of the relation between wind speed and capillary waves. If cr° is in fact related 
to capillaries (and there is much reason to believe this to be so because the cap- 
illary wavelength is comparable to the X-band wavelengths) then some methods 
for measuring capillaries will have to be devised. To date no such measurement 
capability exists. Investigations by Kinsman and Molo Christianson and cur- 
rently underway to measure capillaries. Also, Stilwell photography at low 
altitudes may make the measurements possible. 

Variations of cr° as a function of sea direction are negligible and within 
the measurement accuracy as shown in the representative measurements of 
Table 5-1, This is further substantiated in the extensive data in Appendix A 
where additional bearing angles are covered. 

Variations of cr° as a function of angle from vertical are negligible, as 
shown in Figure 5-4. This agrees in part with the curves of Moore and Schooley 
(Figure 5-5). 
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Figure 5-5. Backscatter vs Attitude 

(Reference: WHOI Report, No. 65-10, 
April 1965, Page 23) 
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Table 5-1 


Flight 

Run 

<r° vs. Sea Direction 
Bearing* 

cr° 

14 

9 

0 

12. 9 

14 

10 

90 

12. 8 

14 

11 

180 

13. 4 

14 

12 

270 

13. 0* 

* Bearing angle with 

respect to sea direction 


1. 1 Pulse 

Limited vs 

Beam Limited 



The basic radar range equation defines a (radar cross section in square 


meters) as 

follows: 


P r (4tt) 2 h 4 L 

a = 

2,2 

PGA 

t 

where 


P R 

- peak received power 

h 

= altitude 

L 

- system losses 

P 

t 

- peak transmitted power 

G 

= antenna gain 

A 

- wavelength 


The average radar cross section (a°) is the radar cross section per unit 
area or 

o a (5-2) 
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When the transmitted signal is a pulse, (see Figure 5-6) the return signal 
will be scattered from the area (A), or footprint,which is a function of time (t) or 
angle {$), where: 



c 


t is here defined from the instant the leading edge of the pulse initiates 
the first return and 6 is the angle from vertical 


The maximum value of t is t (pulsewidth),but the maximum value of Q 
is not so well defined. The 3 dB point has sometimes been used as the maximum 
value of 0, but this is purely arbitrary. 

The data in this experiment was reduced to a values by measuring F , 

P ,and h,and then performing the computation in Equation (5-1). Values of F 
t . -K- 

were obtained by measuring the peak and average voltage of the return pulse 

(see Section 4 for a more detailed description). The a values were then com- 
puted using Equation (5-2) with values of A as follows: 

A = 7T c T h • ' (5-4) 

3 

This in effect assumes an h relation for the radar range equation - a 
pulse limited relation. How valid is this assumption? 

The area (A) can be related to the angle 0 using the simple trigonometric 
relation tan Q - r/h : 


A - TT r 

2 2 

= 7T h 9 (for small angle 9) 

Note here that Q is a function of time (t) as follows: 


(5-5) 


2 

d 


= ct/h 


(5-6) 
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Figure 5-6. Geometric Relations 


If the maximum value of t (t = t) is reached before some maximum value 


of 0 (0 = 0 ), the area then becomes 

m 

A = IT c Th (5-7) 

T 

which is the pulse limited situation, or an h cubed relation in the radar range 
equation. 

If however a maximum value of 0 (0 = 0 ) is reached before t - r then 

m 

the area becomes 

A = 7T h 2 0 ^ (5-8) 

0 m 

which is the beam limited situation or an h squared relation in the radar range 
equation. 

The problem is now one of defining Q 9 or beam limiting. To define 0^ 
as the half beamwidth angle seems indefinite, and some initial calculations of 
the data with this assumption provided inconsistent values. 

The approach used was to work the problem in reverse and determine 

what 0 should be. Flight 4 provides the data for making such a determination 
m 

because on flight 4 nine runs of data were obtained with 100 nsec pulse, and nine 
runs of data were obtained with a 20 nsec pulse. Even though all runs were at 
the same altitude and over the' same ocean area>it is evident from the data that 
the received power is functionally related to the pulsewidth. 

By comparing the average received power for the two pulse widths we can 
determine 0 

m 

If 

p = average received power 100 nsec pulse 

P = average received power 20 nsec pulse 

R Z 
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From the data of flight 4 


P /P 

V *2 


= 4 


(5-9) 


Note here that if the functional relation were proportional to pulse width 

(t), then the above ratio should be five rather than four. This indicates that 

the transistion between pulse limiting and beam limiting at an altitude of 5000 

ft occurs somewhere between 20 and 100 [isec. By letting the radar range 

equations for P and P take the corresponding limiting relations 
R 1 R 2 



2 , 2 
e /h 



(5-10) 

(5-11) 


we then obtain 

9 2 = 4 cr/h. (5-12) 

m 

This gives a value of 

0=7°. (5-13) 

m 

The effective half beamwidth (3 dB) can be computed from the measured 
beamwidth by assuming a Gaussian distribution and summing the exponents. 


and 


HI 


H2 


El 


E2 


0 „ = 


0 H1 0 H2 
0 H0 +0 H2 + 


g El g E2 
0 E1 9 E2 


(5-14) 


where 0„„, 0, 0„„, 0„„ are half the 3 dB beamwidths measured from the 
HI H2 El E2 

antenna patterns and are equal to 5° , 6.5° , 4. 3° , and 6° respectively. 
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Performing this computation^ 


0 = 2.5°. (5-15) 

e 

Comparing this to 0^ = 7°, we see that the effective angle where beam 
limiting -occurs is much greater than the 3 dB angle. 

5.1.2 Maximum Value of cr° at Vertical Incidence 


If we assume that the scattering surface is aflat plate, normal to the trans- 
mission and of infinite dimensions (see Figure 5-7),and that the energy impinging 
upon the flat plate is all reflected to the source, we then can compute the maximum 
value of a°. The radar range equation for the above assumptions becomes: 

2 „ 2 
P G A 

p = - T - - . (5-16)* 

(47 r) 4h L 

The symbols are the same as previously defined. 

This is the equivalent of a transmission and reception one way over a 
distance 2 h, twice the’ "altitude. 


The comparable range equation for a reflected surface is: 

P G 2 A 2 ff°A 
T 

R (4t r) 3 h 4 L 


(5-17) 


Combining equation 5-16 and 5-17 we get 

a° = 7Th 2 /A. (5-18) 

2 2 

The maximum value of A is 7T h 0 for the beam limited condition* and 

m 

therefore 


o 


a 



(5-19) 


For the previously obtained value of 0 
a° = 18. 3 dB . 


m 


= 7°, we obtain a value 


of 


(5-20) 
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Figure 5-7. Computation of Maximum Value of a 



A similar value for the pulse limited condition is obtained,but the computa- 
tion must be based on a plate of finite -dimensions, and the two way path must be 
taken into account. 

5.1.3 Calibration 

The values of received and transmitted power used in the computation of 

< 7 ° were obtained from oscilloscope photographs. Calibrations of scope 
deflection voltage as a function of received power (Figure 5-8) were made for 
the two receivers. Calibrations of detected power level from the transmitted 
output were also made. 

The losses in the system were primarily waveguide and cable losses 
(Table 5-2). An additional loss factor associated with the antenna pattern was 

also included (Table 5-3). This pattern loss is due to the peak of the pulse 
occurring at an angle away from the beam center. In general the values associ- 
ted with the pattern loss are small enough to be negligible when considering the 
variability of the data. 

The other factor involved in calibration is the antenna gain. Antenna 
pattern measurements were made with an aircraft mockup. The results of these 
measurements were used to compute the gain. Table 5-4 is a summary of antenn 
char acteristics. 

An overall check on calibration was obtained by flying the aircraft over 
the Ground Support Equipment (GSE). In this instance the GSE received the 
signal and the value of the received power was recorded. The GSE then 
transponded the pulse and transmitted to the aircraft where it was again 
recorded. Measurements made in this fashion served as a check on the system 
gains and losses which had been independently calibrated. The operational 
requirements of this technique were not fully developed. It is felt that coordina- 
tion and synchronization of the aircraft, radar, and GSE operation must be 
automated to a considerable degree in order to improve the reliability of this 
method of calibration. 
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Table 5-2. Measured Losses 


Transmitter waveguide and coaxial losses 

a. coaxial from TWT to output coupler 1. 2 dB 

b. waveguide loss from attenuation output to 

transmitter antenna 2. 0 dB 

c. insertion loss of variable attenuation 0. 5 dB 

Receiver waveguide, waveguide switch and all coaxial, 2. 0 dB 
flexguide, etc. , in both parallel and crossed 

polarized modes 

Total Losses 5.7 dB 


Table 5-3.' Pattern Losses 

Pulse width Altitude Lo s s 

r (psec) h (ft) (dB) 

10 5000 0.3 

20 5000' 0.6 

10000 0.3 

15000 0.2 

20000 0.15 

100 5000 3.0 

10000 1.5 

15000 1.0 

20000 0.75 
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Table 5-4 


Antenna Characteristics 
£ = 9 GHz 

Transmitting Antenna - Aircraft 

Beamwidth. - H Plane 0.0° 

Beamwidth - E Plane 8.5° 

Gain 25. 2 dB 

Receiving Antenna - Aircraft 

Beamwidth - H Plane 13° 

Beamwidth - E Plane 12° 

Gain 22. 1 dB 

GSE Antennas 

Beamwidth - H Plane 18° 

Beamwidth - E Plane 18° 

Gain 18.9 dB 
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5. 2 Pulse Shapes 

5.2.1 Theoretical 

Before discussing the actual results of return pulse shapes it is necessary 
to discuss what the expected pulse should look like. Pulse shapes for various 
pulsewidths, altitudes and a ° were used to compute the expected waveforms 
(Figures 5-9, 5-10, 5-11). All waveforms .are normalized so that comparisons 
of the wave shape can be readily made. 

Figure 5-9 is based on antenna beamwidth of 6 degrees. The pulse shape 
corresponding to 100 nsec and 5, 000 ft is clearly beam limited, whereas at 
20 nsec and 20, 000 ft the pulse shape is clearly pulse limited. 

• Figure 5- 10 is for a beamwidth of 12 degrees and the distinction between 
pulse limiting and beam limiting is less pronounced. As will be seen, the actual 
pulse shape data more closely resembles the 12 degree beamwidth. 

Figure 5-11 is a plot of pulses of 20 nsec width at altitudes 20, 000 and 
10, 000 ft. The pulses are normalized and a factor is introduced which is 
based on curves by Moore and Schooley (Figure 5-5). The pulse shape changes 
very little as a function of changes in a^. 

5.2.2 Measured Pulse Shapes 

Figures 5-12 through 5-37 are actual photographs of multiple pulse 
returns under various ocean conditions, altitudes, pulse widths, etc. A 
sampling of the more than 10, 000 frames taken is used here to show some 
significant features. 

It first must be pointed out that these are reproductions of reproductions 
and at each step much information is lost. In fact, the reproductions of the 
single pulse traces were not visible at all even though they were visible on the 
original negatives. 
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Figure 5 -9 » Expected. Waveforms at BW — 6°, Wind 3-7 Kts 





Generally the frames were taken one second apart. It was clear that there 
was little variation from frame to frame. The sweep speed in all cases was 50 
nsec/ cm. Selected frames are shown' in the following figures: 

Flight 4 (Figures 5-12 through 5 -17) shows the difference between beam 
and pulse limiting. Figures 5-15 through 5-21 are based on a 100 nsec pulse 
where some evidence of a flat peak is visible. Figures 5-12 through 5-14 on 
the other hand are based on a 20 nsec pulse for the same conditions and here 
no flat peak exists. 

Flight 5 (Figure 5-18 through 5-30) shows true pulse limiting as 
evidenced by the long trailing edge when the altitude of 10, 000 ft is used with 
a 20 nsec pulse. This should be compared with the expected values of Figure 
5-10. Flight 7 (Figures 5-24, 5-25, and 5-26) shows increased pulse limiting 
at an altitude of 15, 000 ft for the same 20 nsec pulse. Flight 7 also shows the 
slight variations in exposure by using three values of pulses per frame (50, 

147, and 278). 

Flight 10 (Figures 5-27, 5-28) shows 10 nsec pulses over very calm 
water at 10, 000 ft. 

Flight 14 presents a comparison between pulse shapes as a function of sea 
direction. Figure 5-30 is taken with the polarization of the transmission in the 
direction of the sea whereas Figure 31, Figure 5-32 and Figure 5-33 are at 90, 
180, and 270 bearings with respect to sea direction. It is clear that wave shapes 
are not related to sea direction. 

Flight 14 (Figure 5-34) also shows the results when reception is cross 
polarized from the transmission. The signal level here was into the noise 
indicating greater than 30 dB between direct and cross polarized returns. 

The returns from ice at low altitude (2, 500 ft) are shown in Figure 

5-35. 
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Figures 5-12 through 5-37. Flight Pulse Data 
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RUN: 
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P/F: 

PW: 

BRNG: 


LEGEND 

Flight Number 

Data Run for That Flight 

N 
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Transmitted Pulse Width in Nanoseconds 
Bearing 
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The pulse shapes of flights 14, 15, and. 16 (Figures 5-33, 5-36, and 5-57) 
show that winds of 12, 6, and 22 knots, respectively, have little effect on pulse 
shape. 

5. 2. 3 Individual Pulses 

Figure 5-38 represents traces of individual pulses copied directly from 
the negatives. This procedure was necessary because the faint traces on the 
negatives were washed out in the reproduction process. Figure 5-39 represents 
10 pulses per frame which are barely visible. 

It is clear that considerable amplitude variations are visible in the in- 
dividual pulses and that multiple traces are actually overlaps of individual 
traces. In general, the leading edges of the individual traces are fairly linear 
(see Appendix E) and peak at the pulse width. The trailing edges, however, 

show wide amplitude fluctuations. 

* 

5. 3 Ocean Truth 

Two methods were used to obtain ocean parameters on all flights - 
observations from a NASA ship stationed in the test area in the area, and 
Stilwell photographs. In addition, on flights 5 and 6, a Cessna aircraft from 
Office of Naval Research flew alongside with a laser profilometer aboard. 

The ship observations were used in reducing the data because they repre- 
sent the most complete available information. Comparisons of ship observations 
with the laser data on flight 5 (see Figure 5-40) show very good correspondence. 

The Stilwell photographs (Figure 5-41) offer the greatest potential for 
ocean truth measurements because, when reduced, they provide a two dimensional 
spectral representation of the ocean. However, the techniques and computations 
associated with reducing the data require further development and evaluation 
before it becomes an acceptable ocean parameter measurement tool. Two 
photographs from flight 6 were reduced to provide the two dimensional ocean 
spectrums shown in Figures 5-42 and 5-43- These were in turn converted into 
line spectra along the dominant wave direction (see Figure 5-44). 
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Figure 5-39. 
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Figure 5-41. Stilwell Photographs 
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Figure 5-41. Stilwell Photographs (Cont. ) 
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Figure 5-44. Energy Spectrum for Raytheon Flight No. 6 


These line spectra agree quite well with line spectra taken by the laser pro- 
filometer on the previous day (see Appendix B). 

The results of the Stilwell process of ocean spectra measurements are 
discussed in detail in Appendix B. The technique offers great promise because 
of the ease of implementation,but further work on processing methods is required 
before it can be used to provide valuable and reliable outputs. 
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SECTION 6. CONCLUSIONS AND RECOMMENDATIONS 


6. 1 Radar Cross Secti on 

The results do not indicate that average cross section (cr® ) is functionally- 
related to the large ocean wave parameters, i. e. , wave height {Figure 5-1), 
wave period (Figure 5-2), or wave direction (Table 5-1). A functional relation 

appears to exist between <r° and wind speed. This seems entirely credible in 

<• 

light of the possible relation between wind speed and the generation of small 
high frequency capillary waves. This is particularly logical when one con- 
siders the dimensions of the capillaries {1. 7 cm or less) and the dimensions 

1 , 2 

of the transmitted X-band frequency (3 cm). Mr. John W. Wright has shown 
that X-band backs cattering at angles other than normal are highly dependent 

3 

on capillary waves. Kinsman (see also Appendix C) has also shown that the 
slopes of capillary waves can be very high, even exceeding 90° , and that the 
average slope of waves is more affected by capillaries then by the large wind 
waves or swells. 

The results also show that variations in cr° over the looking angles of 0 to 
4 degrees are small at the various wind speeds (see Figure 5-4). This agrees 
in part with the curves developed by Moore and Schooley (see Figure 5-5). 

Distributions of pulse amplitudes, and hence c® , appear to be Rayleigh 
and' independent of ocean parameters. 

Cross polarized returns at normal incidence were found to be at least 
30 dB down from the directly polarized returns. 

6.2 Pulse Shape 

Pulse shapes were recorded on film showing the relation between individual 
pulses and the integrated effect of many pulses. The individual pulses showed wide 
fluctuations in the trailing edge, especially immediately after the peak. The leading 
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edge of the pulses, however, showed a consistent ramp with little fluctuation. The 
start of the leading edge did not vary significantly from pulse to pulse. 


The multiple pulse tracings showed an . integrated average effect which 

4, 5 

agrees with calculations performed by others 1 . The pulse returns show an 

integrating effect from the time of the first pulse returns up to the pulse width, 
and then a decay thereafter. The impulse response can then be determined as 
follows: 


h(t) = k 



= k e 


-t-T / T 


0 < t < r 


t > r 


(6-1) 


where 

f(t) 

T 

T 


k 


.= the input pulse 
= pulse width 
= decay time constant 


2 , 

8 T h 

2.77c 

= a constant 


(see Ref. 4) 


which includes cr° and all other radar range parameters. 
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From this we can derive H(CJi) by performing the Fourier transform 


H(cd) 


K 

j (tD-CD ) 
o 


where CD = l/ T 
o 



(6-3) 


(6-4) 
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The above relations assume that the beamwidth of the radar is large com- 
pared to V cr/h , that a pulse limited condition exists . In actual calculations from 
the reduced data it was found that the transition region between pulse limited and 
beam limited condition occurred at an angle which was twice as large as the 3 dB 
beamwidth angle. 

In summary, then, average pulse shapes agreed with that predicted by the 
theoretical computations. The peak pulse amplitude is determined by a ° . The 
rise times are linear and equal to the pulse width and the decay times are related 
to altitude and beamwidth. No changes in o’ with looking angle were found, and * 
whether these could also affect decay time is unknown. Wave shapes were 
found to be independent of sea direction. 

6. 3 Ocean Parameters 

As in all experiments involving reflections from a target, the characteristics 
of the target are never as well defined as they could be* It appears that measure- 
ments of capillaries are needed to fully determine afunctional relationship between 
radar and ocean waves. 

Of the ocean measuring techniques used in this experiment, the Stilwell 

6 

process offers the greatest promise. It may provide the most comprehensive 
quantitative description of the ocean surface and could be operationally easy to 
implement. However, the total process is still in a research and development 
stage, and requires further theoretical and experimental evidence before it can 
be an accepted tool. 

The good old-fashioned 11 eyeball 11 technique of measuring ocean waves 

again yielded the most reliable and consistent results, although lacking in defin- 
ition* accuracy, and resolution. Laser profilometer measurements of wave spectrums 
did provide some valuable measurements,but the instrumentation and data reduc- 
tion is more complex than Stilwell photography. Laser profilometry is, however, 
a more proven technique. 
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6. 4 Recommendations 


Further radar backs catter measurements should be made at X-band to 
verify that a relationship exists between radar cross section and small high 
frequency waves (capillaries). This is not primarily a radar problem, but it 
is an ocean parameter measurement problem. Techniques for making the 
ocean parameter measurements should be further developed, especially the 
Stilwell process. 

An altimeter data processor should be developed which takes into account 
the leading and trailing edge characteristics discussed in this report. Such a 
data processor, or range tracker, should take into consideration the wide ampli- 
tude variations in the trailing edge near the peak, as well as the relatively minor 
amplitude variations on the leading edge. 

The next phase of precision satellite altimetry experimentation should 
be started. This involves accurate range or altitude measurements on a 
ground test range in order to calibrate equipment delays, and in an aircraft 
in order to measure altitude biases. 

Further flight tests should include correlation measurements between 
pulses, waveform sampling techniques, and evaluation of candidate data 
processors. 

The data gathered in this experiment should be further analyzed in order 
to determine rise and fall times as well as the distribution of amplitudes as a 
function of time. 
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APPENDIX A 


COMPUTER PRINTOUT OF FLIGHT TEST DATA 

This appendix contains the results of the data processing program described 

in Section 4. 1.4. It includes a listing of all the ocean and radar parameters for 

all of the flights, along with 0 ° calculations for all of the selected frames. Average 
o. 

a s, standard deviations, and frequency distributions were calculated for various 
groupings of 0 . The results are divided into five sections as described below. 

A-l. o° Per Frame 

This section includes significant ocean truth and radar parameters for 
each flight and for each run of each flight, and for each frame the calculated 
received power, target cross section, and radar cross section per unit area 
(all values in dB) using both the average return pulse peak and the absolute 
return pulse peak ih the calculation of ff°. 

At2. Average q 0 Over All Flights 

The mean o’® and its standard deviation and sample size calculated over 
the ten selected flights are presented here. 

A-3. Average 0 ° for Each Flight 

This section includes the mean (f°, its standard deviation, and sample size 
calculated for cr° based on average power return and for o-°, based in turn on 
absolute peak power return for each flight. The one standard deviation boundary 
limits are included for both values of 0 °. A histogram and cumulative probability 
distribution have also been calculated. The histogram is divided into two parts, 
the first is a frequency distribution over the spectrum of cr° values with each 
bar of the histogram being one half of a standard deviation wide, and the second 
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is a frequency distribution for values of tr° within the one standard deviation 
limits, with each bar being one eighth of a standard deviation wide. The cumu 
lative probability was derived from the first histogram. 

A -4. Average tr° for Selected Runs 

The output presented here is similar to the input of the data of Section 
A-3, except that the results are based on selected runs rather than flights. 
This shows that <r° is independent of sea direction. 

A-5. Average cr° vs. Pulses Per Frame for Each Flight 

Here the average cr° is calculated for each variation of pulses per frame 
during each flight. The values of pulses per frame were 1, 50, 148, or 278. 
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APPENDIX A. FLIGHT TEST DATA ANALYSIS 
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PULSE WIDTH(NS) 

PEAK POWER 

POLARIZATION! 1 1 

C54 HEADING 

1 

5000 * 

1 


10 

3 

1 

180 ‘ 

-2 

-51LQU 

JL 


ID 

3 

] 

90 

3 

5000 

1 


10 

3 

1 

360 

4 

5imiL 

L 


1S> 

3 

L- 

213 

S 

5000 

2 


10 

3 

1 

180 

6 

bOOO 

2 


LD : 

3 

- 1 

90 

7 

5000 

2 


10 

» * 3 

1 

360 

‘ 8 

SMJQ 

_2 


IP 

i — 

L 

- .210 

V 

5000 

10 


10 

3 

1 

160 

1_0_ 

b-QM. 

.,11), 


ID 

1 

1- ' 

90 

11 

bOOO 

10 


10 

3 

1 


1? 

5000 

10 


IQ 

3 

1 

270 

1 J 

5000 

bO 


10 

3 

1 

180 

LA 

5Qi)J L 

bD._. 


ID 

3 

1_ 

. 93 

IS 

500 0 

50 


10 

3 

1 

360 

_Ui_ 



5D_. 


JLO 

2. 

JL 

270 

1 r 

bOOO 

148 


10 

3 

1 

160 

lb 

bOOO 

14 8_ 


ID 

3 * 

1 

90 

19 

5000 

148 


10 

3 

1 

360 

23 

50JUI 

Xhh^, 


JlQ 

2 

1_ 

. .2 IP 
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ERii tl£ KCYLL IDHMLl Ay ,R,cyo_ HWK (L)BH) PEAK S LGMA- AY 5USttA_B£ AK £. IGMAZ. l£EU_AY SlGI 4 AZ.UlJai_P.EAK , 



13 

l 64 9 

-H4.01 -80. B7 

49.44 

52.58 

7.012 

10.953 


1 3 

JJlbiL 

"ddiil? -79.69 



9^aox 

12a. 121. 


13 

l*bl 

-82*33 -79. S3 

51.12 

53.92 

9.409 

12.289 


) i 

1 67b 

-76.70_ -74^83. 

5J>,J5 

5ji._i.2_ 

15.116 

16.989 


]<4 

1690 

-80.80 -77.89 

52.64 

55.56 

11.014 

13.927 


l n 

1 *4 1 

-79.67 .... -76.74 

53.78 

_SfiuJU 

12*148 

--.15a.Q7.7 


1A 

1S4P 

-80.34 -76.92 

53.11 

56.52 

11.480 

14.894 



17^12_ 

-8?. 30 -78*55 ___ 

_ 3XU5- 

._5iufi9- 

^9 ^23. 

L3*264l 


IS 

1 / 34 

-81.14 -78.06 

52.31 

55.39 

10.677 

13.763 



| 73* , 

-8iijd_U -77.46 . _ _ 

53.35 _ 

._55^49_ 

UX716. 

13.862 


)b 

1 77A 

-80.11 -76.79 

53.34 

56*66 

11.708 

15.030 


* 1 b 

1 / /S 

-74.1/ -77.45 

54.28 

56.00 

12*640 

14t367 


16 

177* 

-7 7.99 -7b.7<2 

55.46 

57.73 

13.826 

16.102 


] 7 


-77 f B6 -75.14 

55.59 

_5£*3X 

X3*25ii 

X*jl6JX 


17 

Id 17 

-75.43 -73.31 

57.52 

60.14 

15.606 

18.509 


1 t 

1*18 

-74. S3 -71.54 

58.91 

6J..91 

17*284 

20.282 


17 

1819 

-75.41 -72.0b 

58.04 

61.40 

16.409 

19.768 


1 7 

1 820 

-75.85 -73.32 

57.60 

60.13 

15.968 

16.500 

■ 

1H 

IrtSH 

-79.45 -76.16 

► 54.00 

57.29 

12.368 

15.656 


1H 

Jhb9 

-79. 8U -73.9J1 

5J.65 

53*4.1 

.UjOi* 

17.ft.36, 

u 

IS 

1880 

-74.^8 -76.01 

53.77 

57.44 

12.139 

15.812 


18 

1 H8l 

-78.U -75.16 

55.14 

58.29 

13.507 

16.656 


IB 

" 186? 

-78.78 -75.74 

54.67 

57.71 

13.035 

16.079 


W 

1900 

__.-79.86 -7 7.14 

53.58 

56i31 

1U954 




14 

1401 

-79,49 -75.99 

53.95 

57.46 

12.325 

15.826 



„ 1 40? 

-80,45 rJiL. 66_ 

52.99 


XU364. 

iSa.158 

% 

19 

~ 1403 

-80.29 -77.17 

53.16 

56.28 

11.529 

14.647 


1 4 

_L40A_ 


54.17 

._s:u4 a 

1/U544. 

15_*TB2. 

• 

20 

194? 



-HU.?5 -7b. 6b 

53.20 

S6.79 

11.573 

15.161 


?0 

1 443 

-80.40 -75.82 

53.05 

_57»63 

Ll*423 i— 

lfr«QQ3 


?t» 

1444 

-79.43 -76.08 

54.02 

57.37 

12.388 

15.736 


?0 

_J_94b 

-HO, 43 -76.78 

52._52_ 

56.67 

10.086 

15.038 


20 

1946 

-79.47 -76.39 

53.98 

57.06 

12.349 

15.427 


FL IGHr ShA 

TiT^tclT 

J.\l WAVE hEJbhl(M) wAVt PtK 1 OU ( SEC ) 

WINU SPEED ( KT ) WIND DIRECTION* 

WATER TEMP <F > AIR 

TEMP AIR PRES (MB) 


1 1 494 


W 94 

99 999 


99 99 

9999.90 


RUN 

a] 1 I 1 u )t (h 1 ) PUl Sts PFK FRAME HUl SL 

_WJUJ H(_NSJ PEAK .POWER. 

PpkAH LZ AJJLQtf JULL _ _ _C 54 _HE ApXNfi 


i 

no nAT'\ cot i »• 


7 FUAMfc KC Vt) PW&(i)HM)AV KCVD BWH < L)BH ) PEAK SIGMA AV SIGMA PEAK SIGMAZ(OB) AV SIGMAZ(OB) PEAK 


P CTx^OhT SP“A “TfWirCT TTJm WAVE MtTbHKFn WAVE PEPlOUIStC) mT I ND SPEED < KT ) WIND DIRECTION WATER~7EMP (F ) AIR tEmP SlR PESSIMIST 


lg 94 99 IS 330 99 99 9999.90 

Hij in AL! JLLU jffmjme, pulse jvip_thjn.s_i peak, _power PJ>yw iXAJXcn y _> c^s^jheading^ 


l * 

? 

mono 

1 1>0 Ilu 

l 

l 

20 

20 

3 

3 

1 

1 

355 

J 

1 nooy 

i 

20 

3 

1 

265 

4 

1 0 0 OU 

i 

20 

3 

1 

175 

5 

10000 

2 

20 

J 

1 

85“ 

6 

loooo 

2 

20 

3 ■ 

1 

355 

{ 

loon o 

2 

20 

3 

1 

265 
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H 

1 OOOo. 

2 

20 

3 1 

175 




10 

loooo 

luoru 

10 

10 

ZT5 

20 

J 1 

3 1 

85 

355 



~ 11 
12 

ioot)o 

1 i) Ol)o 

10 

10 

20 

20 

3 1 

3 1 

265 

175 


■ 

1 1 
14 

10 0i)O 
1 00 00 

50 

SO 

20 

20 

3 i 

3 1 

85 

355 



IS 

loooo 

50 

20 

3 1 

“ 265 

- 

** 

lb 

1 OOOil 

bO 

20 

3 1 

175 


'' - — ■ 

17 

10000 

1 4M 

20 

3 1 

85 

<icc 



18 

10 000 

14H 

20 

3 1 

355 

— 

— — — ' ■* 

19 

1 00 00 

140 

20 

3 1 

265 


w5 

20 

1 00 00 

140 

20 

3 l 

175 




21 

10000 

27b 

20 

3 1 

85 


— » 

22 

1 0 C» ilO _ 

2 7H 

* 20 

3 1 

355 

1 

— 

2J 

10 0(10 

2/B 

20 

3 1 

265 



24 

loooo 

2 7H 

20 

3 1 

175 

o 



F*AME 

KCVi) PWPfDBMftAV 

RCVU PWR < DBM 1 PE AK 

SIGMA av sigma peak 

SIGMAZ(DB) AV SIGMAZ(DB) peak 



] 7 

?642 

-6 4.S2 

1 

1 

1 

1 

1 

1 

i» 

If 

1 

1 

1 

1 

61.97' 

64.61 

14.310 

16.955 

1 / 

?64 3 

-62.6b 

-61.21 

62.84 

64.20 

15.193 

16.626 

17 

2644 

-62.74 

-60.73 

62.75 

64.76 

15.097 

r tttho 

1 7 

?64b 

-6J . 1 2 

-60.69 

62.37 

64.00 

14.723 

17.146 

1 7 

264b 

-hi. 12 

-60.58 

62.37 

64.91 

14.717 

17.260 

10 

2b89 

-65.43 

-62.91 

60.06 

62.50 

12.400 

_JU.?25 

1 H 

2690 

-64. «4 

-62.36 

60.65 

63.13 

12.995 

15*477 

1H 

P691 

-64.43 

-62.30 

61.06 

63.16 

13.407 

„__X^53A 

1H 

269? 

-63.95 

-62.12 

61.54 

63.37 

13.888 

lb* 723 

18 

2693 


-62.74 

60.85 

62.74 

13.206 

15*094 

19 

?/2S 

-6h • 3t> 

-62.74 

61.13 

62.75 

13.481 

15*100 

19 

21 2b 

-b6. 73 

-6J.91 

58.7b 

61.58 

11.104 

13.932^ 

19 

2727 

-6b. 26 

-62.67 

59.23 

62.02 

11.575 

15* 168 

JL9 

J>7i8 

— 66 . (If) 

=62^6, 

59.49 

63-33 . 

XI .839 _ 

15* 679 

1 9 

? /29 

-64. 15 

-61.36 

61*14 

64.13 

13,485 

16.477 

20 

2 /67 

-67.42 

-63.19 

58.07 

62.30 

10.419 

14*650 

20 

2768 

-6b. SO 

-62.39 

59.99 

63.10 

12.340 

15.448 

20 

2769 

-6b.si6 

-62.86 

59.53 

6JJ.63 

11^878 _ 

X'tiV*. 

20 

2770 

-bb.78 

-62.03 

59.71 

63.46 

12.060 

15*813 

20 

2771 

^66il>2 

-62.13 

59.47 

63.36 

11. 821 

15*706 

R IGHT SEA 

DIRECT ION 

WAVE HEIGHT (H) 

WAVE PLklODtSeO WIND 

SPEED (KT) , 

WIND DIRECTION 

WATfcR T£MP<F>_ 

AIR TEMP .AIR PRES (MB) . 


** M 17 -h^*76 -60. /3 


o 

o 

o 

o 


13 


22 o 


26 


220 


50 


55 


1005*90 


RUN 

ALTITUDE < F T ) 

PULSEb PER FRAME 

PULSE WIOTH(NS) 

PEAK POWER 

POLARIZATION ( 1 ) 

C54 HEADING 

* 

i 

i 

1 

! 

I 

J 

I 

1 

1 

1 

10 Mil 

loooo 

LDiliUl 

JL 

1 

_I 

1 

1 

1 

1 

1 

1 

i 

oj=> = 

r\J<\j ru 

1 

1 

a 

A 

L 

1 

1_ 

.220 

130 

itP 


4 

s 

loooo 

ft OilOO 

1 

10 

20 

20 

3 

-3 

1 

1 

310 

22 St 


6 

100U0 

10 

20 

3 

1 

130 


7 

UlQlliJL 

10. 


3 

L. 



A 

10000 

10 

20 

3 

1 

310 


V 

IOjMjI 


.20 

A 

JL 

22P 


10 

loooo 

50 

20 

3 

1 

130 


1 1 

1 0000 

b(> 

20 

3 

I 




l? 

loooo 

SO 

20 

3’ 

\ 1 

310 


13 

iimui 

X4H 

dP 

.3 

X 

2 20 


• 14 

1 0 00 0 

148 

20 

3 

1 

130 , 


IS 

JLOillLil 

J_4M 

£0 

A 

X 

40 


16 

* loooo 

148 

20 

3 

1 

310 


17 

1 00 0 0 

2 M 

20 

3 

l 

—220 


10 

10000 

2/0 

20 

3 

1 

130 



liUiiLy 

?/8 

jLU 

-- — A 

1 

40 




A- 14 




20 

ioooo 


er* 

20 


3 

1 

310 





21 

iOuon 


50 

20 


3 

2 

220 






22 

Fudoo 


so 

20 


3 

2 

130 





pi 

t nono 


so 

20 


J 

2 

40 





24 

) 000 0 


so 

20 


3 

2 

310 



— 

— 

— 

FKAMfc 

"kcviT 

UJHM) AV 

KCVl) PWR l DbM ) HtAK 

SHjM/TaV ~~~ 

SIGMA PEAK 

S1GMAZ <0B1 AV 

SIGMAZ(DB) PEAK 

— 


FI I6HT 

SFA 

UTRElCl In* 

"mve 

ME | liH 1 (F r> 

WAVE PLR 10U ( SE C ) 

W I NO 

SPEED (KT> WIND 

DIRECTION 

WATER TEMP(F) AIR 

temp air PRES C MB > 



14 

_ JSO. 


b 

* 

5 

_ J_2 

360 


55 43 

1020*30 




RUM AI-niuMhfH) PlM.ShS PhR FRAME PULSE WinrH(NS) PEAK POWER POLAR 1 Z AT I ON ( 1 > HEADIN6 


1 

, 2 

3 

. 4 

1 0 0 u <f 
100o0 M 
1000 0 
j o o n u 

1 

1 

1 

1 

20 

20 

20 

20 

12 

12 

12 

12 

l 

1 

1 

1 

— 

350 

260 

170 

fio_ 


l > 

lOOtjO 

10 

20 

12 

1 


, 350 

‘ 

* f> 

10 0 0 0 

10 

£0 

IZ 

1 



260_ 


( 

lOOo-P 

10 

20 

12 

1 


170 


. H 

LUMlL 

IP 

2P 

IZ 

1- 



&p 


4 

j 0 0 0 0 

bO 

20 

12 

1 


350 


10 

10000 

50 

20 

12 

1 


2frQ 


11 

101)00 

bO 

20 

12 

1 


170 


12 

1 uOOO 

Si) 

20 

12 

1 



80 


13 

1 Ooon 

J AS 

20 

12 

1 


350 


« 1* 

1 00m) 


iLO 

JiL 

JL 



_g60 


15 

IOOImj 

1 46 

?0 

12 

l 


170 


. \ h 

l on oo 

14H 

20 

12 

1 


&0_ 


u 

10 (MM) 

so 

20 

12 

2 


350 


t 1* 

looou 

Si) 

2U_ 

1-2. 

Z 

.. 

_260 


1^ 

1 OCfliO 

so 

20 

12 

2 


170 


, 20 

1 Oj.MMl 

bO 

20 

1J> — 

2 



8 0_ 


21 

10000 

50 

20 

12 

1 


350 


22 

lft/10 0 

bO 

20 

12 . 

1 


260 


?J 

1 0 0 0 0 

SO 

20 

12 

1 


170 


fi* 

lJJOOjL 

so 

'dP 

IZ 

l 



80 

— 


FRAME 

KCVU PrfR <L>HM) AV 

RCVU HWR (L)hM)REAK 

SIGMA AV 

SIGMA PEAK 

SIGMAZ(OB) AV 

SIGHAZ(DB) PEAK 


l 

l v 

-h2.P7 

-62.27 

57.20 

57.20 


9.549 

9.549 

1 

44?0 

-h 4.44 

-63.44 

55. 98 

55.98 


8.331 

8.331 

1 

~V~? l " 

-M. M 

-63.' 73 

55.74 

55.74 


8.089 

8.019' 

1 

WdZ 

-ns.s/ 

-6b. 57 

53.90 

53.90 


6.250 

6.250 

1 

3423 

-63. s/ 

-63.57 

55.90 

55.90 


8.245 

8.245 

l 

^425 

— h3 • 

-6 j.64 

55.83 

55.83 


6.178 

0.176 

1 

1426 

-ho rn 1 6 

■" -66.16 ~~ 

53*31 

53.31 


"57555 

5.658 

\ 

34^/ 

—64 • 54 

-64.54 

54*92 

54.92 


7.274 

7.274 

1 

142 H 

-h~2 » h 7 

-62.67 

56.60 

56.60 


8.950 

8.950 


*424 

-6b • 04 

-65.04 

54.43 

54.43 


6.782 

6.782 

1 

3430 

-hb*4rt 

-66.46 

52.49 

52.49 


4.837 

4.837 

1 

3431 

-hd.TU 

-66.30 

51.16 

51.16 


3.513 

3.513 

1 

34jp 

-68.^4 

-66.44 

50.53 

50.53 


2.875 

2.875 — 

1 

14 43 

-66.36 

-66 . 3b 

5J.11 

53.11 


5.463 

5.463 

1 

1434 

-hb.?U 

-65 • 20 

54.27 

54.27 


6.619 

6.619 

1 

3435 

-63. 4l 

-63.41 

55.55 

55.55 


7.904 

7.904 

i 

3434 

-64 • S3 

-64.53 

54.94 

54.94 


7.291 

7.291 

1 

3437 

-6s. .14 

-65.09 

54.36 

54.38 


6.732 

6.732 

l 

34 3B 

-h 1747 — ~ 

-63.47 

s575o 

5S.50 


7.845 

7.845 

1 

3434 

-6h .48 

-66.46 

52.99 

52.99 


5.337 

5.337 

1 

1440 

-64.54 

-64.84 

54.58 

54.58 


6.929 

6.929 

\ 

3441 

—63 .70 _ _ 

-63* 70 

55.77 

5£i_77. 

— 

8^121 

, 8.121 


A- 15 • 



l 

l 

1942 
394 3 

-6b, 12 
-64. 79 


-6b. 12 
-64. 79 

54.35 

54.68 

54.35 

54.68 

6.701 

7.026 


6.701 

7.026 


y 

41RO 

-61 .93 


-59.07 

57.54 

60.40 

9.080 


12.746 


9 

4 1 H l 

-6«i * 46 


-Srt .bO 

.58.62 

60.97 

10.971 


13.316 


9 

Mb? 

— 6 l • “>K 


-59.07 

57.88 

60.40 

10.234 


12.746 


10 

4?21 

-61 .'64 

< 

-bb . 44 

57.78 

61.03 

10.129 


13.383 


10 

4 22? 

-62.63 


-5b. 9 7 

56.64 

60.50 

8.905 


12.846 


10 

4221 

-62 . 5b 


-S9.36 

56.92 

60.11 

9.268 


12.457 


11 

4263 

-6U .4.1 


-58.05 

59.04 

61.42 

11.392 


13.770 


11 

4264 

-60.60 


-bH .12 

58.87 

61.35 

11.215 


13.700 


1 1 

426b 

-60.44 


-57.96 

58.53 

61.50 

^ 10.083 


13.654 


1? 

4 JOS 

-62.17 


-bH. 73 

57.30 

60.74 

9.649 


13.087 1 


1? 

4306 

-62.0b 


-5b. 49 

57.42 

60.98 

9.773 


13.327 


1 2 

4J07 

-62.4b 


-59.1b 

56.99 

60.29 

9.337 


12.637 


1 1 

4347 

-57.46 


-54.9b 

61.99 

64.52 

14.342 


16.66E> 

* 

13 

434H 

-5b. 0 1 


-55.67 

61 .46 

63.79 

13.807 


16.144 


13 

4349 

-bb.30 


-55.43 

61.17 

64.04 

13.519 


16 .387 


n 

43S0 

-5b .44 


-55.10 

61.03 

64.37 

13.377 


16.715 


n 

43b 1 

-57.31 


-54.05" 

62,16 

65.42 

14.507 


17.772 


]4 

43H9 

— b 7 • 4 7 


-55.07 

61.99 

64.40 

14.344 


16.750 


14 

4390 

— S 7.0b 


-54. bb ' 

62.41 

64.59 

14.757 


16.937 


14 

4391 

-Sb.H7 


-55.82 

60.60 

63.65 

, 12.952 


15.997 


14 

4 392'"' 

-55. lb 


-55.49 

61.32 

63.98 

13.666 


16.328 


14 

4i93 

-57.47 


-54.99 

62.00 

64.48 

14.351 


16.827 


lb 

4431 

-58.50 


-54 • 54 

60.97 

64.92 

13.321 


17.273 


)b 

4432 

-bd. 73 


-bb.B2 

60.73 

63.65 

13.084 


15.997 


lb 

4433 

-58. so 


-56.81 

60.97 

63.66 

13.321 


16.012 


15 

4434 

-5H , 0 7 


-55 . U6 

61.40 

64.41 

13.746 


16.7S7 


IS 

4435 

-57.6b 


-5b • 35 

61.82 

64.12 

T4. 166 


16.469’ 


16 

4473 

-b4.09 


-5b.b4 

60.38 

62.92 

12.728 


15.274 


is 

4474 

-54.3J 


-56.52 

60.16 

62.94 

12.511 


15.293 


16 

44 75 

-5b. 91 


-bb. 72 

60.56 

62.75 

12.909 


15.101 


16 

4476 

-59.66 


-57,08 

59.81 

62.39 

12.161 


14.7^ 


_l6_. 

4477 

-59.59 



-56.60 

59.88 

62.87 

12.232 _ 

— 
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-h2.s0 

-Sy.ys 

56.97 

59.52 

9.316 

11.865 

lh 


-h 1 . MA 

-sy .63 

57.53 

59.64 

9.882 

11.991 

1 MOI 

Av! / Ml 1 ll'*l 

CUUt- < l L'lKtC 1 *2 CROSS 






NOTI-JPU^ H U.»n S^b PI-k hrfAric v/AKitn UUKlNU KUN With VALUES UP 5U*14B*278 



A- 19 


A-2 OVFRALl FLIGHT* 


MFAN STfiA {IJH ) STn.DEv. S^iRLt Skt MtAN SIGP(DB) , SID. DEV. SAMPLE SIZE 


16*115 


14.6^ 44/>.i» 


19.171 4M.0 


^ FLTGHT SltW hjt<> AV SAMPI E Si/'t bTU.UtV. tOfcU MEAN S1GZ(DB)PK SAMPLE SIZE STD.DEY. (DB) ' 

’ — “<T hTtT i 7 iToT ir.88 la - 

— f7,/A V^+ ~ l M ~"iTcV^I J H TlGZAV - i b I I) UEV UH SIGZPK ♦ l STD UEV OB SIGZPK - 1 STD DEV DB 


FLIGHT HISTOGRAM (NO. OF Mh.UtV.FKUM MEAN | 


histogram .from - 1 _rn_ UH _ . W / ri J.U_ ^JU 1/A _t D_ J_ _S IP_±J l£Vj 

00 1 0 U 0 1 2 1 10 0 12 2 0 

CUMULATIVE PKOHArt I L II Y 

0^0000 (j fU 0 _ 0 _i> n_u ^l_h(^62 i 4JJLT3 •.55556 • 13333 1 #P 0 0_0 0_ 


1.00000 1.00000 




FLIGHT HISTOGRAM (NO. OF SfO.UtV.FPOM MEAN 


£LJ ght _( y-^1 _<.=/_ IIL.-J -Pi _J riUl. JfJ> _r*>.L -!=.<£ _jo..o L io_I9_* 5>_ i- 5 jo. 1L u_IP_i. s> _u .5_to _?i _JGE_2i _ . 

... j a 1 4 1 5 O Q 


Jj LST QGQ-AM t£yJUl _=X AiL-=JjH±s=Ju 
1 1 V 0 l 


t^-^a_icL,x^mai£X. 

0 0 0 2 1 11 


CUMULATIVE MKOMAmII ITy 
0 . 0 0 0 no . Oiiitl 


. Silt . 


FLIOmT ^ftN Mh/(uM)AV 


“7^ix“7zt bTD.OEV.IDB) MEAN SIGZ(0B)PK SAMPLE SIZE STD.DEV • (DB) 


mr./AV ♦ 1 SM UhV ni» 


TlUZAV - l blU DEV Ob SIGZPK ♦ 1 STD UEV DB SIGZPK - 1 STD OEV DB 


FI IGHT Ml SI d()PAi« (NO .Oh Si li.UEV .HPOi ‘1 ME AN 


Fi i qjx j i_J .-p ) t -s In - t i rJ X=l_ LlX^J. JL XtJ. _T.il _r _I X il L_ i_0 — ID *jS.L -L»-^- -1 0- J-l— U — 1-0 — l^Sl — (2 *3_JD. -21 — LGE — 2J — 


_UI — Ial ^ Um.xjsJJjX -lil_s J/JtA==XZH _[ U_ X Ul^llE r 5 T f 

? [ J | 1 O 1 o 0 2 0 1 o 1 l 


CUMULAl 1 Vf I I v 

it . tin (i'll. f. - jh, i«hj 



A-20 


FLIGHT h l ('Jo.iJF m.mkvpK^om mpan 


5 jl ] l A 3 4 1_ 3 10 


HISTOGRAM. t-upM -j hj w/d.-7/h' f0 7/8 TO I Sin.PEv. 


1 2 


l 


i 1 


0 


1 


CUMULUtVf P^tw'.Mll If/ 

ruMiiiUL J jJi_U i 'A32-S3— ^50 £M *72222 *17273 sSJthhiL L*_0QO£P JUCUPfiPiL 




PI IGHT 1 C jhi/HH)AV SA*P LL SIZE STD.UK V. <DB > MEAN SIGZCDBIPK SAMPLE SIZE STD.DEV.(DB) 


12.36 


22.38 


16.16 


Sfh/A\/ ♦ 1 DhV \U* SK’ZAV - 1 b1l> DfcV DH SlGZPK ♦ 1 SID UEV UB SlGZPK - 1 STD DEV DB 


1 ~*rt 1 7 , fi2 

FI IGHT HIsTon^A f r i*l«O i * « I ii.ii)- .h ^ii*n ME/V'J 


23.31 


21.19 


. £t_I fiHL -1 Uj=S. L > *. X _S =.U“*_ -abtl _S - .£ a W. Ji. >_ JJ>_LL>_.£ LL.£J0.1LU_LO_Uy _U,5_T2 2) <5E g>. 

t _LL i l 2 1 3 0 1 1 0 

HIST0G'JA4.i*^o<* -i , n -j/*. 7 /m fu i blU.Dtv. 

0« 1 I 10 0 0 3 0 0 0 0 0 0*0 


CUMIILAJT Vh Ppoh,vHll [ 1 i 

0.«0ii«i» m <i,oi ,MlO 


_-_LUL l . J3J33 .44444 .77778 .77778 .88889 1.00000 1.00000 




FLIGHT * '(> «.i»r ** I »h *> .i* ^u«'i Ai\‘ ~ 


1 


-H-.1-QH IT H 1-21 <-• ' » t -uM Mo 10 -l ) (-] 10 (- t b TO 0) (0 TO .5) US TU 1> (1 TO l.b) U.S TO 2) (GE 2) 


31^JSL c i^Iil^2iA—ZiJ2 — UA\il/2 J 2.2 i£«j2r2J23 3.2. A. y_._ 

11 / 2 i o o 1 1 o 


CHMUlATTVh Ilf 

O.OiMHM* u.oOOOO 0.00000 


^.A4A44 .66667 .77778 .77778 .88889 1.00000 1.00000 


*4 44 *4 44 444444*41 44444444 


PI TGnT ^hI>-"s|!w |i'n)AV SAMPLE MZt b 1 D.IJE V. (DH ) MEAN SIGZ(UB)PK 


"stoplf siz* snnuFrrmr 

78 TeT&7 


20 . 


76 


lb. bl 


24.01 


ST67AW + 1 DH SIG^AV - 1 S f D i)(£V UH SIOZPK ♦ 1 STO U£.V DH SlGZPK - 1 STD D£V DB" 


21. Mg iVf.OS 

FLIGHT hISTOGHum <ko.oI» s I D.nLv7rw(7M~^ AN 


25.10 


22.55 



A-21 


FtIRHr M T— ??> (-<■» TO -]♦•>) <-l.S lu - 1 ) (-1 10 -.S) <-,b Tu Qi (Q TO .5) (.5 T<i> 1) (1 TO i,5> (1,5 TO 2) <CE 2) 

7 1 3 J , J_6 13 17 0 6 0 

_ht sr on*? I/iiLTi/ii _i <L rrZ/Ji _IX ± oty* 

455*813360564410 

WhulattvV p^ohakh n’y 

IKUP009 .0 3H4o .07^1^2 bJh ._32Ub_i ♦487,1^8 .70513 .82051 .92308 L«JMOO_ KOOOOO 


FLIGHT HISTOGRAM (MO.UF S f U.Ut V .1- HUM m£an 

F.U QHX _1 LL£1 _J^g_ i_ JLtJ _J_U_r -XL J b. J0_ J) L, JLP_IP_- 5 )_ i_. b. _T0_ 1 )_ iJLlO_3La.5> _UiB _T0 _2Jl JGE_?JL 

7 ? 3 f 12 15 17 9 ; 6 5 2 

• r 

HI H$C££^1<Q£Xjl 1 jn rvn *^ / / h. _jcl -a^Az^rl/Jt J±L± J> Mn OF v^> z 

244,>644l'75231242 

— — — j — — — • — Yr77i ,—jYy 

0.00 000 .02584 .08410 .15385 >30769 >50000 >71795 >83333 .91026 >97436 1 >00000 

* uisjt 

rTt ghT 7757 ?j 7,F<7)8) “v SAtfQ”i7e s tl>7d7 v^Tob > hl7n“ si &z7 dbTpk sample “sTze “sfbTDEV - , To § ) 

8 14.34 161 " 9.20 " ” 1 0 > 19 135 * ‘ 13.70“ “ * 

Sfb7AV ♦ 1 <Tt) UfcV Do SIG2AV - 1 STD UEV Db TlQZP^^~SjiTDEy~OB Sl7zpfT7 7 ~s7o~ DEV ~03 

— — - - 17.75 \9 . b i ~ * 16 V 20 

FI IGHT HISTOGkAHLNO.OL •>JU,otV.h>aiM MEAN * • 1 " 

EXIGJHXai^r/X^^^lD^XoXJLril J ^J^_-Xl_JLnlJiLjifii_X^5.XQ-ilJ_iiL_m^5JUL.aXU-JLJ_U_aiJ^ii^J^_TD J21_AGE_^ 1 


js h l : u xi an a* aa n 6 s 

ID 

4 5 / ✓ h 410 0*1111 6 6 9 7 4 0 

aww7f?vt pH<777nTr ^ _______ 

PiOttim dLaJjLL JJtXlU *1X1X1,.., n47ti2fr iflflM .06335 .93160 „ >96094 __1«0P0Q0 


iil WiMUlAiUii ‘ 

FLIGHT h J 5Tni,vAJ» (ivM.iih 8 ([>.*)£ v >£ ^l)to •'IlaN _ ~ _ ”™* 

XUfL4l~JJ-X=2J._ C^li_-^>iX^U.^XJL=JLLX=a_m^^LJl^aU.JlXJLft_UL^XJ^JiLXL-lJL^Ii)„jL i 53_a>5_lD^XXQL_ai 

H jo i» 8 39 49 19 3 3 J> U __ ___ 

MISIQGaArtJfrtO*- "I In I D , 1 . blDiflhYi 

9 6 4| 5 •* U 18 9 6 8 4 1 1 2 0 0 

CUMUl AtTvF pwiihiiJ U IT r ~~ ‘ 

Jl tJllUliLU 1» iLLLiQo li AliiianiL aM 4*4 *23443 iii.31Q4 jJ*59£4l ItOMflX 





A-22 


RIGHT SAM^Lt Sl/t S TO.Ut V ., ( Urt ) MEAN 5IGZ<UB)PK SAMPLE SIZE STO.OEV. (OB) 

Tfi |““ *3 fu 730 1 SV99 " 3 J 13716 — 

s r u/ ft 7 ~7 7 TiT ur V~7i .< m<vav - i sro uev dh sigzpk ♦ l stu dev ub sigzpk - i sto oev db 

— ^7j7 17781 1277*8 

FLIGHT Ml STGGW/* »Mim >.{>r S I V .h *.’UM . Ht 

_FL_! gh_l JL^rL - 1 -_L *_'_L _1 zJ L _L-J _Jii »_>J_ J J±L A L_ U>_ IP_.5 Ll_4Jy.lL <J_ IP _1_* Si _(_1 .5_J 0 JJ]_ JGE _2_t 

_1_0 0_ 0 i * lb 3 1 | 2 2 

HlSTOfiPflM.fr tVUM - 1 TO 7/o H» -4/A. f/H lU 1 SlU.OEV* 

3 0 bill) 111000 

cum ul a fi 7“ i^onTiTTl T } ’ 

L'syJiyiilL tUJjMp *27273 *72.727 •JJLBld *93939 1 *_00000 

RIGHT HiSTofiKAMl'n.Ur > I lUDf. v • ! n>U n l£ AW 

* 

RIGHT II 1-2) f i (- 1 .*- fu -1) <-l fu -.b> (-.b lU 0) <0 TO -b> C.S TO 1) <1 TO 1*5) U.S^.TO 2) IGE 2) 

f ' ■ 

_J£ 0 <- < £> Iji 6_ _U L _2 _2 

HISTOfiPAK.fr *u<* -1 In -7/fr«-//rt m ; - j/t». 7/d lu j a HUjiL V_. 

ll<’2‘*j2~'iJ0210000 ' ~ 

CUMULATIVE PHmHiU li TT7 : ' ’ ' 

_ 0 iPiiPJ) 0 1* 2fiJ 'ii' 0^‘K 110 o. ._y 9091 .2727 3 *6666 7 ^84848 ^^848 .87879 • 93939 1 ^00000 

/ 

i> VhS Vb VfcJ ** *> ** 

FTTght mh “ mTT7 iihmv TTTTT’slTF sTu.utv. <dh> mean sigz(ub)pk sImplFTITE sToIBETTTBBT 
17771 2o 7777 TbVSS ”20 9717 

"s PH 7 a 7“ “i “>77 “tv ~i 7 ”i7,7aT”7i~" 71 717tv"uh “s i“g7p7 ““VstcT dev“ub s’ iG7p*K ““sto"*oev”db 

__ TTTTb lb. 68 U. 78 

FI IGHT HlSTiHiP'N iMO.Ofr •* I ij . . jfr v . r MfLAf* 


RIGHT ( I 7 -S ) {-s in -|,n) <-t.o lu -1) (-1 U) -.S) t-.’D IQ 0 ) (1) 10 .b) (.S TO 1) (1 TO l.S) (1.5 TO Z) <OE 2) 


1? U *• 2 z' 2 3 

2 2 

2 0 

V 


HISTOGPAfUf POM -l fu - r/K.-z'/n 7u -V'4. 7/rt 10 i SMl.UhV. 

13I20U1 iUJOO 

1 0 0 

1 



cdmul at i ve; huohah h, i r » 

0 • OCi f>0 0 0 « Oirti (ill n.iUMHKt « 1 (HI 0 0 .AbOUO • bbOOO 

.70000 .80000 

•90000 1.00000 

1.00000 





FLIGHT HlSTnGPA'MMi.oh <; I u .i»i- v .r pOH A\i an 



TOTAL f K 


£UGHT JLTzSi Jz** Jit JLL1°_I9 .j£L !•* JtlaSJtfiiSL JGEJM. 

12 0 1 l 5 7 0 2 1 3 0 

.it Hi stoqp a** £hom_ —i—Itl—J-/* zrl/A^JA j/Wr zrl/h J it J [ttv • 1 

02.4ui<*2200000101 

II ■ ■ _ 

CUMULATIVE PMiihAhU JTr 

tttjiiui amum asuM amum amum .damp iflSgaa itmoo IiOpp.p.p. 

1 nw yysy m s 


Fl IGHT w AM SIG/(lto)*V Sfti^Lt SIZE sriJ.OEV. 108) MF.J^ SIGZIDBIPK SAMPLE SIZE sfD.DEV.T08) 

_ 7o7v5 sb vTzl is. 38 32 T5755 

s F g Tav ~ - i~.vi if 7i Fv ~7) T sigZav _ -~i~siFTe"v _ ob 'sf6zi»K“i — std dev"1 5b slGlPTi^nrsTD Devdb 

_ n.~fi * lT.Tir ll .“aa 

Fl IGHT "ilSTOGHi'MNO.OF S l U.llEV .F t>OM MEAN 


, FUfi HT (| T-? ) l-P TO -I. a) <-].i> TO -U 1-1 H> -.S ) I -.5 TO 0) <0 TO .5) (.5 TO 1> <1 TO 1.5) II.S TO 2) 10E 2) 


14 0 0 

B IB 

10 

4 

5 , 

7 

. 5 . 

.. 1 . .. . 

■ 

■ 

_H_L5vlQr»£AM*Ftf flM -1 in in -3/4.-.-7YK iG lSTU^ilFV^ 








7 4 5 2 

4 3 12 

1 1 

0 2 

0 2 

2 

1 




CUMULATIVE PROBABILITY 

n.noQoti o.Qfldftu _ 

(KDiHMJO .10714 

-42H57 

AA714 

A7Aq7 

.7B7H6 

-A92B6 

.90214 

1.00000 

, 









FUfiHT mSTOM'A’Hun.uF <l" 

•UtV.FHUM MEAN 

— 

— 


— 



: 

— 

EULGUL. JJJLzZX A -2 TO _r_L*h L 

_ -L-A jJ?___LU_rJ»l_.LnL 

-i3i_JLuib__L(JLJD_ 

_L0_j_a_»5_) 

.U5.J.Q_JL) 

Ll_Ul _1 a 5L _1 J_.S JtL.2L.16E ^2L 



14 0 0 

V 3 

5 

3 

5 

5 

1 

1 



HI STOOP 4^ ♦ F POM -1 TO -?/rt 

*-//M TO -J/4* 7/ B 

fO 1 STU.DEV* 








12 0 0 

1 0 1 J 

0 U 

0 3 

1 2 

1 

1 

i 



CUMULATIVE PROBABILITY 

0* 00000 0,000 00 

IK.IUJ 000. ^g.lP.b 

J 47 b GO 

1-53.L2S 

_?i3.250i> 

_j.7J9.L25_ 

19375.0 

*968.75 

1.00000 



I 



FL T (>H T MfcAM Mb7(l>H)AV 

SAMPLt SIZE 

bTU*Ut V * (DB) 

MEAN SIGZIOBIPK 

SAMPLE SIZE 

STD. DEV. (DB) 


4 15 lfi.17 

il / 

10*05 

4* 

<5 

• 

1 

99 

11.12 


5 1 (i/av ♦ 1 Slii utv \)n 
lY.i a 

FLIGHT hJST06KA”(N0*0F S1D.DF.V. 

MG2AV - 1 
1*4.95 
FPO^ HE AN 

Sfli UEV Ob S1G2PK ♦ 1 STL) (JEV D8 

19*95 

SIGZPK - I STD 
18.83 

1 i 
1 1 
1 i 
1 1 
1 1 
1 ( 
1 I 
1 1 
i 1 
1 1 
1 1 
J 1 
] 1 
j 1 
j 1 
1 1 
1 1 

s | j 
si i 





A-24 


HISTOGPAM.FROM -1 TO -7/B«-7 TO -3/At 7/8 TO 1 STO.nEV, 

3 3 3 *+ 4 3 611 7 5 


CUMULATIVE PHORAHlI.I TV 
p _ j_00 OjJO _ .0 3 * 1V - 


* il±W* 9 .^9060 _*49S73 ._6 58J.2 . B46 15 i?^872 *9^145 lj»00000^ 


R IGHT HIS rOGH/U* (NO. 01- SlU.UtV.I* JUM ME an 


■■ 

\ A \ -3^9 286^624440 

LL *ji u j_v c_Vj? L^-L — .J 6. 1Ji< 1 f 26jdfe3 •JLQl Q3 * $££££ Q — * 26S ZD 1 *Afl ftD Q_ _ 

_ __ _ _ „ A £>_* £*>* *!£>** ' 

FLIGHT ..Fill s I <»/ ( J*H ) A V SAMPLE Si/t b f U • DEV . (OH > MEAN SIGZ(Db)PK SAMPLE. SIZE STO.DEV.(UB) 

— - 7\7}i) Ad i .68 11.41 32 3*83 

M" 77T T ~ f TT” OR/ Ijh SIu/ aV - l Sll) UtV UH SIGZHK ♦ 1 STD DEV OH S1GZPK - 1 STD DEV DB 

7777 7774 ' nrm ^ruTss - 

- fi"7 gmT "h TsTt ig 7a“ in f77 nT — 71 17 .71777177 g7 “ 7 e an 

Fl.lGHl (If -2) (-2 To -i.b) fu -1) ("1 In -,b) (-.b liJ 0) (OTQ.b) («b JO 1) (1 TQ_L» 5) (1.5 T0._2J <UE 2) 

16 0 i i 4 it 7 7 _0 2 1 

HISTOGRAM. r" vo?* 1 -1 F<» _/VH __| u_ J/ ‘±AZT—Uli .J, il _L ^ XU^ut Vj*. 

? ij * 1 j'lii 1 IJ1212J1 1 


CUMULATIVE P'JiMAHli. I T * 
o . K'ofjA 


. i^su . iS^b « 3437S ^4 6^87 5 rMI^° •22£25_ 


1.00000 


**'*** ****** sbi ****** t>4 $ 1 ) b * T >* S * S S 


~| IGHI Hlsfng ^ r » ( f ,' r.'ir -> 1 ' / • I »h V . h ^ U . Ht AN 


hi stog^am . y /(jm -i r ii !li * JJJ_ 4? z_"ir Zi *1 JJL i_ J5.1U *1^ )Lr 

oi i > i 3 i i”J o 4 o l l J 0 


CIIMIM A r I VH WMhAMtl \ 1 T 
0 . UOlSMt . U \ \d-> 


1,00000 


r>**±5,* fc******s*i*bS*£^b%^_^^*S*_ 



A-25 


A-4 I DISTRIBUTIONS F ow SCLfCltli KUNb 

x . ' FLIGHT ?UN ^KaT* VuVaIuh) STD . HT V . SArtTblTt MEAN SIGP(OB> S7U.UEV. SAMTsTZE 

* 8 11 8“.7v" ” To ?OlF3 7T728 TO 

ii — _ _ _ _ _ 

RUN HISTOGRAM (isin.OF b fu.UtV . h MEAN ~~ ~ "* ' "* 

* 

lUM _l*S_[o J9L JQ_IP_i£L J.«b_J_0 JL jJ_LQ_i»S>_ a«5>_Tp. 2J. jtGE_2> 

§ 11 o '1 ± j 1 __2 A 2 1 1 

HISTOG^AN«FROm -1 H> ro -3 /h. 7/8 TO I STO»m;V. 

•7021 122OU20O1O1O 

CMMULA fl l7lT v * 

^iitSLOofi LLr-XX ily ixx Xii ._*oojjo . j>_oo_o •_b_&ooo ,_/ooo_o .aoooo ^qooo _ . 95ooo 1.00000 

Rt7N~HTsTot>7 <li7f'io 7\)F ~s t7jT1 TVcT*” SlSn 

FLIGHT uun ( L I - P ) <-P lu-l.s) (-1.5 fO -1) (-1 10-.5) WS Tu 0) (0 TU .S> (.5 \0 i) (1 TO 1.5) <1.5 TO. 2) <GE 2) 

ft. X7 0 1 2 4 3 A 2 2 ^ 0 

U X XQ_ zJXXi =^i/iL XX X X LOjtiih.'Lt 

111 1 0 11120*01001 

CUMULATI Vb P^iMAiiu 1 Tv ' _— 

Jlal'Mim aX<X>JL U^lWLQ *.3£XPX »XOfl<L0 ^MMJL ^SU_PiLQ UfiiULQJtt l*J)SLQS)SL 

ili iiiU'iiiii iijmiilti 


DISTRIBUTION* Hh M-irLftl/ RUNS 


FLIGHT 

KU*>J 

mP aN 1 in* { 1 H> ) 

s nxotv. 

SAM.blZt 

MEAN SIGP(DB) 

STD. DEV. 

SAM. SIZE 

8 

"”T4 

1 ^.0! 


*0 

1 f. 15 

9.09 

20 


PtIN H ISTCKWa * (NO. OF ST.Jiiitv.HlJ* Mfc AN 


XJX GbX A+ LLJ — (_L Lr 2X- XrX _L no X X-X - jJLL. IrJL _l U_ .=.01 X XtX _Tll X X XlXXX*5 X XXa_TiLX X XX XXL. 5J__ (1.5 TO 2) _iGE_Tl 

3 lit 0 I d £ H 5 0 2 1 1 

tu SLtunxai’ii. L^nr—ii _j±l _=• £/-c. _j ll Jtixji nxait y_, 

02 1 l 12100. J U <10000 

CUMUlATlVb |U ^ _ - — — — - “ 

g^MinUl O.vn ,!,:*» IXJLim aizMil i3bU0Q »bb 0 Qu __.rioooo ,80000 .pqoqq .9500q 1.00000 

Sl* ^ is ; 

RIIM HISTOG* v*(i|fi.nF srj.M-V.H.;-. Mt — ~ 


*- 


« 


iUGiXjJlJjt_lLL^XL"X.J.x-loLl-J Jl J5_W_-lLL-lJsi -OX_L^JJlJLl(L XX^?LX^?XXaXXXXXXbJ_U^5_X>_XlXii£X2J 

ft 1 ft [> I | A 6_1 3 1 J 0 1_ 



A-26 


H I STUfiR AN* *-*»<>■» - I Ml TO - * 7/rt IU I STU.lihV* 

llliPI2lO210 0 00l 

CumTTj"a ri vh" --~\T~ ~~ 

LirlL'± *J.iJ_00 il • jOu OU .60 QOJJ. ._7 500,0 . BOO 00 .95 00 0 .^500_0__ _ Jl* 00 000 

AlJb l^fa****^***** S1?> *** Afrfcfe 

niSTPlUuTIOMS F t »► SrfllliU »«m-* 

Flight FITnj oJJiJTi s7J7,77. sA,T.~i7t mean siopiuh) si u. dev* sam.size 

_ — — To HTTob 5727 Jo 

RUN HISTOukhm (h'n.ur Si UUl’ S»r<' , i> < i**t M'M 

FI-IfiHi nH (I.I-M (-f I -I .^) 1 - 1.1 Lii. -I) t-1 m z*2l Lth hi m <u II* .S> 1.^ ro li il TO 1,5) (1,5 TO 2) (6E..2) 

__H lb. I £_ 4_ 4 4 b 1 0 1 

HI SJ fKV « Af< * \ IM ~ i | ' r « ~ //^ *j_ *j/u t «■«>- _//_rt 1 O | S I O.l l)h. _V» 1 « 

1 ) •> ' <. 2 I 1 " * 2 0 0 1 13 0 

I 

CIIMUI A I 1 VI- 'll 1 \ ' 

. ZJ-3LL 1 1 5-1*112 j£P .6bj^0U ^9U_000 *_95p^0 *^000 1^0000 

_________ __ _ _ vy_* ^ VbV b»£b_S *** Sj> VM» *± Vbjfc 4* 

\ . ,-- 

PUN HlSTObkfM ’'*.•■** M - - ,r s -r / 'J 


fl_t_ght_ ih ij_ tj; i_ i_ -j* tj_> _rJ_)_ lr_L _±\L Jil 
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I. INTRODUCTION 


The variable amplitude of the sea surface waves constituting a sea wave 
profile can be recorded on a photographic film if the sky light providing the 
illumination of the sea surface varies continuously in a linear manner. 

The amplitude of the sea surface waves may be defined in terms of a 

( 3 ) 

potential energy spectrum of the wave surface. If one thinks in terms of an 
infinite summation of sinusoids defining the sea surface, each sinusoidal com- 
ponent must be characterized by a frequency and vector wave number. The energy 
spectrum of the sea surface will then, in general, depend upon surface position 
and wave number vectors and the time and frequency scalars. The photographic 
record of sea surface amplitude variations or slope variations is also a repre- 
sentation of the potential energy spectrum of the surface. 

The optical density variations on the photographic film, referred to later as 
the scene negative, are in a one-to-one correspondence with the sea surface dis- 
placement or slope,pr oviding the linear sky light illumination condition is satisfied. 
The optical density characterizes the variable light transmission property of the 
negative and is defined by the equation I = 1^10 D , where I is the light intensity 
transmitted through a medium of optical density D, and I is the incident light 
intensity. For a photographic film, the optical density D is related to the light 
energy which exposes the film by the linear relationship 

D = y log E (1) 


(3) Kinsman, Blair, "Wind Waves, Their Generation and Propagation on the 
Ocean," Prentice Hall, Englewood Cliffs, N.J., 1965. 
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for long enough exposure times. The quantity E is the exposing energy density 
and depends on the product. 


E 

where 

I 

T 

k 


= klr 


j oule s 

2 

m 


( 2 ) 


= exposing light intensity. 


watts 

2 


m 


- exposure time, sec, 

= relates to film sensitivity to light 

= the slope of the film characteristic curve (log E vs D) in the 
linear region 


A camera located above the surface of the sea and looking down in the 
direction of the water wave propagation will record a two dimensional optical 
density pattern - D(x, y). As indicated by equation (1), the numerical value of 
the density at each point on the scene negative will depend upon the light intensity 
exposure time, film gamma, and sensitivity. 

The light intensity (I)received by the camera lens will depend upon: 

(1) the geometry of the sea surface and camera position, (2) the manner in which 
the skylight intensity varies with zenith angle and azimuth angle, and (3) the polariza- 
tion and frequency sensitive water reflectivity. 

A. Linear Sky Assumption 

For a basic understanding of how the slope angle <J> of the sea wave profile 
is recorded on the scene film, the assumption of a linear sky is a good starting 
point. If the sky light reflected off the ocean surface varies linearly with zenith 
angle 0 (see Figure 1), then the light intensity that exposes the camera film I(x’) 
will have a functional dependence which is linearly related to the slope angle <|>(x) 
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or wave height T|(x). This important concept is illustrated below* Assume the 
sea wave system is such that no azimuth angle dependent light reflection is 
recorded. 


Z 



Figure 1* A One ^Dimensional View of Ray Geometry for 

Determining I(x) Under Calm Sea Condition* Distance 
x! is Measured on the Scene Negative. 


In Figure 1, the light intensity I(x ! ) will have the form 1^ + I^x 1 provided 
the sky illumination a(6) is linear with zenith angle 0* The light ray a will 

C* 

undergo specular reflection at x = 0 and expose the film at a point corresponding 
to say x' = 0. Light ray will reflect at x = x^ and expose the film at point 
x 1 = x’^ with greater light intensity. A measurement of optical density D(x‘) along 
the x' axis of the scene negative will be that shown in Figure 2 (microdensitometer 
tracing) 


B-5 



D(x’) ~l(x') 



Figure 2 . Scene Negative Optical Density Trace for Linear Sky Condition 
If the sea surface is now wind driven the camera illumination I(x') is of the 

form 

I(x') = Ij + I 2 (x* + c^ (x 1 ) ) (3) 

where <j>(x') is the slope angle and c^ is a constant relating angle to distance. 

Comparing Figures 1 and 3, it is seen that the tangent plane for specular 
reflection defined by slope angle <{>(o) i- n Figure 3 will provide illumination at the 
camera only by sky ray a .and not a (Snell's Law). In terms of the unperturbed 

D & 

sea state of Figure 1, this corresponds to specular reflection from point x^. 
Hence, the term c^4>(x') must be added to x’ to account for the increased light 
intensity exposing the film at x’ =0 (see Figure 4). 

A scene negative will thus contain an optical density pattern D(x’) 

D(x') = y log £l A + I 2 (x’ + c 1 <j)(x , ))| kr ( 4 ) 


B-6 





where if a is the amplitude of laser light incident on the scene negative, a 10 
o u 

is the transmitted laser light amplitude. 

B. Optical Computer 


-D(x* 




( 5 ) 



where £h is a calibration constant, and K is the vector wave number in radians 
per meter. F is the position vector. The vector wave number corresponding 
to the point u^, v^ specifies the wavelength, direction, and frequency of one 
sinusoidal component of the sea wave spectrum. 


C. Results of Calculation 


The energy spectrum is calculated from the optical density pattern D(u, v) 

provided by the Fourier transform negative. A two dimensional integration of 
-D(u, v)/v 

10 , the optical density transmission function, and evaluation of appro- 

priate calibration constants determines the energy spectrum. The calibration 
constants are calculated from microdensitometer tracings of the scene negative, 
transform negative, laser beam cross-section negative,Ronchi grating transform 

negative (fA product), and camera focal length and position. The microdensitometer 

-D/y 

basically measures the optical transmission property of a negative, 10 


A one -dimensional energy spectrum calculation for Raytheon Flight No. 6 

ocean wave data indicates that the peak spectral component of the wave system 

2 

is in the range of 0. 188-0. 226 ft -sec at a wave length of approximately 20 meters. 
Assuming 20 meters to represent the dominant wavelength of the system, the 
corresponding period is 3.6 seconds. Shore observation recorded a wave period 
of approximately 5 seconds for this flight. Also, for frequencies greater than 

two radians per second, the spectral amplitudes decay according to the well 

-5 

known (JO dependence where o> is the radian frequency. 
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II. SEA PHOTOGRAPH ANALYSIS 


— 13) {x^ y M 

The general two-dimension transmission function 10 v * 3 ' represented 
by the scene negative is now derived for the two-dimensional case. The analysis 
approach follows that given by Stilwell. 


Consider the illumination at the camera from light reflected from some 
arbitrary point on the sea (see Figure 6): 


1(A) = Ka{0) p(9) (6) 

where 

K = a constant and varies inversely as the square of distance 

from the point of reflection to the camera position in the sky. 

a(0) = the luminance of the sky and is a function of zenith angle, 0. 

p(0) = the reflection coefficient of the seawater and depends also on 

the zenith angle. 

A = camera angle 


In Figure 6, the position vector r is drawn perpendicular to a differential 
area in the focal plane of the camera. & is normal to tangent plane. 

Equation (6) implies that the illumination of the camera is governed only by 
light rays that lie in the plane defined by the unit vectors r and n. 

The change in camera illumination for a small perturbation of the surface 
normal unit vector n (small slope angle <f>) is 


da(0) 

d0 


dp(0) \ d0 
d0 J d<j> 


dl 

d<}> 


K 


p (0) + a(0) 


(7) 




de_ 

d<j> 


1 d .A a A A a 

(r * n) where r • n = cos 0 


sin 0 dc|> 


sin 0 d<|> 


[- 


y cos A + z sin A) • (x sin 4> sin + 


y sin cj) cos ip + z cos 4>) 


»1 



cos ip for 0 + <j> + A = 90“ em.d cos 4> = 1, sin <f> = 0 


(9) 


So 


dl_ 

d$ 


K 


da(6) 

d0 


p(e) + j cos ip 


Integrating equation (10) gives 


( 10 ) - 


1(9, ft « - K (^L p(9) + a (e) ) cos ft • * + I o (11) 

The constant of integration I may be looked upon as the illumination I for ip = 90° 
or <(> = 0. Equation (11) is simplified by defining 

k * p ( q ) + a ( 0 ) -^p-) = X ' = £ ( 0 ) . w 


.*. l{Q,ip, <f>) = I + 1* 4> cos ip (13) 

o 

Equation ( 13) is the two-dimensional illumination function originally derived by 
Stilwell. 

In terms of general coordinates on the sea surface. Equation (13) is changed 
to read 

I(x', y*) = I (x 1 , y') 1- I*(x»,y») ^(x’, y*) cos vji (14) 

o 

The illumination terms I (x, y), and I’(x, y) are assumed to be constants 

q 

in Stilwell 1 s optical analysis . This would imply that the sky light intensity a(0) 
(see Figure 6) changes linearly with zenith angle 0 and all light rays reflected off 

If this assumption is true, then 

= constant. (13) 


the sea surface see the same water reflectivity. 
I 1 defined by equation (12) would take the form 


r = K y • p(6) 


The validity of the assumption I (x, y) = constant may be examined by inspection 

o 

of the optical density plot of the scene negative. (In Figure 2, the average change 
in optical density is related to I (x‘). ) 
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III. OPTICAL ANALYSIS' 


A. Average Power Spectrum of Photograph 


In general, then, for small slope angles 4>(x, y) (cos tf> = 1, sin cj> = 0) 


i(x, y) = i (**» y‘> + y')4> ( x '» y') cos 4* 


A scene negative will record the above illumination function in the 


form 


D(x', y 1 ) = v log [k I(x>, y*)r ] . ( 17 ) 

The primed coordinates refer to distances measured on the scene film, kj 
is a constant related to the film sensitivity. Tjis the film exposure time. 

A laser beam transmitted through the scene negative, which may be thought 
of as a diffraction grating, will experience attenuation as shown in Figure 7. 


y* 



Figure 7. Transmission Through Scene Negative 
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The transmission function for light amplitude is 


T'(x»,y*) = 10" D/2(x, ’ y,) = 


( 18 ) 


where "Yjkjj refers to the scene film. To put equation equation (18) into a 
form useful for analysis, a binomial expansion is made. 


T'(x*, y«) 


V,/ 2 


y* ) 1- r(x', y' )4>(x', y») cos +] 1 

tWo***' y *^ 1 [1 f Ux'ly 1 ) 4>(xt ’ yt) cos 


-Vj / 2 


(19) 

( 20 ) 


-yj 2 


£ [TUfx-y-)] 1 [1 - -j- g(x', y* )+(«'. y') cos 4» ] 


where 


, , A r X 1 , y’) 

g(x', y 1 ) = 7 

I (x 1 , y‘) 


( 21 ) 


( 22 ) 


and the assumption is made that 


[g(x‘, y'^x 1 , y 1 ) cos ijj ] « 1 . 


(23) 


The transmitted light amplitude a^x', y 1 ) is passed through the lens 
shown in Figure 8, and exposes the film placed in the back focal plane. The 
light intensity registered by the film is the Fraunhofer diffraction pattern 
(Fourier transform) of the light amplitude of the front focal plane. Each 
optical density point on the exposed film is uniquely related to a spatial fre- 
quency composing the transmission function of the diffraction grating (scene 
plane). The value of D(u, v) indicates the relative intensities of the frequency 
components. 
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V 


film 


Figure 8 


negative 


The light amplitude a^/u, v) exposing the film is therefore 


a^u, v) = const 


f iK 

P p - ^- 2 ■( ux, ' , ' v y , ) 

\ \ a^x 1 , y 1 ) e dx ! dy* 


= const ^ ^ a(x®, y 1 ) T ! (x f , y 1 ) e 


■(ux^vy 1 ) 


dx 1 dy ! (25) 


TOO . 

pp -Yj/2 Y 

= const J J a(x', y') [t^jI^x*, y*)] [l - — g(x', y'Hte'* y') cos 


- — - — (ux'fvy*) 



The constant of proportionality is fA (see appendix). 


+00 

a T (u, v) = (f A) JJ a(x', y*) 10 


-D 1 (x , ,y')/2 


-oo 


iko 


(ux'+vy 1 ) 


dx' dy* 


[l g(x', y'J^x 1 , y*) cos 4 / ] 


(27) 


where A is the laser beam wavelength, f is focal length, k = Zir/\ - and 

— P) ( x t y') h — y /2 ^ 

jq 1 ’ y. j ^ x , ,yj 1 • The film in the u, v plane is sensitive to 

1 l'O 

\ I 2 

the light. intensity [ a (u, v)| ! So , the optical density D (u, v) is given by 


7 D ? (u,v)/ Y 

k 2 T 2 I a T (u, v) | =10 


(28) 


where k^ is the transform film sensitivity factor, r is the exposure time, and 


Y^ the gamma of this film. 

From equations (27) and (28), we have 
+oo 


jj e . ( 1 / ^ )(x , V)io -^'.y'^ ^, yi 

(ux'+vy*) 


2 
iko 


)4>( x '» y') cos i^] 

2 

I 

dx* dy* 


i (d 2 (u ’ v >/V \ 


10 


k_ t a 
2 2 o 


/ 


(fA) 


(29) 

(30) 


. 2 f 2 ,2 

For a Gaussian beam, a(x', y') = a e V X + y Now let (Stilwell 

o 

assumptions) 


-D 1 (x',y«)/2 

10 = constant 


(31) 


and 


A D <j> 

Yj g(x', y 1 ) = — — be constant => g(x‘, y') = constant. (32) 

° g 10 
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Hence 


(fM 


2 (D (u, v)/y~)+ D (x‘, y‘) 


< k 2 Vo > 


(10 


2* ' " & r~ " ' 4 


lo S io e 


2 — 2 

cos ip D 

4> 


(33) 


+00 


ik 

f* P /i / ^\/ — ~~ (ux'+vy 1 ) 

( ( “t 1 / 0 " )( x f y ) ./ i n f • , , , , 

\ \ e <M X > y 1 ) e dx' dy* 


-oo 


The Fourier integral of the constant term for a finite aperture goes like 

and will contribute, to the frequency spectrum at and near center of the 
transform (co = 0). In the squaring process, it has been assumed that terms 
multiplied by the Fourier transform of the constant term do not contribute to the 
desired power spectrum. 

The average power spectrum recorded by the transform photograph is * 


( sin x > 

k x j 


P(K , K ) = 
x y 


2ir 


•foo 


ik 


- 2 , .2 


e -tr (X ‘ ^(xSy*) 


(ux 1 +vy* ) 


-oo 


where 

K 


—■u = K 
f x 


K 


■ v = 


= K y 


dx* dy* 
(34) 

(35) . 

(36) 


and u = x’, v = y* for the geometry of Figure 8. 

Equations (35) and (36) relate density points (spectral amplitudes) in the trans- 
form to corresponding wave numbers. Equations (33) through (36) are the basic 
optical analysis formulas. 

To gain some insight into the meaning of equation (29), take as an example 

Y, 


1 - g (x 1 , y*) 4>( x *» y‘) = 1 - x' 


(37) 


* The factor 2 t r/<j^ is related to an elemental area in the transform plane 
and is derived in StilwelFs paper. 


B-17 






IV. SEA SPECTRA 


A. Basic Power Spectra 

The sea surface displacement T|(x, y, t) may be expressed as a linear super 

position of sinusoidal components (wave packet) provided that statistically 

J (4) 

independent, small amplitude components are assumed . 

n n r> +i{K x-KK y- COt ) 

r,(x, y, t) = JJj B(K x ,K yjW )e X 7 dK x dK y dC0 (40) 

R K co 
y x 


where the amplitude spectrum B(K ,K ,oj) is given by 


x 


i ^ c c c -uk xfK 7 _G ^ 

B(K x ,K y ,co) = (— ) J J jr,(x.y) e X 7 dx dy dt. (41) 

x y t 

The instantaneous power spectrum of the sea surface displacement may be 
derived from the above relations and is 


% (K X K y x> Y> ^ = ^ ) J j *1 ( x >y> t)T)(xf x o> y+y t) 


x y 


-i(K x +K y ) 
x o y o 

c dx dy 

o o 


(42) 


where the function 


H(x, x q) y, y Q t) ^ T) (x, y, t) TUx+x^ y+y Q > t) (43) 

is the autocorrelation function of the sea surface displacement. (The bar 
indicates mean value. ) Its inverse is 


H(x, x , y, y , t) 
o o 


K K co 
x y 


+i(K x +K y ) 

- x ° 7 ° dK dK to. 

x y 

(44) 


For a homogeneous sea (one independent of the reference position for meas- 
urement) and stationary wave process, we have, dr opping the x, y and t dependence. 


r n n H(K x fK y ) 

H <VV = ® ,K x’ K y“ |e X ° y ° d W“ 

K K co 
x y 


( 45 ) 


• (4) Ibid (3) 
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dK dK 
x y 


(46) 


where 


- tf 


$ (K , K ) e 

x y 


i(K x fK y ) 
x o y o 


K K 
x y 


$ (K , K ) = \ ® (K , K , w) dw (47) 

x .y J X y 

co 

(K , K ) is the two-dimensional wave number spectrum. The function 
x y 

$ (K K , oo ) is the three-dimensional spectrum for a spatially homogeneous 
n x y 

and temporarily stationary wave field. Kinsman states that 1 ^ (K^, w ) 
represents the waves generated on an ocean of infinite extent over which 
a statistically uniform wind has been and will be blowing forever 11 . “More 
reasonably, it isn*t a bad model for the interior of a three- or four-day 
North Atlantic storm covering an area of 500 by 500 nautical miles. n 


The spectrum that is easily measured in practice (Spar Buoy) is the 
one -dimensional frequency spectrum given by 


•<«> = Jj 1 *(K x , K , «) dK x dK y (48) 

K K 
x y 

It gives the contribution to the potential energy of the wave coming from 
each frequency co, irrespective of the vector wave numbers associated with 
that frequency. The frequency spectrum gives no information about which 
directions the various waves come from. (They add up at a point to define 
the potential energy of the composite wave. ) 

A typical frequency spectrum is shown in Figure 11. 

\ m 



Figure 11 
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Thus, for the ocean wave process, there are three basic power spectra: 

(1) The three-dimensional spectrum -$ (K x , Ky w) 

(2) The two-dimensional wave -number spectrum - $ (k.,k y ). 

(3) The one -dimensional frequency spectrum - 

B. One -Dimensional Frequency Spectrum 

A relationship between the wave number spectrum and frequency spectrum 
• for infinitesimal waves which have a unique connection between wave 
d frequency. For small amplitude, deep water waves, the relationship 


, . . 2cd 
$,(cd) -■ — — 
ly 2 

g 


2tt 

f 

0 


^ (K, ip) df 


(49) 


where g is the gravitational constant, if/ is the azimuth angle (see Figure 2) and 

2 „„ „ „ meters radians / Kn \ 

cd = gK where g = 9. 8 — , K = . ( bU ) 


sec‘ 


meter 


Kinsman points out 11 that finite amplitude waves do not have a unique relationship 
between wave number and frequency but yet the application of equation (49) seems 
to give very good results. 11 

For each wave vector K = x + y there is a energy given by 


-Jk 2 +': 


$ (K.K ). , All wave vectors of the same length VK +K and corresponding 

Tl’ x* y x y * 

2 

frequency CD = VgK lie on a circle of radius | K| = ” in the K , K plane (see 
Figure 1? K 
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^g^ATIQNSHlP BETWEEN SEA SPECTRUM AND PHOTOGRAPHIC SPECTRUM 



Figure 13. Geometry for Calculating Relationship Between 
$ (K , K ) and *AK , K ) 

T| x y q> x y 


The wave height spectrum 
as follows: 


is related to the slope spectrum 3? (to) 


For a wave propagating in the r direction in x-y plane (see Figure 13), 


8TL (x, y) 
9r 


tan (j> (x, y) 


<i>(x, y) 


(51) 
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and x 


r sin 4> 


so 


&J. (x, y) 


9x 


= 4> sin i}/ <v 


(52) 


For a stationary wave process, equation (40) gives 


Tit* 


• y,= i y v“ )e 


i(K x + K y) 


K K CD 
x y 


dK dK dcD 
x y 


(53) 


= iK x ^^ B(K^K ,cd) e x 7 dK dK (to 


x y 


(54) 


= iK T) (x, y) 


= <t>(x, y) sin f 


or 


iK 

' Mx ' y) = TIT# 11 * x< y * ' 


(55) 


<}>(x, y) = iKT|(x, y) . 
Therefore, equation (41) 

2 


l r r "UK X + K y) 

B(K^,K^,cd) = (— ) ^ ^ T)(x, y)e X dx dy 


K K 
x y 


becomes 


(56) 


(57) 


B(K x ' v m) = ^ ( K> J I +<*• y) 


i(,K x + K y) 
x y 


dx dy 


(58) 


and the power spectrum given by 
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2 


%(K , K ,©) 
! I x y 





-i(K x + K y ' ) 

(x, y)Tl(xH-x , y+y ) e dx x dy - y ° 
o o o o 


becomes 


( 59 ) 


$ (K , K , to) = (— r 
T1 x y 2 tt 


K 




x y 


(x, y)4> (x + x , y + y ) 
o o 


-i(K x + 
x o 
dx dy 

n J r 


K y ) 
y o 


(60) 


or 


K 2 $ (K , K , to) = (K , K , w) 
T1 x y 4 x y 

and 

K 2 ^ (K , K ) = * (K , K ) ) 

T) x y 4> x y ; 

(both for small <j>(x, y) . 


( 61 ) 


(62) 


The power spectrum defined by the equations (33 ) and (34) is proportional 

to * (K , K ): 

9 x y 


*< K *' K y »‘[{t) 


(fM 


(Ja^(u, v)/y 2 )+D 1 


(K- T a ) 
2 2 o ■ 


- - - , 2 a log ]0 e -| 
10 .4 sec <j> . — 


D 


'4> 


(63) 




The constant of proportionality is ( ~ y ■ ) where 6 , A are scale trans- 

W 3 I x y 

formation factors relating the sea coordinates x, y to tbe transform coordinates x 1 ’, y 


r (5) 


(5) Ibid (,1). 
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VI. CALCULATION PROCEDURE 


The calculation of wavelength is based upon Equation (5) which relates 
various points in the Fourier transform to corresponding vector wave numbers. 
The fA product is determined from either the optical bench or the Ronchi grating 
transform, f is the focal length of the lens and A is the wavelength of the laser 
beam. Knowing the focal length and position of the camera used in taking the 
scene picture, a scale factor relating a distance on the ocean surface to a dis- 
tance on the transform can be determined. A calibration equation relating a 

(3) 

distance in the transform x to a corresponding wavelength A^ is therefore 
determined. (Propagation along the orthogonal axis has been ignored. ) 

The computation of the energy spectrum function $(to) is based upon 
Equations (49, (62) and (63). Required for this evaluation are the calibration 
equation and microdensitometer tracings of: 

(1) laser beam cross-section 

(2) scene negative 

(3) transform negative or isodensitometer of transform negative. 

The transform density tracing shows rapid oscillations with frequency. An 
averaging process has been applied to the transform for a reasonable estimate 
of the spectrum function. 

(A) Wavelength Calibration 

The vector wave numbers and frequencies comprising the sea state energy 
spectrum 3? (K , K ) are calculated from the relationship between vector positions 

$ * y 

in the transform or isodensity tracing of the transform and wave number K. 

A typical Fourier transform and its corresponding isodensity tracing is shown 
in Figures 14 - 17 for f = 0 and f = 90° . The major axis of the pattern indicate 
the general direction of wave propagation. 
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Figure 16. Fourier Transform for Flight No. 6, ^ = 0 
Camera located in y - z plane 



Figure 17. Fourier Transform for Flight No. 6^ f =0‘ 
Camera located in x - z plane 
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The following relations apply for all points in the transforms: 


and 





K 



( 64 ) 


(65 f 


f is the focal length, of the lens and A is the wavelength of the monochromatic, 
collimated laser beam. 

Alternatively, the fA product can be determined by analyzing the far field 
pattern (transform) of the Ronchi grating. The basic equation is 


A 


Af 

*b = 2b 


( 66 ) 


where 2b is the distance between grating lines and Ax^ is the distance between 
two maxima in the transform recording. (See Figurel9 ). 

The first step in the wavelength calibration procedure is to determine 
the scale factors relating a distance on the sea surface Ax, Ay to a distance in the 
transform or isodensity tracing of the transform, and a wave numberE^, of 
the sea state to a wave number in the transform. 

In Figure 18 is shown the camera position relative to the ocean surface (ij/=0). 
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Figure 19. Ron chi Grating Isodensity Tracing for Measuring 
fA. f7\ = 2b Ax-h and Line Spacing is 2 lines/mm 






Let 


Ax 


n incremental distance on sea surface 


and 


= incremental distance on a 9 inch film in camera (scene negative) 

ss incremental' distance on a 35 mm film (scene negative) 

- incremental distance on isodensity tracing of transform 

K K ^ K K ^ are the corresponding wave numbers, 

x* x 9 x 9 x 

It will be assumed that the power spectrum is highly directional along the 


x axis and therefore all K =0. 


Now, for flight No. 6 and follow-on film processing work. 


Ax 


( 1 ) 


Ax cos 45° 


focal length of camera 
range 


1 

1500 VT 


so 



Ax 

3000 


(67) 


The 9 inch film is now reduced to 35 mm film: 


Ax 


( 2 ) 


Ax 


( 1 ) 


1 

14.4 


( 68 ) 


The isodensity tracing of the transform is enlarged 10 times over the 35 mm 
transform. 
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10 


Ax 

Ax 


( 3 ) 

( 2 ) 


The scale factor 6 is 

x 


1 


P 


x 



Ax 

Ax 


( 2 ) 

U> 


Ax 

Ax 


(3) 

( 2 ) 


or 


(69) 


Ax 


(3) = P x = 4320 


Ax 


and since wave number K is dimensionally 


radians 
unit length 


K 


(3) 


K 


= 4320. 


x 


The position -wave number calibration equation is 

K (3> 

Ax ( 3) = (fA) X 


2ir 


, . 4320 v 

= (2b Axb) K . 


27T 


x 


(70) 


(71) 


(72) 


The line spacing Ax^ is measured on the Ron chi grating isodensity transform 
(see Figure 19^ as approximately 0.50 inches. The grating is specified as 
having 2 lines per mm, so, 2b =0.5 mm. The result is 



27.4 

A 

x 


-2 2 
x 10 meters 


Calibration Equation 


(73) 
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(B) Frequency Spectrum Calculation - $ (co) 

- • .. . 


The two dimensional wave number spectrum in terms of the photographic 
parameters is 


P 8 / , (D /v >+ D, 

, k , k ) = -3 y- m 2 /v 2- 1 

♦ x y 2 2 L 2 10 

a tt \k r a 
'220 


sec f , 2 e 

2 108 10 (74) 

D , 


and using 


* (K , K ) = K * (K , K ) 
4> x y I] x y 


, , . 2cd 
§ (to) = — — 

T1 2 

g 


I V K x’ K y W 


we have 


vH-feKf-) 

l « wn « ,k 2 V 0 


1 2 

10 4 log e 


(D 2 (u,v))/y 2 2 


Equation (77) is the basic equation for the energy spectrum evaluation. The 
following numerical values apply: 

(1) 6 = 4320, (3 =1 

(6 =1 because a one -dimensional spectrum is calculated) 


2 £0 

(2) K = — - 


1 00 

2 2 
K g 


where 


00 = Kg for small amplitude, deep water waves. 


B-35 








(3) 


a = gaussian beamwidth of laser beam. 

From the one -dimensional density tracing of the film negative of 
the laser beam exposure (see Figure 2^) (which is not gaussian! o 
is approximately 0.875 inches. 

On the 35 mm negative, this corresponds to 
-4 

cr = 22. 2 x 10 meters . 

2 

(4) f 7\ = 0. 65 mm (Ronchi grating data) (see Figure 19) . 

(5) The factor k t a ^ is equivalent to the term 10 ^\in which D/y is 

2 2 o 

the optical density arising by exposing film in the scene plane to 

2 

laser intensity a . From the gaussian beam density trace of 
o 

Figure 20, the peak value of D is approximately 0. 62 . 


Z 0 6Z 

k r a = 10 =4.16 (Assume y = 1 in calculations) 

2 2 o 

(6) The terra 10 is related to the I q (x, y) of the illumination function 
I(x, y) and corresponds to the average optical density of the scene 
negative. Figure 21 is a one -dimensional density trace of the scene 
negative used for this example. 

_ D 1 • ,«0.74 _ _ 

10 =10 =5.5 

(7) The term 


_ _vl'fe',y 1 

D <|, I(x',y') og 10 

o 

is related to the slope of the average optical density trace of the 
scene negative (see Figure 21 ). 
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If we assume (Stilwell linear approximation) that 



dD dx 
dx d<J> 


then 


D ± = 0.47. 
9 


The frequency spectrum, $(cd), is therefore given as 

2 


9 (CD) = 


(2)(4320)(0. 65 mm 2 ) (5.5) (4) log 2 Q e 
(22.2xl0 -4 ) 2 tt 2 cd(4.16) (0.47) 2 


2tt 

s 

0 


D 2 (u, v)/y 2 2 

10 sec ^ 


meter - sec 
radian 


(78) 


2ir 


3.69 

03 


x 10 


-3 


I 


D 2 (u, v)/v 2 2 2 

10 sec ^d^FT -sec. 


(79) 


CD = V K g = 
x 


From the calibration equation (78), 


(9. 8)(6. 28) Ax 


(3) 


27.4 x 10 


-2 


= 15 


forA* (3) in meters. 


.*. CD = 15 




(80) 

(81) 
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Equations (79) and (81) and the density tracing of the Fourier transform 
shown in Figures 22 and Z 3 were used to calculate both the one dimensional 
energy spectrum (K = 0) and wavelength calibration scale for the isodensi- 
tometers shown in Figures 14 and 15 • A 21 point calculation of the energy 
spectrum was based upon the smooth curve construction shown in Figure 23, 
The smooth curve values and equation (78) resulted in the energy spectrum 
plot of Figure 24 (positive wave numbers only). An experimentally determined 
energy spectrum for a different Raytheon flight is shown in Figure 26. 
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(FEET -SEC) 


EXPERIMENTAL DATA: 

SIGNIFICANT 

WAVE HEIGHT IS 

4-6 FT 

WIND 16 KNOTS 

WAVE PERIOD 7 SEC. 



FREQUENCY IN Hz 


Figure 26. Laser Profilometer Data. The average optical density 
that appears in the scene negative density tracing and 
manifests, itself as a dc term (cJ6=0) in the Fourier trans 
form does not appear in Figure 26 





VII. CONCLUSIONS 


The calculations presented in the preceding section outline the basic 

r 

computational technique developed by Stilwell for analyzing a sea energy spectrum 
from photographs. The particular sea photograph analyzed was taken on Raytheon 
Flight No. 6, and the resultant calibration photographs and transforms were done 
by D. Stilwell at Naval Research Laboratories. The microdensitometer tracing 
work was done at Raytheon, Autometrics Division, under the supervision 
of S. Hendrickson. 

The method of energy spectrum analysis of sea waves as outlined by 
D. Stilwell has been applied to Raytheon Flight No. 6 ocean wave data. 

The calculated wavelength calibration equation and one -dimensional 

energy spectrum §(cd) indicates that the peak spectral component of the wave 

... ft2 -seconds . , - 

system is in the range 0. 188 - 0. 226 radian at a wavelen g th o£ approx- 

imately 20 meters. (See Figure 24 ). Assuming 20 meters to represent the 
dominant wavelength of the ocean wave system (gravity waves) of the photograph, 
the corresponding period is 3. 6 seconds. 

As pointed out by O. M. Phillips^, the consideration of wind generated 
waves when the duration and fetch of the wind are large, justifies the considera- 
tion of an equilibrium range in the energy spectrum for large values of 

-5 

frequency (D. The spectrum obeys the relationship «&(<&) =cdi , where c is a 
constant. Application of this asymtotic rule to the calculated spectrum is shown 
in Figure 24. 


(6) O.M. Phillips, The Equilibrium Range in the Spectrum of Wind 
Generated Waves, J. Fluid Me ch., 4, 426-434, 1958. 
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Shore observation data recorded a wave period of approximately 
5 seconds and a significant wave height within the range of 3-5 feet. The time 
average wind velocity was 8 knots with 16 knot gusts. Aircraft observation 
reported that the seas for this flight, unlike that of previous flights, appeared 
somewhat confused. Even though there appeared to be a well defined wave 
propagation direction, other wave directions were observed. 

Experimental measurements of the energy spectrum were made with a 
laser profilometer. Results for a different flight are shown in Figure 26. 
Raytheon Flight No. 6 data was 'not available. 

The basic theory of sea wave energy spectrum evaluation from photo- 
graphic data as developed by D. Stilwell rests on assumptions that require a 
good deal of physical insight into ocean wave properties. It appears that the 
really fundamental aspect of the entire process is the requirement for under- 
standing what sky- sea-state is necessary to obtain scene photography that 
indicate an accurate measure of the instantaneous power spectrum of the 
-local sea state. The sky- sea- state referred to above should be related to such 
parameters as: 

(1) Type of wave system (capillary, gravity, infragravity, etc. ) 

(2) Position of sun relative to camera position. 

(3) Optimum sky condition (uniform overcast, partial overcast, etc. ) 

(4) ' Wind velocity and direction (related to type of wave system). 

(5) Wave profile function rj (x, y) (white caps, sinusoidal, etc. ) 

(6) Reflectivity characteristics of ocean surface. 

Condition (3) related to the average sky illumination with zenith angle 
can be determined from examination of the optical density tracing of the scene 



negative. A linear change in density along the x or y axis of the scene 
negative constitutes what Stilwell refers to as a "good picture. " If the linear 
sky condition is not met, the spectrum evaluation does not represent 
the power spectrum of the two dimensional slope angle cj)(x, y) characterizing 
the sea wave profile. The illumination function I(x, y,A)* that characterizes 
a scene photograph is really the most important function in the entire optical 
analysis. Its mathematical form must be understood before detailed analysis 
of the sea wave system energy spectrum can be made. With I{x, y. A) known, 
one could then give an accurate interpretation of the scene negative density 
tracing as far as identifying the local sky- sea- state. Consequently, the 
Fourier Integral of the height function ij (x, y) or slope angle function 4>(x, y) 
could be clearly identified in the transform isodensity plot, or in other words, 
the region in K- space that corresponds to the energy spectrum of the slope 
angle function would be known. 

Several observations based upon the analysis and calculation are given 
below: 

(1) The scale factors that relate distances and wavelengths of the 

sea surface to distances and wavelengths in the transform are 

two dimensional functions, p (x, y), (3 (x, y). 

x y 

(2) The transform isodensity tracing shows an asymmetric spectrum 
plot and one that varies rapidly with wave number. It is suggested 
by Cox, Monk, and Stilwell that an averaging process be applied 
to this density function so that the general trend in the spectrum 
behavior be ascertained. It is this trend which is subjected to 
optical analysis in order to evaluate real sea energy spectra, 

* The illumination function represents the light intensity reflected into camera 
from sea surface. 
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( 3 ) 


The collimated monochromatic laser beam is not truly Gaussian. 

As indicated in the optical analysis a Gaussian beam is assumed and 
is used to calculate the Fourier transform of the scene plane trans- 
mission function T(x', u') modified by the beam function 


-2 2 2 
exp a (x 1 + y' ) • 

The Fourier transform is 


{FT T' (x 1 , y') e 


-l/a 2 (x' 2 + y' 2 ). 


One is actually performing the Fourier transform of the slope 
angle function multiplied by the beam shape function, re suiting 
in a distortion of the slope angle power spectrum recorded in the 
transform plane. 

(4) The theory as developed by Stilwell does not allow the energy 
spectrum evaluation for all azimuth angles 


2tt 

${.<*>) ~ ^ 10 A2 sec 2 if/ dtp 

0 

2 7T 

since sec was f — . The transforms generally 

indicate a non-zero wave number (K) near the center of 
the transform. 
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(5) The density tracing of the scene negative could be enlarged and a 
detailed analysis of slope function made. A mathematical description 
of this function could then be employed in a digital computer, Fourier 
transform calculation. 

(6) The question of accuracy of results in the optical analyses of 
energy spectrum remains to be answered. For example: 


(a) finite granularity of film 

(b) grain characteristics of film 

(c) camera movement and wave movement during exposure time 

(d) the number of photographs necessary of a particular sea 
state to completely characterize its energy spectrum. 

(e) aperture size of microdensitometer slit used for plotting the 
transform 

(f) lens field of view as related to calculating the scale factors. 

(g) Statistical limitations of calculated energy spectrum §?(w). 
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VIII. RE COMMENDATIONS FOR FUTURE WORK 


(1) Extend mathematical analysis of entire Stilwell process, i. e. , 
linear sky assumption, non- Gaussian laser beam effect on 
calculated spectrum, camera field of view as related to scale 
factors, etc. 

(2) Develop theoretical understanding of the statistical limitations 
of calculated energy spectrum. 

(3) Exa min e optical bench assemblies that would provide efficient 
calibration of scene negatives and transforms. Use photoelectric 
cell or microdensitometer for calibration equipment. 

(4) Optically process existing Raytheon ocean wave data based on 
part (3), and present knowledge of optical evaluation technique. 
Select those scene negatives that have a linear variation of 
optical density and experimental data (laser profilometer and shore 
data). Compare in' detail, experimental data with optical analysis 
data for a large number of scene negatives. 

(5) Study the Stilwell process for application to radar scattering 

measurements. What, if any, correlation exists between energy 
snfictrum and the corresuondine radar cross-section a ? 



IX. APPENDIX 

A. Relationship of Optical Power in Scene Plane to Optical Power in 
Transform Plane ^ 


In the back focal plane of the lens (transform plane), the light amplitude 
v) is the Fou; 
of the scene plane 


a^(u, v) is the Fourier Transform of a^x 1 , y ! ), the light amplitude distribution 


’too 


a T ( u » v ) = (const 


•> II 


-{iK /f)(ux'tvy') 


a’fx', y 1 ) e 


dx 1 dy 1 . 


-oo 


The corresponding light intensity is 

2 


u T (u, v) = |a T (u, v)| 


= c 


too 

II 


-(iK /fXux'fvy* ) 


a'(x', y 1 ) e 


dx' dy* 


-oo 


where c is the constant to be determined. 


The light amplitude distribution a^x', y’) is 


a*(x',y') 


-l/(r 2 {x» 2 fy' 2 ) 

a o 6 T'(x', y*) 

laser beam Transmission 

(Gaussian) function of 

scene negative 


In the absence of the scene negative (T^x 1 , y 1 ) -*■ l) and evoking the require- 
ment of equal optical power flow in the input and output focal planes, 
we have 


2ir oo 

II 
0 0 


2 -2/r 2 r' 2 . . , , 

a e ^ r'dr’d©' = 

o 


2 

c 


^ , , 2, ,2, ,2 -(iK /fHrf'V 2 ) 

e -l/<r (x- ty* ) e o ^ dy , 


II 11% 

0 0 -oo 


r" dr" d0" 


output plane powe r density 


B-52 



where the independent variables are cylindrical coordinates about the optic 
axis. The single prime refers to the scene plane and the double prime 
refers to the transform plane. 


The power density 

foe _ _ -(iK /f)(x" 2 fy" 2 ) 

-l/o- (x* fy 1 ) e ° dx* dy' 

-00 


2 


2 2 4 
a i ff 
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2 2 2 
tt o- r" 8 
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2 2 2 

r M = x M fy" , f = focal length, and X = beam wavelength in angstroms = 2nl¥.^. 
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APPENDIX C 


THE FINE STRUCTURE OF THE SEA SURFACE 
Dr, B. Kinsman** 

Radar altimetry over the ocean confronts difficult problems -- not the least 
being how to describe the target. In a search for an answer it has turned to the 
oceanographers and has found little comfort. Waves are a primary response of 
the ocean to the wind -- and extended regions of calm seldom exist for very long. 
The waves of the sea are highly irregular, they run rapidly across the surface, 
and they constantly change their forms. It is only since the second World War 
that oceanographers have made a start at understanding waves as they exist on 
the ocean and have begun to create usable forecasting methods. Understandably, 
their attention has been focused on n pr actical 1 1 problems, which in this case means 
the larger waves -- those big enough to interfere with amphibious landings, 
refueling, and air -sea rescues. Now radar altimetry tells us that the small 
waves, those with lengths of the order of 1 to 10 centimeters, are of intense 
"practical 11 interest. They are the ones which have the greatest effect on radar. 

For waves on the water surface with lengths in the range 0{lcm) to 0(10 cm) 
both gravity and surface tension may be important restoring forces. Waves with 
lengths greater than 5.47 cm are primarily gravity waves; the effect of the surface 
tension on their phase speeds is less than 5% of the effect of gravity. Waves with 
lengths less than 0„54 cm are capillary waves; gravity contributes less than 5% 
to their phase speeds. For waves with lengths between 0. 54 cm and 5.47 cm both 
gravity and surface tension are important. 

A great deal is known about the physics of these small waves. The 19th 
century hydrodyn ami cists did not neglect them. Lamb ! s 1 ■hydrodynamics, 11 for 
example, records this extensive body of results. Our theory even provides us in 

* Oceanographer 

Chesapeake Bay Institute 
John Hopkins University 



Crapper 1 s wave that very rare commodity: an exact solution to a nonlinear 

problem. Further, theory has been abundantly confirmed by experiment. Small 
waves are relatively easy to house in a laboratory. 

What clues, then, do theory and experiment offer us? The linear first- 
order, small -amplitude theory is satisfied by a sinusoidal wave and, because of 
linearity, by sums of sinusoids. The theory tells us that each component wave 
will run at a phase speed or celerity, c, 

c 2 = (pg/k + kT)/ (p coth kh) 

where 

p = the density of the water, 

T = its surface tension, and 

h = its depth; 

while 

g = the acceleration of gravity and 

k = the radian wave number which is 2ir divided by the wave 
length, L. 

When the water depth h is greater than half the wave length, the waves are 
in "deep water, " coth kh a 1, and the phase speed is very closely approximated 

by 

c 2 « (g/k) + <T/p)k . (C-l) 

For waves on the range 0 (1 cm) to 0 (10 cm), even if we interpret 0 (10 cm) 
as 50 cm, any depth h> 25 cm will be deep water so that (1) applies to most 
of the ocean. 
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Equation C-l has a number of interesting features. The effects of gravity 
and of surface tension on the phase speed are additive. It is on the basis of 
(I) that a simple calculation yields the wave lengths for 5% residuals previously 
given. Thus with an error of no more than 5% in the phase speed c we may use 
for gravity waves 

c » g/k whenever L> 5.47 era' (C-2) 

and for capillary waves 

« {T/p)k whenever L< 0.54 cm. (C-3) 

If 0.54 cm < L < 5.47 cm the waves may be called ripples and the full approx- 
imation, equation (C-l), is required. 

The phase speed for gravity waves is a monotonic increasing function of 
the wave length, c ~ L. For capillary waves it is a monotonic decreasing function 
of wave length, c ~ l/L. Clearly, within the ripple band there must be a minimum 
phase speed which can be found by searching equation (C-l) for extrema. One 

1/4 

finds that no wave can travel slower then c = (4g T/p) which corresponds 

1/2 m -2 -1 
to a wave number k = (gp/T) ' . Using g = 980 cm sec , T = 74dynescm , 

_o m -1 

and p = 1 g cm one finds a slowest speed of c^ = 23.2 cm sec associated 
with the wave of length L = (2TF)/k — 1.73 cm. Waves either shorter or longer 

than 1. 73 cm travel faster. From Figure C-l it is clear that there is always a 
pair of these waves, one shorter than 1.73 cm and one longer, which travel at a 
common speed. For such a pair the shortwave will appear as a fixed roughness 
on the longer wave. There will be other short waves whose speeds are near 
those of the long wave which will appear as a slowly changing roughness on the 
long wave. Such phenomena do occur in a wind driven sea and are clearly 
apparent to the eye. That they are precisely described by the linear theory may 
be doubted since the linear theory was developed for small -amplitude (math- 
ematically infinitesimal) waves. If you can see it, its not small -amplitude. 
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When the small -amplitude restriction is removed theory tells us that the 
sinusoids become deformed. For gravity waves one solution is an inverted 
prolated trochoid. Such curves can be generated by tracing the locus of a point 
on the radius o± a circle which is being rolled along the under side of a straight 
fine. Figure C-2 shows a trochoid which is improbable as an ocean wave because 
6 and k were selected to be 1, Ocean waves theoretically never become that 
n pointy fr , the maximum angle just before they break being 120° . If you look at 
the ocean or at plunger generated gravity waves you will see that the trochoid, 
while too regular to be a good description of the sea surface, is still much more 
realistic than, the sinusoid. The sinusoid always has smoothly rounded crests 
while the crests of the trochoid have an edge. A sinusoid is symmetric about 
a level halfway between crest and trough. You can’t tell whether it is upside 
down or right side up. The trochoid has troughs which are long and flattened 
and crests which are short and angular so that mean water level is not half way 
between crest and trough. Further, its wave steepness, i.e^, the ratio of its 
height H to its length L, 6 s H/L has some sizable value. It is not infinitesimal 
as 6 must be for the small -amplitude sinusoid. Thus the trochoidal wave can be 
visible to the eye. The effect of this finite steepness is to make the trochoidal 
wave run slightly faster than the speed given by equation (021 as shown in 
Figure C-3. 

For capillary waves the exact solution of the nonlinear problem which 
makes no obeisance to the small amplitude assumption yields a water surface 
shown in Figure C-4. As drawn, the uppermost line is the water surface for a 
capillary of maximum steepness, 6 = 0.73 in the sense that locally, near the 
troughs, the slopes have become vertical and the surface recurved so that the 
surface is about to coalesce. The lower lines are streamlines of the flow for 
this extreme capillary. However, if a capillary has a steepness of 6 = 0.53 the 
second line is the shape of the wave and the lower lines remain streamlines. So 
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also for the others. What you see in Figure C-4 may seem startling since most 
of the waves you have looked at have probably been gravity waves that resembled 
Figure C-2. Here it is the crests which are broader and rounded while it is the 
troughs which are narrower and more pointed. However, in any process governed 
by surface tension, as, for instance, the formation of droplets on a non-wetting 
surface, the shape taken is always as near spherical as the other forces acting 
permit. Strange as they may seem these wave forms have been shown to exist 
in laboratory experiments. I think the bearing of such wave forms on radar 
should be obvious. The larger waves of the sea are like trochoids whose maxi- 
mum steepness is 6 = 1/7 = 0. 143. Observed values even in young seas which 
tend to be steepest seldom show values greater than 6 = 0. 12. Even locally they 
don’t become very steep. However, these longer waves may be and frequently 
are covered with short waves. If the short waves are capillaries then the overall 
slope may reach 6 = 0. 73 and the surface can locally approach the perpendicular. 



Figure C-4. Streamlines for Crapper’s Wave 
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As with the trochoid al wave there is a correction to the phase speed given 
by equation (C-3J which involves the finite steepness as shown in Figure C-5. In 
this case the effect of the finite amplitude correction is to slow the wave a bit. 
The combined effect on Figure C-l of the corrections shown in Figures C-3 and 
C-5 is to raise the curve a bit when L/L < 1. It may also change the value of 
k used for scaling but only very very slightly* 



Figure C-5* Phase Speed Correction Factor for Pure Capillary 
Waves as a Function of Wave Steepness 

Even casual observation of wind driven seas shows that a fine structure of 
small waves riding the longer, larger waves is the rule rather than the exception 
in reasonably strong winds* As recently as 14 years ago we had no remotely 
satisfying deduction from physical principles to tell us how the wind might make 
these small waves; those with lengths less than 5*47 cm for which surface tension 
must be reckoned with* Certainly if surface tension is important then viscous 
dissipation will probably be important and these waves will have to be continually 
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regenerated by the wind. Further, when- a wind strong enough to make waves 
flows over a calm water surface the first thing that always happens is the forma- 
tion of waves with lengths of a centimeter or so running, not with the wind, but 
symmetrically at an angle to both sides of the wind. They form a typical diamond 
shaped or rhomboidal pattern on the water surface. In 1957 Phillips offered the 
first mechanism deduced from hydrodynamic principle that could account for some 
of what we see. This is hardly the place to attempt an account of the Phillips 
generation theory. Suffice it to say that the generation of small waves by his 
process depends critically on a resonance between the waves and the pressure 
field on the water surface created and swept over the water surface by the 
turbulent wind. The resonance is selective for the shortwave lengths and the 
off -wind angles of wave propagation. However, if the turbulent wind does not 
contain appreciable energy in the small pressure structure it can be strong and 
still create no small waves. The locus of resonantly excited wave numbers in 
wave -number space is shown in Figure C-6 which may puzzle you a bit. Actually 
it has been borrowed and contains more information than we need here. The 
wave-number component axis is lined up with the mean wind so that the k^-axis 
is cross-wind. This is only half the plot. There is a mirror image in the 
k > 0 k < 0 quadrant. Further, it is assumed that the pressure field on the 

J- Lt 

water surface expressed in wave -number space has no appreciable energy (for 
whatever reason -- viscous dissipation in the turbulent wind, perhaps) at wave- 
numbers outside the circular sector AOB . Under these conditions Phillips has 
shown that there is a critical angle to the mean wind which defines a locus in 
k -space and a narrow band associated with it (the shaded part) in which resonance 
will strongly excite and continue to feed energy into surface waves with wave- 
numbers within the band. Note that beyond the arc AB the band has been left 
unshaded. This represents waves with wave numbers which could be excited but 
are not since the surface pressures necessary to excite them are lacking. You 
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can see from Figure G-6 how it may be possible that a strong wind may not 
generate waves by this mechanism. If the' wind is uniform and lacking small 
structure itself, the radius OA of the sector AOB shrinks and, if O A becomes 
less than OD it falls below the resonance locus. Should that happen no small 
waves will be created by this mechanism. The mechanism is probably not unique 
but so far it is the only one which has been worked out that accounts for small 
waves and it is an excellent start. Theory could tell us many other things about 
waves with lengths 0 (1 cm> to 0 (10 cm) but I doubt that they would be any more 
useful in the real world than what weTiave already. 



Figure C-6. The Locus of Resonantly Excited Waves in Wave 

Number Space 
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The real ocean can be a very confusing place and, lacking good data, 
perhaps the best I can do is to describe a few of the things I, myself, have seen. 
One thing 1 can hardly overstress in the patchiness and evanescence of the fine 
structure. When a wind begins to generate waves on a calm surface you often 
find the rippled surface in cat's paws separated by regions of still calm water. 

The cat's paws may be small, 10, to 12 wave lengths, or relatively large, 500 to 
1000 wave lengths. They may persist, spread, and coalesce until the entire 
surface is rippled. However, if the wind is puffy and dies momentarily, the cat's 
paws may disappear again in a few seconds. If the wind sets in steadily and 
freshens, the longer gravity waves soon appear and the fine structure, while 
still there and wide spread, is usually very different in character on the different 
parts of the longer waves. Often the down-wind faces of the large waves are 
covered with much more prominent fine relief than are the up -wind faces. 

Sometimes the small waves on the down wind faces are moving at almost 
the same speeds as the larger waves and seem to hang fixed on them. I remember 
watching something like this in the Miles River for an hour or so. The wind was 
blowing about 30 kts over a fetch of perhaps 5 NM. The large waves were around 
3 ft high with characteristic lengths near 30 ft. Over and over again a packet of 
5 to 10 very steep small waves with lengths around 5 to 10 cm would suddenly 
pop up on the forward faces of the large waves. They ran just a little slower 
than the large wave so that ultimately they were overtaken by the crest and they 
disappeared as the crest went by. Before this happened, however, the small 
waves seemed to catch each other up; steepening on their forward faces until they 
looked like the overlapping tiles of a tile roof. At that time their forward faces 
were a centimeter or two high. The entire process from formation to vertical 
forward faces to obliteration when overrun by the crest occupied no more than 

10 seconds. The phenomenon w as just as local spatially occupying an area of 

, 2 

perhaps 10 m . 
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One more situation may perhaps be of interest as illustrating the need 
for fine wind structure if fine wave structure is to be .excited. My pier is located 
on an arm of South River. The water runs north-south and is perhaps a quarter 
of a mile wide. A twenty -five knot wind had been blowing from the south for 
some hours with an unobstructed fetch of 3/4 mi. The usual slope of gravity 
waves with periods around a second was running. What struck me when I looked 
at it was that it was abnormally smooth. For a wind that strong there were very 
few small waves. As I looked up to windward I saw a pier that had been toppled 
by the winter ice and was lying half submerged across the wind. In its lee the 
water was entirely covered with a plume of rhomboidal fine structure which 
extended some 30 ft to leeward where it suddenly disappeared. Looking further, 
every obstruction to the wind that I could see had its own pendant plume. I feel 
sure that this was an instance of the OA radius in Figure C-6 being less than OD. 
As the wind tripped over the obstructions small scale turbulence was generated 
in it and OA momentarily became greater than OD which permitted the small 
off-wind waves to be created. However, energy dissipation in the air flow 
rapidly sapped the energy from the small scales and, once it was gone, the 
small waves no longer had an energy source. With nothing to sustain them, 
they were rapidly destroyed by viscosity. I cannot prove that the story I have 
told is true but it certainly provides a rational explanation of what I saw: no fine 
structure on the open water even with a strong wind and plumes of rhomboidal 
fine structure behind each obstacle to the airflow. 

Considering the spottiness, rapidity of change, and variety of form of the 
fine structure it is interesting to turn to the results presented in the body of 
this report in Figures 5-1 to 5-3. Figure 5-1 plot average o° against wave 
height while Figure 5-2 plots it against wave period. Both the wave height and 
the wave period are visual obserations which notoriously fail to take proper 
account of the small waves. It is well known that if you have a record of water 
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elevation made at the same time and place a visual observation is made, you can 
usually duplicate the visual height observation very closely by reading all the 
heights on the record, putting them in order of size, and then taking the mean of 
only the highest one -third. The period estimate also ignores the small waves. 
But' the practically important waves for radar, and consequently for o , are the 
small ones. Thus one can say a priori that no connection is to be expected. This 
is precisely what the plotted points in Figures 5-1 and 5-2 confirm. 

Figure 5-3 plots average o against wind speed and this is something else. 
The wind speed, whether estimated or computed as a mean from cup anemometer 
data, is just as "gross" a parameter as the wave height or wave period. It 
tells us nothing about the turbulent structure of the airflow. However, the higher 
the wind speed the higher the Reynolds number and the higher the Reynolds 
number the more turbulent the airflow is likely to be. There is at least some 
reason to anticipate a connection between small waves and wind speed since the 
higher wind speeds are more likely to contain appreciable energy in the small 
eddies necessary to excite small waves. It isn‘t a tight connection but something 
ought to be there. Figure 5-3 suggests that there is and further, that you aren't 
going to be able to say anything about fine wave structure very precisely if all 
you know is the mean wind. 

The problems of the hydrodynamics of the fine structure of the sea and 
how to characterize it usefully urgently posed by radar altimetry over the ocean 
needs and deserves close study. 
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APPENDIX D 


FLIGHT PLANS 

This section includes the flight plans for all 16 flights. These plans were 
formulated before the scheduled flights. Comparing the flight plans with the 
flight summaries (Section 3.2) will indicate where and for what reason(s) the 
flight plans could not be followed. 
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Flight 1 - Nov-24-1969 
Low Altitude Shakedown 


1. Pre-Flight Briefing 

2. Take-Off 3:00 P.M. 

3. Fly to Test Area (Coastal Waters Near Wallops) 

4. Fly Stilwell Pattern at 6000ft. 

5. Fly Rectangular pattern at 6000ft. with the following parameters 

Pulse Width = 100 nsec 
Pulses/frame = 50 
Polarization = cross 
Peak Power = 12 watts 
Sweep speed = 50 nsec/cm 
Attentuation = 0 dB 
No. of Runs = 4 

6. Fly Stilwell pattern at 6000 ft. 

7. Return to Wallops 

Flight 2 - Dec-11-1969 

A) High Altitude Shakedown 

1. Briefing (Pre- Flight) 

a. Desire course rather than track 

b. Define Sea Direction 

c. Reference names 

d. Required data 

2. Take-Off and fly to test area 

3. Fly Stilwell pattern at 100 ft. 

4. Fly flight Pattern 2 (6 runs) 

B) Radar Parameters 

Pulse width = 100 nsec; runs 1, 2, 3 
20 nsec; runs 4, 5, 6 
Pulses/Frame = 50 
Polarization = cross 
Peak Power = 12 watts 
A/C Attitude = 0° 

Sweep speed = 50 nsec/cm 
Frame rate = l/sec 
Attentuation = 0 dB 
No. of Runs = 6 
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1. Fly Stilwell at max. altitude 

2. Fly Stilwell at 1000 ft. 

Flight 3 - Dec-12-1969 

A) 2nd High Altitude Shakedown 

1. Pre- Flight Briefing 

2. Take-Off and fly to test area 

3. Fly Stilwell pattern at 1000 ft. 

4. Fly to max. altitude 18K ft. 

5. Fly flight pattern 2 {6 runs) 

B) Radar Parameters 

Pulse width = 100 nsec 
Pulses/Frame = 50 

Polarization = cross (3 runs); Direct {3 runs) 

Peak power = 12 watts 
A/ C Attitude = 0° 

Sweep speed = 100 nsec/cm 
Frame rate=l/sec 
No. of runs = 6 

6. Fly Stilwell pattern at max. altitude 

7. Fly Stilwell pattern at 1000 ft. 

8. Return to Wallops 

Flight 4 - Dec-15-1969 

1. Briefing (9:00 AM) 

Discuss: 

a. Course vs. track (we desire course) 

b. Definition of Sea Direction 

c. Reference name of those on flight 

d. Inform crew of data required during flight 

2. Take-Off (10:00) 

3. Calibrate GYRO Compass of both A/C 

(Cessna and the DC-4) by flying together to the ship 
(Ship will be at Lat. 37 07' , Long. 73° 38 ie. Section 9, 
Quad. D. ) 

4. Pilot contact ship and request to start recording data 
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5. Start Stilwell pattern at ship with first heading in same 
direction as sea. DC- 4 at 1000 ft. , Cessna flying same 
pattern at 700 ft. 

6. Both A/C repeat Stilwell pattern; DC-4 at 5000 ft. , Cessna 
at 700 ft. 

7. Both A/C fly flight pattern 1 (12 runs) DC-4 5000 ft. , Cessna 
at 700 ft. 

A) Radar Parameters 

Pulse width = 100 nsec 
Pulses/Frame =1; for 1st 5 frames 
50; rest of run 
Polarization = cross 
Peak Power = 12 watts 
A/ C attitude = 0° 

Sweep Speed = 50 nsec/ cm 
F rame rate =1 /sec 
Attenuation - 0 dB 
No. of runs = 12 

8. Repeat step 6 (Stilwell pattern). 

9. Return to Wallops 

3.2.5 Flight 5 - Dec-17-1969 

1. Briefing (8:30) 

2. Take-Off 9:30 A. M. 

3. Fly to Wallops sec 9-D; Lat. 37° 07 T ; Long. 73° 38' 

4. Pilot contact ship at earliest time and ask for weather 
data (sea direction etc. ). 

5. At test area; DC- 4 fly Stilwell at 1500 ft; Cessna fly 
profilometer pattern at 300 ft. (single pulse radar runs) 

6. DC-4 fly to 10, 000 ft. and fly pattern 2 (6 runs) with 
following parameters: 

Pulse width = 20 nsec 
Puls e/ F rame =50 
Polarization = direct 
Sweep speed = 50nsec/cm 
Frame rate =l/sec 
Attenuation = 0 dB 
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7. Fly pattern 2 again (6 runs): Same parameters as #6 . 

except: Polarization = cross 

8. Fly to 15000 ft. and fly pattern 2 (6 runs) with same 
parameters as #6 

9. Fly Stilwell at 15000 ft. (or below ceiling): Also fly 
radar pattern #2 with parameters of #7 

10. DC-4 fly to -1500 ft. and fly Stilwell pattern. Cessna 
fly profilometer pattern at 300 ft. 

11. Return to Wallops. 

Flight 6 - Jan-5-1970 

1. Pre-flight briefing 8:30 A. M. 

2. Take-off 9:30 A.M. 

3. Fly to Wallops Sec. 9-D; Lat. 37° 07* ; Long. 73° 38' 

4. Pilot contact ship at earliest time and request weather 
data. 

5. At target area: DC-4 fly Stilwell pattern at 1500 ft. 
Cessna fly profilometer at 300 ft. 

6. DC -4 fly to 15000 ft. and fly pattern 2 (6 runs) with 
following parameters: 

Pulse width =100 nsec 
• Pulse/frame =50 
Polarization = direct 
Sweep speed = 50 ns/ cm 
Frame rate = l/sec 
Attenuation = 0 dB 

7. Fly pattern 2 again (6 runs); same parameters as #6 
except: Polarization = cross 

8. Fly to 20, 000 ft. and fly pattern 2 (6 runs) with same 
parameters as #6. 

9. Fiy Stilwell at 20, 000 ft. (or below ceiling): also fly 

radar pattern #2 with parameters of #7. 

10. DC- 4 fly to 1500 ft. and fly Stilwell pattern. Cessna fly 
profilometer pattern at 300 ft. 

11. Return to Wallops 

Flight 7. Jan-,6-1970 

1. Pre-flight briefing 8:30 A.M. 

2. Take-Off 9:30 A.M. 

3. Fly to Wallops Sec. 9-D; Lat. 37° 07* , Long. 73° 38' 



4. Pilot contact ship at earliest time and request weather 
data. 

5. At test area fly Stilwell pattern at 1500 ft. 

6. Fly to 15, 000 ft. and fly pattern 2 (6 runs) with following 
parameters: 

Pulse width = 20 nsec 
Pulse/Frame = 50, 100, 500, 1000* 

Polarization = Direct 
Sweep speed = B0 nsec/cm 
Frame rate = l/sec 
Attenuation = 0 dB 

* 4 frames for each pulse/frame 

7. Fly Stilwell at 15, 000 ft. (if weather permits) 

8. Fly Stilwell pattern at 1500 ft. 

9. Return to Wallops 

10. Fly over hanger at 5000 ft. GSE will be ready for 
. calibration. 

11. Land 

3.2.8 Flight 8 - Jan- 8-1970 

1. Pre-flight briefing 8:30 A. M. 

2. Take-Off 9:30 A.M. 

3. Make 3 passes overy runway (parallel with runway) at 
5000 ft. ; make another 3 passes over runway at 10, 000 ft. 
All 6 passes at minimum speed. Calibration using GSE 
will be made during passes. 

4. Fly to test area, Wallops Sec. 9-D; Lat. 37° 07' , Long. 

73° 38* 

5. Pilot contact ship at earliest time and request weather data. 

6. Fly Stilwell pattern at 1500 ft. 

7. Fly Stilwell pattern at 10, 000 ft. 

8. Fly pattern 2 (6 rims) with following parameters (10, 000 f 

Pulse width = 20 nsec 
Pulse/Frame =1 
Polarization = direct 
Sweep speed = 50 nsec/ cm 

Attenuation = 0 dB 

9. Repeat Step 8 with Pulse/Frame = 2 

10. Repeat Step 8 with Pulse/Frame =50 

11. Repeat Step 8 with Pulse/Frame = 500 
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12. Fly Stilwell pattern at 1500 ft. 

13. Return to Wallops £ 

14. Repeat step 3 

15. Land 

* Pilot notify Wallops plot when starting back to Wallops 

Flight 9 - Jan-9-1970 

1. Pre- flight briefing 

2. Take-Off 9:30 A.M. Q 

3. Fly to Wallops Sec. 9-D; Lat. 37° 07' ; Long. 73 38‘ 

4. Pilot contact ship at earliest time and request weather data 

5. At test area, fly Stilwell pattern at 1500 ft. 

6. Fly pattern 2 (6 runs) at 10, 000 ft. with following parameters: 

Pulse width = 20 nsec 
Pulse/F rame = 10 
Polarization = Direct 
Sweep speed = 50. nsec/cm 

Attenuation = 0 dB 

7. Fly pattern 2 again 
Pulse/Frame = 20 

8. Fly pattern 2 again 
Pulse/Frarae =148 

9. Fly Stilwell pattern at 10, 000 ft. 

10. Fly Stilwell pattern at 1500 ft. 

11. Notify ship that we are leaving the area. 

12. Return to Wallops 

Flight 10 - Jan-20-1970 

1. Pre-flight briefing 8:30 A.M. 

2. Take-Off 9:30 A.M. 

3. Fly to Tangier Sound; Lat. 38° 04* , Long. 75 00' 

4. Pilot request weather data from range recover ship 

5. Fly Stilwell pattern at 1500 ft. 

6. Fly pattern 3 (4 runs) at 5000 ft. with following parameters: 

Pulse width =10 nsec 
Pulse/Frame =1 
Polarization = Direct 
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Sweep speed = 50 nsec/ cm 

Attenuation = 20 dB 
(40 frames /run) 

7. Fly pattern 3 with: 

Pulse/Fr ame = 2 

8. Fly pattern 3 with: 
Pulse/Frame =10 

9. Fly pattern 3 with: 
Pulse/Frame =20 

10. Fly pattern 3 with: 
Pulse/Frame =50 

11. Fly pattern 3 with: 
Pulse/Frame =148 

12. Fly pattern 3 with: 
Pulse/Frame =278 

13. Fly Stilwell pattern at 5000 ft. 

14. Fly Stilwell pattern at 1500 ft. 

15. Return to Wallops 


Flight 11 


21 - 19^0 


1. Pre-flight briefing 8:30 A. M. 

2. Take-Off 9:30 A. M. 

3. Fly to Tangier Sound, Lat. 38° 04'; Long. 75 00' 

4. Pilot request weather data from ship 

5. Fly Stilwell pattern at 1500 ft. 

6. Fly pattern 3 (4 runs) at 10, 000 ft. with following parameters: 


Pulse/width = 20 nsec 
Pulse/Frame =1 
Polarization = direct 
Sweep speed = 50 ns/cm- 

Attenuation = 0 dB 
(40 frames/run) 


7. Fly pattern 3 with: 
Pulse/Fr ame = 2 
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8. Fly pattern 3 with: 

Pulse/Frame = 10 

9. Fly pattern 3 with: 

Pulse/Frame = 50 

10. Fly pattern 3 with: 

Pulse/Frame =147 

11. Fly pattern 3 with: 

Pulse/Frame = 278 

12. Fly Stilwell pattern at 10, 000 ft. 

13. Fly Stilwell pattern at 1500 ft. 

14. Return to Wallops 

Flight 12 - Jan-22- 1970 

1. Pre- Flight 8:00 A. M. 

2. Take-Off 8:30 A. M. 

3. Fly to Long Island Sound (Middle Ground) 

4. Determine Sea and Wind direction 

5. Take 10 frames of transmitted pulse 

6. Fly six rectangular patterns with: 

Altitude =10, 000 ft. 

Pulse/width = 20 nsec 
Polarization = direct 
Sweep speed = 50 nsec/cm 

Attenuation = 0 dB 
Pulse/Frame =1 30 frames - 

- 2 30 frames - 

= 10 30 frames 

= 50 40 frames 

= 148 40 frames 

= 278 40 frames 

7. Take 10 frames of xmitted pulse 

8. Land at JFK 


Flight 13 - Jan-26-1970 


1. Pre-Flight briefing 8:30 A. M. 

2. Take-Off 9:30 A.M. 

3. Fly to Wallops Sec 9-Dj Lat. 37° 07* , Long. 73 38* 
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4. Contact Range Recoverer ship and request weather information 

5. Fly Stilwell pattern at 1500 ft. 

6 . Take 10 frames of transmitted pulse 

7. Fly six rectangular patterns with: 


Altitude = 10, 000 ft. 

Pulse Width = 20 nsec 
Sweep speed = 50 nsec/ cm 

Attenuation = 0 dB 

Peak Power = 3 watt 

r un # puls e/frame pol ar 

^ 1 direct 

2 10 direct 

3 50 direct 

4 147 direct 

5 278 direct 

5 50 cross 


8 . Fly Stilwell pattern at 10, 000 ft. 

9. Take 10 frames of transmitted pulse 

10. Fly Stilwell pattern at 1500 ft. 

11. Return to Wallops 

Flight 14 - Jan-27-1970 

1. Pre-Flight briefing 8:30 A.M. 

2. Take-Off 9:30 A.M. 

3. Fly to Wallops Sec 9-D 

4. Take 10 frames of transmitted pulse 

5 . At test area fly Stilwell pattern at .1500 ft. 

6 . Fly six rectangular .patterns with following radar parameters: 

Pulse Width = 20 nsec 
Peak power = 12 watts 
Altitude = 10, 000 ft. 

Attenuation = 0 dB 

Sweep speed = 50 nsec / cm 


Pattern # 

Pulse/Frame 

Polarization 

Roll Angle 

1 

1 

direct 

0 ° 

2 

10 

direct 

0 ° 

3 

50 

direct 

> O 
O O 

4 

148 

direct 

u 

_ o 

5 

50 

cross 

0 

6 

50 

direct 

3° 
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7. Fly Stilwell pattern at 10, 000 ft. 

8. Fly Stilwell pattern at 1, 500 ft. 

9. Take 10 frames of transmitted pulse 

10. Return to Wallops (fly over Tangier Sound to check on 
ice cover). 

Flight 15 - Jan-28-1970 

1. Pre-Flight Briefing 8:30 A. M. 

2. Take-Off 9:30 A.M. 

3. Fly to Wallops Sec 9-D 

4. Take 10 frames of transmitted pulse 

5. Fly Stilwell pattern at 1500 ft. 

6. Fly six rectangular patterns at 10, 000 ft with: 

Pulse width = 20 nsec 
Peak power = 12 watts 
Attenuation = 0 dB 
Sweep speed = 50 nsec/ cm 

Polarization = direct 

Pattern # Pulse/Frame Roll Angle 



1 

1 

0° 


2 

10 

0° 


3 

50 

0° 


4 

148 

0° 


5 

278 

0° 


6 

50 

5° 

7. 

Take 10 frames of transmitted pulse 


8. 

Fly Stilwell pattern at 10, 000 ft. 


9. 

10. 

Fly Stilwell pattern 
Return to Wallops 

at -1500 ft. 



Flight 16 - Jan- 29- 1970 

1. Pre-Flight briefing 8:30 A. M. 

2. Take-Off 9:30 A.M. 

3. Fly to Wallops Sec 9-D; 37° 07' , 73 38' 

4. Take 10 frames of transmitted pulse 

5. Contact range recover for weather data 

6. At test area, fly Stilwell at 1500. ft. 
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7. Fly six rectangular patterns at 10, 000 ft. with: 

Pulse width = 20 nsec 
Peak power = 12 watts 
Attenuation = 0 dB 
Sweep speed = 50 nsec/cm 

Polarization = direct 


Pattern # 

Pulse/Frame 

Film 

1 

1 

B & W 

2 

10 

B & W 

3 

50 

B & W 

4 

148 

B Sc W 

5 

50 

Color 

6 

148 

Color 


8. Take 10 frames of transmitted pulse 

9. Fly Stilwell pattern at 10, 000 ft 

10. Fly Stilwell pattern at 1500 ft. 

11. Contact Range Recoverer when leaving area 

12. Return to Wallops 
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APPENDIX E 

SHARP LEADING EDGE OF MULTIPLE PULSE RETURN 
FROM OCEAN SCATTERING 


E-l. Introduction 

It is interesting to note that the leading edge of the multiple return pulse 
signal is more sharply defined then the top of these pulses (see Figure E-l). 
There are several reasons why the observer should expect this to be true. The 
variance of the signal strength at any time during the pulse is proportional to the 
amplitude at that point so that lower signal values on the leading edge results in 
less signal variance. In addition the density variance of leading edge traces is 
measured in a direction perpendicular to that leading edge; while the variance 
of the pulse at the top of the pulse is measured perpendicular to the time base. 
From Figure E-2 we see that when the pulse has a slope to the leading edge this 
results in the variance of trace density being reduced by a factor of the cosine 
of the slope angle; therefore, the full variance appears on the top of the pulse 
while the variance on the leading edge in trace density is reduced by two factors, 
the slope of the leading edge and by the lesser expected amplitude of the signal 
on tlie leading edge. 

A formula is derived for the variance of a square law detected Rayleigh 
fading signal corrupted by additive Gaussian noise. In the special case of the 
ramp -like waveforms that arise in satellite altimetry, the formula for the time 
dependence of the variance of the detector output is developed and sketched for 
comparison with flight test data. 
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Figure E-2. Expected Waveform and Variance 

E-2. Analysis 

A narrowband Gaussian (i. e., Rayleigh fading) signal s(t) with power profile 

E[s 2 (t)] = P(t) (1) 

is corrupted by additive zero mean stationary Gaussian noise n(t) of 
variance N. The combination y(t) = s(t) + n(t) is subjected to square law 

2 

detection. We seek the variance of the square law detector output v(t) = y (t), 
namely 

cr 2 (v(t)) = E[v 2 (t),] - E 2 [v{t)] (2) 

where E denotes the statistical expectation operator. 



From Eqs. (9) and (A-2) of Reference 1 we have, respectively. 


E [v(t) ] 

= P{t) + N 

(3) 

E[v 2 (t)] 

= 3[F(t) + N] 2 - 

(4) 


(In deriving (4) we have made use of the fact that p (t, t) = p (0) = 1. ) 

s n 

It follows that 


0 2 (v(t)) = 2[P(t) + N] 2 


(5) 


Observe that the ratio of the standard deviation of v(t) to the mean of 
v(t) is 


q(y(fc)) 

E[v(t)] 


(6). 


independent of t. Accordingly, the RMS value of the fluctuation of the 
signal is proportional^ to the expected value of the signal at each instant. 

E-3. Results Applied to Satellite Altimetry 

In satellite altimetry one receives a signal whose power profile 
is a ramp, 

P(t) = -Pt/T , 0<t<T . 

when a square pulse is transmitted and a wideband receiver is used. The 
standard deviation of the square law detector output therefore grows 
linearly with time during the interval [0, T], Accordingly, under low noise 
conditions (P/N »l).we expect the square law detector output pulses to 
have the general appearance sketched in Figure E-3. 


The constant in (6) is \/2~rather than 1 because in Reference 1 we assumed for 
simplicity that the signal was comprised solely of an in-phase component. 
Actucally, it has both in-phase and quadrature components, each with peak 
power P/2, and similarly for the noise. When this is taken into account, 
the *J~2 in (6) reduces to 1. 
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